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LONG-RANGE ELECTROSTATIC AND ELECTROMAGNETIC 
INTERACTIONS BETWEEN ATOMS 

SUMMARY 


The object of this calculation is to examine the contribution to the 
interaction energy between two neutral hydrogenic atoms in their ground 
vState because of their electrostatic interaction and the presence of an electro 
magnetic radiation field in its "vacuum" state. The separation distance R 
between the atoms is large compared to atomic dimensions, and is of the 
order of the transition wavelength X associated with the is-2p atomic 
transitions. Thus, the atoms are sufficiently separated so that the over- 
lapping of their respective charge distributions is neglected along with "spin. 
The problem is then handled using nonrelativistic quantum electrodynamics. 

The Hamiltonian of the system consists of the following parts: , 

which corresponds to the sum of the Hamiltonians of the respective isolated 

atoms plus the "vacuum" photon state of the radiation field; H^^\ which 
corresponds to the A(p) • P(p) type interactions between the atoms and the 

_ — ( 2 ) 

radiation field through the electromagnetic potential A (p) ; H , which 
corresponds to the A(p) . A(p) type interactions between the atoms and the 
radiation field; and H , the electrostatic interaction between the atoms. 

q 

Stationary state perturbation theory is then used to obtain the interaction 
energy of the system to fourth order in the electron charge e. By retaining 
hi^er order multipoles (up to octupole orders) in the expansion of H and 

/ i /c • r ^ ^ 

expanding the retardation factor ye y in A(p) in a power series, the 
interaction energy is obtained, accurate to quadrupole-quadrupole orders. 


The results are given as corrections to the various electrostatic 
multipoles corresponding to the dipole -dipole, dipole-quadrupole and 
quadrupole-quadrupole interactions between the atoms. The inclusion of 
the radiation field in the Hamiltonian of the system is found to give rise to 
retardation effects in the interaction energy which are a function of the 
separation distance R. The dipole-dipole results are in agreement with the 
calculations of Casimir and Polder [Ij and others. The dipole-quadrupole 
approximations can be written as a sum of three quantities proportional to 
R“^, R~®, and R~®. The R~® term corresponds to the purely electrostatic 
interaction and the others result from the inclusion of the radiation field in 


the system. In the limit of large R (R»A) , the dipole-quadrupole electro- 
static interaction energy (proportional to R“®) is found to be reenforced by 
a factor proportional to R~® and diminished by a factor proportional to R”'^. 
The results go over into the electrostatic case for small R (R«A) , showing 
that the retardation effects are unimportant for small separations as in the 
dipole-dipole case. The corrections to the quadrupole-quadrupole electro- 
static interaction (proportional to R”^®) are much more complex. The 

resulting expressions consist of terms proportional to R~^ R”^^ 

The above results are expressed in terms of integral functions over 
b = kR, where k is the magnitude of the wave vector k associated with 
the radiation field. 


INTRODUCTION 


The object of this calculation is to examine the contribution to the 
interaction energy between two neutral hydrogenic atoms in their ground 
state because of their electrostatic interaction and the presence of an 
electromagnetic radiation field in the system. The inclusion of the radiation 
field in the Hamiltonian of the system gives rise to retardation effects in the 
interaction energy [ 1 ] which are a function of the separation R between the 
atoms. These effects are unimportant for distances smaller than the wave- 
length X, associated with the radiation being considered, and increase in 
importance as R approaches X, In the case of atomic systems, the smallest 
wavelength associated with atomic transitions is large when compared to the 
atomic dimensions (\ » a.Q, where ag is the Bohr radius); hence, the 
problem can be treated using nonrelativistic quantum electrodynamics. This 
allows us to neglect the particle "spin" and to use simple product eigen- 
functions for the atomic and photon states. In addition, the above simplifica- 
tions make it possible to treat the problem using standard nondegenerate 
perturbation theory. This feature is crucial since the calculation must be 
pushed to fourth order before obtaining non-zero corrections to the interaction 
energy due to the radiation field. 

This problem was first treated to first order in both the electrostatic 
interaction potential and the radiation field (dipole-dipole approximation) by 
Casimir and Polder [1] . These authors used an elegant by asymmetrical 
method consisting of first treating the interaction between one atom with 
the radiation field and then considering this distorted field and the second 
atom as the initial configuration for subsequent approximations. Leech [2] 
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attempted to solve the problem (to first order) in a more systematic manner 
using nondegenerate perturbation theory. His results were in disagreement 
with Casimir and Polder, and, later were reported to be in error by Aub, 
Power and Zienau [3]. Subsequent calculations by Power and Zienau [4] 
using a different method verified Casimir and Polder's results. They applied 
straightforward perturbation theory to a reduced interaction Hamiltonian 
consisting of only the transverse component of the electric field vector. In 
this reduced Hamiltonian the electrostatic interaction is not given explicitly; 
hence if the interaction between the electrostatic interaction potential and the 
electromagnetic potential is desired, this method cannot be used. Subsequent 
calculations [5] using field theoretic techniques have since established the 
correctness of Casimir and Polder’s results. 

In this calculation the Casimir and Polder results are verified, using 
a systematic application of nondegenerate perturbation theory and a conven- 
tional expression for the interaction Hamiltonian. The calculations are 
extended beyond the dipole approximations of Casimir and Polder, by including 
terms up to octupole order in the electrostatic potential, in order to obtain 
interaction energies accurate [6] to quadrupole-quadrupole orders, and by 
retaining the first five terms in a power series expansion of the retardation 
factor of the electromagnetic potential. In this way, corrections due to the 
radiation field are obtained for each of the dipole -dipole, dipole-quadrupole, 
and quadrupole-quadrupole electrostatic interaction energies. In addition, a 
systematic analysis is performed to determine how the electrostatic inter- 
action potential is coupled with the electromagnetic potential to produce the 
retardation effects to the electrostatic interaction energy of the system. 


FORIVIULATION OF THE HAMILTONIAN 


Consider a system of two one-electron neutral atoms separated by a 
large distance R, (R » a^, whex'e ao is the Bohr radius) , interacting 
through their electrostatic potential H^ in the presence of a radiation field 

described by an electromagnetic potential A(p). In addition, assume the 
respective nuclei to be at rest (Born-Oppenheimer Approximation) and 
neglect all "spin" interactions between particles. The Hamiltonian for this 
system is then given by 


H = H, + H^^ + H + H 
I II q r 


( 1 ) 


3 


where and are the Haniiltonians for atoms I and n, given by 



H is the electrostatic potential [6] between the two charge distributions 
Q. 

defined in a coordinate system such that R is along the Z axis (Fig. 1) 
and given by 


y y (-l)^r^r^ 


M=t-L 


X E 


M"'' 

(01,0i)Y 
^ ^ 


4 7T 6 (Li-[ + ! 

[(2Li+ 1)(2L2+ 1)] ^2 

(02> ^ 2 ) 


-M* 


Zj _ 

=-L^ |_(Li + M) ! (Li - M) (L 2 - M) ! (Lz + M) I ' 


(3) 


is the electromagnetic field Hamiltonian given by 


2 

“r = E E 

«><0 X =1 


(4) 


The creation and destruction operators (? ^(/c) and Q (iO are defined 

A A 

through the quantized vector potential^ A(p) given by 


A (7) 



2 


E 

\=1 



iK • p 
e 




P 



1. See Heitler’s Quotum Theory of Radiation for the definition of the various 
quantities in A(p). 
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Expanding equation (2), one obtains 


H 


I 





Pi+Pj 


0 


2 ^ 0 “ 


( 6 ) 


where the first term in braces is just the unperturbed Hamiltonian for a 
single atom. If one defines the terms in equations (3) and (6) as follows: 




e 


2 




(7) 

( 8 ) 

(9) 


3^ = H 


( 10 ) 


then the total Hamiltonian given by equation (1) may be written as 


H = H<0> . + e 




II q / r 


Further simplification is accomplished by considering the fact that the 
radiation field Hamiltonian, as given by equation (4) , has as eigenvalues 
the number of photons present in any given state. Thus, if one picks the 
vacuum state defined by | , 0, 0, . . .> as the ground state of the radiation 
field, there are no photons interacting with the atoms in the ground state, and 

H^ may be written as H^*^^ and combined with H^^^ inequation (7). 

Incorporating this modification in the above expression, the total Hamiltonian 
may be defined as 
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If one further redefines the preceding quantities as follows. 


. (h<^> . ) , 


( 12 ) 

(13) 

(14) 


one may use perturbation theory techniques to solve the problem by considering 
the total Hamiltonian as a series expansion in powers of the electron charge £ 
given by 


H = 




1 „( 1 ) 


+ e H 


e^ 


(15) 


EIGENFUNCTIONS AND EIGENVALUES 
System Eigenfunctions 

A systematic study of the interaction between two atoms requires the 
knowledge of their respective states given in terms of eigenfunctions of the 
system. When perturbation theory is used to treat a given problem, the 
unperturbed state eigenfunctions must be determined completely^. 

In this problem, it will be assumed that the atoms of the system are 
sufficiently separated so that the overlapping of their respective charge distri- 
butions can be neglected. In this case, the unperturbed nth state of the system 
can be expressed as products of eigenfunctions given by 


(1) ip^^^ <U) ip^^^ (r) 


(16) 


2. This state is sometimes referred to as the ground state of the system. 
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If one examines the system such that the unperturbed state corresponds to the 
ground state, then equation (16) is replaced by 




(17) 


where (I) represents the unperturbed eigenfunction for atom I in its 

lowest state, and (r) represents the state of the radiation field in its 

ground state (sometimes referred to as the "vacuum" state) . Since hydrogenic 
eigenfunctions are used for the atomic states, and the electromagnetic potential 

is expressed in terms of creation and destruction operators (k) , (? (/c) , 

A A 

alternate definitions for the eigenfunctions given in equation (16) are necessary. 

Let us define as follows: 

n, (N) 






( 0 ) 

n, (N) 

( 0 ) ^ 
a 


, ... (k) , • 


... 


I(n) n(n) 


I(n,i , m) II(n,i , m) 


(18) 


(19) 


If one wishes to express an intermediate state in which atoms I and II are not 
necessarily in the same state, and the photon state contains more than one 
photon having different parameters k,X then one denotes this state by 


jp , = 






I(n) H(m)^ 


( 20 ) 


where atom I is found in the state having quantum numbers n; atom II has 
quantum number m; the photon state contains N photons having parameters 
x,A; and N' photons with parameters k',A'. 


The hydrogenic eigenfunctions are defined by [7] 


,p (I) = R 

n j j XXX 


n.f 


(T;)y^(0i,0i) 





where the R . (r) 
n,i 


and (6, (p) are also defined in Reference 7. 


(21) 


8 



The specific eigenfunctions ^ used in this calculation are those 

corresponding to the ground state for which n=i, i = 0, m = 0; 


^1,0,0 ^ 


..>!❖ ■ if (t) 


and for the 2p states for which n = 2, i = 1, m = 0, ± 1: 


'2,1,0 




1 / Z \ 2 / 2r \ “ 2an 


4^] 2n 


cos 6 , 


' 2 , 1,+1 


\ Y+\ 1 /zy^^/zr\ 

'> ^ W / \ ^» / 


sin e cos ( p , 


2 , 1-1 


" 2 , 17^1 


1 I Z 


W^TT 


V / V - — 

^ / Zr \ 2an 


e ” sine sin<^> . 


The above functions, together with the photon states, redefined as 


= .... \{K), .. 


and whose vacuum state is given by 




.. . 0 , 


completely define the system under consideration. 
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System Eigenvalues 


Having defined the eigenfunctions and their various forms, one proceeds 
to describe the eigenvalues corresponding to the above eigenfunctions. First, 
the hydrogenic eigenfunctions are solutions for the unperturbed state of the 
system. The photon states are eigenstates of the radiation field Hamiltonian 
H^ given by equation (4) ; that is. 


YjYj 

K A 




K A 


\(k), 


(28) 


In addition, these photon states satisfy the usual relations given below; 


a |, 

a 

Q a"*" 



= ^^(N+l) 


, ... (N+l), ... 




, ... (N-l), ... 


= ^^W) nT(N) 


, ... N, ... 


= \/(N+r) (N+i) 


... N, ... 


(29) 

(30) 

(31) 

(32) 


It is clear from equation (28) that the vacuum state is an eigenstate of H^ 

with eigenvalue N = 0; hence, one is justified in including H with the 
unperturbed portion of the Hamiltonian. 


Other relations between the photon eigenstates and the field operators 
which will be used later on are given below for future reference [8]: 


aA(«) 


= \jw (N+l) 


, N^(k), ... N^,(^')'^ 


(N+1)^(k), ... (N-l)j^,(K') 


(33) 
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= \l (N+l) « ^/N 


. N^(k), ... Nj^.U’), 


(N-1) (7). ... (N+l)' (/c') 


(34) 


fl^(K) ^31^1 ('f') 


, N^(k), ... Nj^,(K'), 


= ^y (N+l) ' (n+1) 


(N+1)^(k), ... (N+l)j[^,(K') 


(35) 


(S^(k) a^,(K') 


. N^(k), ... Nj^,(K'), 


= ^yN^ 


(N-I)^(k), ... (N-l)j^.(K') 


(36) 


The energy eigenvalue for the system may now be expressed in terms of the 
unperturbed eigenfunctions using nondegenerate perturbation theory reported in 
an earlier NASA publication [9]. Thus, the solution to the perturbed eigen- 
value problem, 


H4> = E 4- 

n n n 


(37) 


may be obtained by expressing H as a series expansion of the form 


„ „(0) „(1) 2„(2) 3„(4) 4„(4) 

H = H +eH +eH +eH +e H +... , 


(38) 


whose eigenvalue may be expressed as 


n n n n n n 


(39) 
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are given in terms 


where the energy corrections to ground state energy 


n 


of matrix elements of the perturbation Hamiltonians H ^ , i # 0. Since one 
desires to obtain results to the fourth order, one needs all the terms indicated 
in equation (39) . But since H as given by equation (38) contains only 

and one can obtain the terms corresponding to equation (39) 

(3) (4) 

directly from Reference 9, by letting H and equal to zero. Doing 

this, the resulting expressions are: 


n 


h(1) 

nn 


E 


( 2 ) 


n 


= H. 

nn 


(3) 


n 


= - H 


( 1 ) 


nn 


- S 

n¥n 


Z 

n¥n 


Z 






nn' n'n 


n'^n ^ 

\ nn' n'n nn' n'n / 


- z z 

n'i^n n"=i^:n 


( 1 ) (!) ( 1 ) 
xi g Jfl A * ^ xl * * 

nn' n'n" n"n 



(40) 


(41) 


(42) 
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=■ I 

n H 


nn’ n*n 


- I 

n H 


h“! hW 


(e^-E^) “" n^„ (Er-E‘?>)' 

■*' ^ / ( 0 ) ( 0 ) \® 

n¥n ( e' ^ - e'V ‘ 

y n n 


I 

H<‘;h‘J>„H<,V(2e“'>-E<?»-E<»>) 

(1) y y nn* n*n" n"n \ n n* n*' / 

““ n*?in n'¥n (e^®^ - E^®^ Ve^®^ - E^JJ^V 
y n n’ y y n n’* y 


-S E 

n*#n n'¥n 


( H^,V + H^,V + ) 

\ nn^ n"n nn* n"n nn^ n*n” n^*n / 

/'e‘'’'-eWWeW-eW\ 

y n n* / y n n" y 


E E E 


H<‘) „») „( 1 ) „( 1 ) 

nn* n*n" n"n"’ n***n 


n*#n n'¥n n'*'#n - E^[^^ \ (e^^^ - E^^^Ve^®^ - E^?,^, \ 

y n n* / y n n" yy n n J 


(43) 


The interaoUon energy E__, given In eqnaUon (39), may be evaluated to the 
desired fourth order using the above results. 


CALCULATION OF GENERAL MATR IX ELEMENTS 


Introduction 

Analysis of the terms in equations (40) through (43) shows that the 

corrections to the interaction energy E^^^ consist of various matrix elements 
of the form 


„(i) „(i) „(i) „(i) u a 9 

nn nn' n'n" n"n”' ’ 
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These terms need to be evaluated and then combined in order to obtain the 
interaction energy between the two atoms. The above operations are simplifed 
considerably if one neglects those terms which correspond to interactions 
between the radiation field and either one of the atoms (10] . Thus in sub- 
sequent discussions will refer to only the interaction energy between 

atoms I and II, either through electrostatic interactions or through the 
radiation field. In the course of evaluating the various corrections to the 
interaction energy, an analysis is made on all the terms which make up the 
overall interaction energy to fourth order. 


First-Order Corrections 


to the 
by 


One proceeds to evaluate the various matrix elements 
energy corrections. The first-order correction 


corresponding 
_ (i) . 

E IS given 
n 





(44) 


Using the unperturbed eigenfunctions given by equation (16), the definition 

for given by equation (13) , and since affect only atoms I and II, 

equation (44) simplifies to ’ 


E 


( 1 ) 


a 




H 


( 1 ) 


/ ( 0 ) , ( 0 ) , > 


+ < (r) 


H 


( 1 ) 

n 




(45) 


Before evaluating equation (45), one simplifies the expression for as 

given in equation (8) be recalling that in obtaining A{^) , the condition 

V • A(p) = 0 was used. (See References 10 and 11 for details. ) Since 

r t— ► >^"1 — ( 1 ) 

I^A, PJ = - iii V • A, one notes that A and P commute and may now 

be expressed as 
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'4 


^ i/C • Pi /“T * Pi 

(? (k) e ^ - a (k) e 

A A 


(K) • Pj(ri) 

r(l) 


(46) 


The specific position dependence in the above expression for Hj is 


included because the operators P and A refer to two different coordinate 
systems. Their relationship is best described by Figure 2 which also gives 
the following relations: 


Pi = Rj + rj , P2 = Rjj + rj 


(47) 


where Rj and R^j are related by (Fig. 1) 


R = Rjl-Kj 


(48) 


Analysis of terms in equation (46) shows that consists of products 

proportional to 


■ Pj(ri) 


a^(K)e^(/c) • Pj(ri) 


(49) 


Since the momentum operators act only on the atomic states, and the field 
operators (J , (/c) , a^(K) affect only photon states, one obtains 

A A 


a^(K) • Pi(ri) 




- €^( 7 ) n/N 

A 


♦Jn- 1) 


(I) 

^n 
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Using this result, the first term in equation (45) is proportional to 


(r) 

’^(N) ^ ^ 


''■(N-l) 





Pl(ri) 


(I) 


Terms of this type vanish, since one must satisfy the relations [8] 


N» 

, X'l , • • • 


a 


N, . . n/N 


and 


N' 

, , ... 


a 


N, = '^ «N',(N^1) 


between initial and final photon states. Hence, the terms in equation (45) 

vanish and the first order correction does not contribute; that is, 

n ’ 


E<‘> = = 0 

n nn 


(50) 


Second-Or(Jer Corrections 


The second-order correction to the interaction energy is given by 


a aa 


I 


h(1) „U) 


q:¥q! - E^?^'\ 

^ O' / 


(51) 


( 2 ) ( 1 ) ( 2 ) ( 2 ) 

Since E' contains both ' and , one first evaluates H' in 

a 

( 2 ^ {21 

terms of A(p). From equation (14), one notes that H' ’ consists of H ' \ 

( 2 ) ( 2 ) ( 2 ) — — ^ 

Hjj and . Since is given in terms of A (p) , one uses equation 

(5) to get 
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X 






^ /T, -i(«-'^')"Pi 

+ (J^(k) a^,(K') e ' ' 




(52) 


( 2 ) ( 2 ) 

Using the above definition for , and the definition for H given in 

equation (14), the first term in equation (51) becomes 


aa 



1 

2jj.c^ 


A^iPx) • Aj(^i) 



+ 



II 


1 

2^ 


Ajj("^ 2) ‘ Ajj(^) 





(53) 


Evaluation of the above terms is simplified by noting that the first two terms 
give the same result so that only two terms need to be evaluated. Rewriting 
the first term, one obtains 



(I) ^ 


( 0 ) 

(N) 


(r) 



A(pi) • a(px) 





(54) 


Since the operator A • A consists of products of the form (2^(5 , (J^ (2'*', 

(2(2"*^, (2(2, the above expression may be solved by applying the relations 
listed in equations (29) through (36). These relations indicate that only the 

term containing (2(2"*" gives results which are nonzero; that is, 
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N^(k) 




N^(k) 


= ^ (N+1)^(k) 


W K/c' 


(55) 


Using this expression in the matrix element, expression (54) yields 





(N+1)^(k) 


\,X' 


6 , e 

k,k' 


i(K 





Separating the sums over k,X and and summing over \,x' and 

gives 



/ 2irKc^ \y> / 2 \ 
\2^c2vo1./^ \ck / 

K 


(N+l) (K) 



(56) 


The above term is significant because it shows explicitly how the R 

dependence contained within the exponential e^^*^ k )• Pi vanishes when 
calculations are carried out to the second order, making the quantity in 
expression (56) independent of the interatomic separation. Terms having 
this property do not contribute to the interaction energy between atoms I and 
II and only give rise to self-interaction energies which are neglected in this 
calculation. Hence, the first two terms of equation (53) do not contribute. 
The last term of equation (53) is given by 



(I) 


( 0 ) 


n 


(II) 



(57) 


This term corresponds to the electrostatic interaction between the two atoms 
and, as such, contributes to the interaction energy. Since all other terms 
vanish, equation (53) is given by 


h ( 2 ) 

aa 




(58) 


Expressing in terms of and the numerator of the second 

term in equation (53) becomes 
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where a ' denotes all the quantum numbers for ttie intermediate states. 
Expanding the first term in equation (59) yields 


,(i) 


<i) 






,( 1 ) 




i)(r^ 

^n,n" 


(60) 


Referring to the expression for H. 


.+ 


1 


one notes that it consists of two terms 


involving (k) and (2 (k), which, when operating on the photon eigenstates, 
A. A 

as previously stated, require that the initial and intermediate photon states 

nates of atom II, a photon must be exchanged between the field and atom I to 
obtain nonzero results. This is illustrated below, using the following inter- 
action diagrams [10]. 


differ by one photon. Since 


is independent of the coordi- 



According to diagram (a) , atoms I and II are initially in the same 
state (single lines), and the photon state corresponds to the vacuum state. 

The interaction operator "carries” the system into an intermediate 

state where atom I is found in a different virtual (double line) state, 

atom n remains in its initial state, and the photon state contains one virtual 

photon [10]. The next interaction operator "brings" the system back 


to its original state. Diagram (b) illustrates the case in which the initial 
photon state is not the vacuum. If the initial system is such that both atoms 
are in the same state and the photon state corresponds to the vacuum, then 
only diagram (a) applies. Similar analysis on the fourth term of equation (59) 
yields identical results. 
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The second and third terms of equation (59) give similar results, as 
may be seen by considering the second term in its expanded form; that is, 


<a I I I I ^ | 




(1) .,.(0) 


JO) 




^ ^n,n” 

Vn' 


(61) 


To have nonvanishing results in equation (61) , the photon states must 
differ by one photon. Hence, there must be a photon exchanged between atoms 

I and n because both and appear in the above matrix product. 

The only possible nonzero interactions are illustrated below. 



Diagrams (a) and (b) correspond to 


correspond to 


( 1 ) 



H 


( 1 ) 


(1) 



n X ’ 


(c) and (d) 


Analysis of these diagrams shows that the 


conditions of equation (61) cannot be met since the atomic states change in the 

transition from initial to intermediate state during photon exchange. Hence, 

the second and third terms of equation (59) do not contribute either. Therefore, 

(2) (2) 

the second term of E ' is identically zero. Since the second term of E ' 

a a 

vanishes, only the term given by equation (58) contributes to the second-order 

correction; that is. 


E<^> 

a 




(62) 


Later, it will be shown that this term vanishes when the unperturbed states 
are s-states. 
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Third-Order Corrections 


The third-order correction ^ is given by equation (42) . Since 

has been shown to vanish, is given by 

nn a 


- Z 

a^+a. 


^ aa' a^a aa 

V O' a' 


a' a'a / 


(1) (1) (1) 

Y Y Q?g * a^a" a"a 

Q!¥q! a"^a ( e ^^^ - 

y a a' J y a a” 


Considering the factors in the first term, one notes that only one of them needs 
to be examined in detail, since the other factor is similar in form. Using 

equations (13) and (14) and expanding , one obtains 

aa'a'a 




"(N-)"' “I ri. 


Further expansion of the first term of equation (64) gives 

-<(“i «j" hi><“'| «r ”i“'| |“)> 

.<^| <^|«U 

The condition 6^ requires that atom n remain in its initial state 

during the interaction between atom I and the electromagnetic field. In 

addition, the operator couples only photon states whose number of 

^ (2) 

photons differs by unity, and ^ couples photon states whose number of 
photons differs by zero or two. Therefore, the only way to satisfy both 


requirements is to allow the initial state (r)> to contain one photon. 

(0) \ I (^1 / 

This is not allowed since must be the vacuum state; consequently, 

the above term cannot contribute to the interaction energy. 
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Referring to equation (64) , one notes that the preceding discussion 
also applies to the fifth term by simply interchanging atoms. Similar argu- 
ments are used to show that the second and fourth terms of equation (64) do 
not contribute either. Finally, the third and sixth terms of equation (64) , 

which contain matrix products of the operators and do not 

contribute because, as noted before, H 4 couple states whose number of 

(2) ^ 

photons differ by unity, and H' ' can only couple identical photon states. 

*1 12 ) ( 1 ) 

Applying the same arguments to H ' H ' , , one finds that the first term in 

=> Q/Q/t aiQ-’ 


equation (63) does not contribute to E 


(3) 


a 


Possible nonzero configurations 


corresponding to equation (64) are illustrated in Figure 3. These diagrams 
are introduced here only to illustrate the method used to analyze more complex 
situations later on. 



Figure 3. Possible nonzero interaction diagrams of equation (64). 

The last term of equation (63) consists of matrix products involving 
only. Using equation (13) to expand this term, one obtains 
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In the above equation, three types of products involving and are 

( 1 ) ( 1 ) ^ ( 1 ) 
found, products involving only or and products coupling H' 

( 1 ) 1 11 1 

and in various ways. If one requires that initial and final states be 

identical, one notes that all the above terms violate this requirement in one 
way or another. This can best be seen by analyzing the diagrams in Figure 4, 
which depict typical nonzero configurations of the terms in equation (65). 



Figure 4. Some interaction diagrams of equation (65). 
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1 


Since the last term is zero and the previous terms of equation (63) 
do not contribute either, one finds that the third-order correction to the 
interaction energy vanishes; that is. 


= 0 . (66) 

n 

Fourth-Order Corrections 

The fourth-order correction to the interaction energy is given by 
equation (43). Using the fact that as well as the products of the form 

given by equation (59) vanish, the expression in equation (43) becomes 



Expanding the first term of the above equation yields 


M, 

p{4) _ ^ aa' a'a 




Z Z 


Qtt* 0>'Qi" O' "a OQ;* Of"g gg* a*Qf'* a''Q! 


(eW - E<»n (e*") - E<?M 
y Of Of’ y y a q;" j 


(1) (1) (1) (1) 

■ II II ■ a a ■ a 

ctcc^ a’cx * a *Qf"' 


( 68 ) 


aL a (e “» ■- E '“>1 (E‘'>'-E<?,')fE<“»-E«>M 

y o a' /y o; Of" y y or a”' y 


The three types of terms corresponding to various interactions in the 
above equation are analyzed as follows: 


25 



a. Interactions between field and either of the atoms: The first and 
fifth terms in equation (68) correspond to this type interaction. These terms 
do not contribute to the interaction energy as may be seen by analysis of 
diagrams (a) and (d) of Figure 5. 


b. Interactions between field and both atoms: The second and fourth 
terms of equation (68) correspond to virtual photon exchange between atoms I 
and II and, as such, contribute to the interaction energy. (See diagrams (b) 

( 2 ) 

and (c) of Figure 5. ] The terms in equation (68) involving both H' and 

( 2 ) ^ 

H ' are also in this category; except that the interaction between atoms 

( 2 ) 

takes place through the operator H , and the field interacts with either 
( 2 ) 

of the atoms via H . Hence, these terms contribute only when matrix 

1, ii 
( 2 ) 

elements over H are nonzero. This group of terms is illustrated in 

^ ( 2 ) 

Figure 6. Note that in this case H gives rise to instantaneous emission 

1 , 11 

and absorption or absorption and emission of two virtual photons. This requires 
that the initial and intermediate atomic states be the same. Therefore these 

( 2 ) 

terms contribute only when the matrix elements of H' over initial states 

q 

are nonzero. / 



Figure 5. Interaction diagrams corresponding to first, second, 
fourth, and fifth terms of equation (68). 


c. Interactions between atoms I and II only: The only term of this type 

( 2 ) 

is given by the last term in equation (68) . Since H does not affect the 

q 

photon states, the initial, intermediate, and final photon states must be the 
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Figure 6. Interaction diagrams corresponding to the third, sixth, 
seventh, and eighth terms of equation (68). 


same. The diagram corresponding to this term is given in Figure 7. Note 
that, when electromagnetic interactions are neglected, and the interaction 
energy is computed to fourth order in the electron charge e, only this term 
and that given by equation (58) are obtained. This is the reason for singling 
out this term in Figure 7. 



Figure 7. Interaction diagram 
corresponding to the atom- atom 
coulomb interaction term given 
by the last element of 
equation (68). 


The second term of equation (67) 
is given by products of the form 


H 


( 2 ) „( 1 ) ..( 1 ) 


H 


H 


cva' a'a" a"a 


H 


(1) ..(2) ..(1) 


aa' 


H 


a' a' 


H 


a^'a 


and H 


(1) ..(1) 


aa' 


H 


a' a' 


H 


( 2 ) 


Using pre- 


(1) 


( 2 ) 


vious definitions for H'*' and H 
the expansion of the above products is 
quite lengthy. These expansions are 
needed to select the nonzero terms 
which contribute to the interaction 
energy. If the indicated expansions 

„ , , „( 2 ) „( 1 ) „( 1 ) 

are performed on only H ,H , ,,H ,, 

ao!’ a'a" a"a 


and the results are expressed in terms of interaction diagrams, the task is 
simplified considerably. For the other two terms, only the resulting expan- 
sions and their corresponding interaction diagrams will be shown. To obtain 
the interaction diagrams, one uses previous results, illustrated in Figures 3, 
4, 5, 6, and 7. 
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Expanding one of the terms yields the following: 



( 69 ) 


In associating the above terms with an interaction diagram, one finds 
that some terms can have more than one possible configuration which gives 
nonzero results. This occurs especially for those nonzero terms which violate 
the requirement that final and initial states be the same. In other cases, one 
of the possible combinations corresponds to the term that contributes to the 
interaction energy. These situations are illustrated in Figure 8, where (a) , 
(b), and (c) correspond to the first category, and (d) , (e), and (f) to the 
latter cases. 

Note that diagrams (a), (b) , and (c) of Figure 8 do not contribute 
to the interaction energy and only (d) of the remaining group contributes. 

The photon states corresponding to diagrams (a) , (b) , and (c) in Figure 8 
are: 
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Figure 8. Interaction diagram corresponding to various nonzero 
configurations of terms in equation (69) . 


(a) 

(b) 

(c) 


<(o| sIm I ^ < 0 1 a^(«) I | . 

h-'x'j a^(“) I I a^,(iT') I ^ , 


o| aj^(«) 3 .^,(. k ') k \ 




k'X' 


a^(«) I kx, 


, k'^<(^X, k'X' a^,(K') K^<^\ a^(Kj 0^ . 


As one can see, these photon states may be associated with various inter- 
mediate states in the diagrams. The ordering of these states may be obtained 
by reading the matrix elements from right to left as one follows a diagram 
from the bottom up. 


The interaction diagrams corresponding to the terms of equation (69) 
are listed in Figure 9; only one configuration for each term is shown, whether 
it contributes to the interaction energy or not. The listing of these diagrams 
is the same as the ordering of the terms in equation (69). Analysis of 
Figure 9 shows that only six diagrams contribute to the interaction energy: 
the fourth, fifth, and the last four. The remainder of the diagrams in 
Figure 9 correspond to either unacceptable situations or to self-interactions 
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Figure 9. Interaction diagrams corresponding to 
terms of equation (69). 

between either of the atoms and the radiation field. Of the six terms which 
contribute to the interaction energy, four involve the electrostatic interaction 
(2) 

operator in two different ways. This can be seen by referring to the 

last four diagrams of Figure 9. 
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The next term to be considered, when expanded, gives 


aa' 


a'a" a 


a 







(70) 


The corresponding interaction diagrams are listed in Figure 10 where the 
rules outlined previously are applied to this case also. 

Analysis of Figure 10 shows that six diagrams contribute to the 
interaction energy as before: the fourth, fifth, and the last four. The 
remainder are classified as in Figure 9. 
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Figure 10. Interaction diagrams corresponding to terms in equation (70). 
The last term in this group to be considered is given by 
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The corresponding set of interaction diagrams is given in Figure 11. Analysis 
of this figure shows that six more diagrams contribute to the interaction energy; 
the second, seventh, and the last four diagrams. 





Figure 11. Interaction diagrams corresponding to terms in equation (71). 


The importance of the preceding diagrams lies in the fact that one now 
finds nonzero terms which contain both the radiation field operators and the 

( 2 ) 

coulomb-type operator, H' \ coupled together within individual factors. 

Q 

The last four diagrams in Figures 9, 10, and 11, respectively, show the form 
of the interactions in a systematic and informative way. Later, it will be seen 
that these terms play a very important role in the evaluation of the interaction 
energy to various approximations. Instead of collecting the contributing terms 
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just obtained, one evaluates the last term of equation (67) and then collects 
all the nonzero terms which contribute to the interaction energy. 

Using the definitions for the last term in equation 67 expands 

into 




(72) 


The terms in equation (72) contain various types of factors, correspond- 
ing to interactions between each individual atom and the electromagnetic field 
or between the electromagnetic field with both atoms. In order to illustrate 
those terms we again pick one representative interaction diagram for each of 
the terms which does not, contribute to the interaction energy, but list the 
various combinations which do contribute. This is necessary since we have to 
consider all possible configurations as indicated by the sums over k,X; k',A' 

of the operator in the various intermediate states. Until now, only one 


interaction diagram per term has been sufficient to take care of the nonzero 
terms. This is no longer sufficient when listing the possible configurations 
for this term. This point is illustrated by considering a specific term of 
equation (72) and listing the possible configurations of this term in Figure 12. 
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Figure 12. Interaction diagrams of some of the possible configurations 
corresponding to the 13th term of equation (72) . (The photon 
parameters k,\; are included to show the variations in 

intermediate photon states. ) 


Analysis of Figure 12 shows that the first two diagrams are two 
possible nonzero configurations for the 13th term of equation (72), which 
contributes to the interaction energy. 


The interaction diagrams corresponding to equation (72) are listed in 
Figure 13; one includes the various configurations for the terms which con- 
tribute to the interaction energy, and only a representative diagram is 
included for those terms which do not. When a given term contains both zero 
and nonzero configurations, the configurations which contribute are used. In 
Figure 13, the various interaction diagrams are numbered; those which 
correspond to the same term of equation (72) are denoted by primes [e. g. , 
(4) and (4) ' correspond to the fourth term of equation (72) ] . Analysis of 
Figure 13 shows 12 diagrams corresponding to 6 different terms, in equation 
(72), which contribute to the interaction energy. This is by far the largest 

,(4) 


number of nonzero terms resulting from a single term in E 


a 


One can see 


that out of a large number of possible combinations indicated by the sums over 
Q!’, a", a"', a number of terms do not contribute, solely because of the 
restrictions on the various photon states. Collecting all the terms which 

,(4) 


contribute to the interaction energy correction E 


a 


one finally obtains 
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1 


= z 


,( 4 ) 

a a* 





+ C H. 


+ < H 






( 73 ) 
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Combining E 


(4) 


a. 


the interaction energy E 


with all the other nonzero corrections, one obtains 

to the fourth order; that is, 

a 


E = E^®^ + e^ E^^^ + e^ E^'^^ 
O' O' O' a 


(74) 


(2) (4) 

where E^ ' and E^ are given by equations (62) and (73), respectively. 

It is easy to see why one must consider fourth-order corrections to the 
interaction energy when considering both electrostatic and electromagnetic 

( 2 ) 

interactions between atoms. In cases where E vanishes, the only con- 


a 


tribution to the interaction energy comes from E 


(4) 


a. 


(4) 

In this case, E is 


a 


simplified somewhat because several terms in equation (73) vanish. 


INTERACTION ENERGY BETWEEN HYDROGENIC ATOMS 
IN THEIR GROUND STATE 

Interaction Energy of the System 

The unperturbed system in this case consists of atoms I and II initially 
in their ground state configuration (Is - Is); the electromagnetic field is in 

its vacuum state ^ I j • • . 0, • , ^ . From the previous discussion, the 

energy eigenvalue to the fourth order is given by 


E 

a 


a 


2 

+ e 


e‘"> 


4 

+ e 


e(4) 

a 


( 75 ) 


Since one is interested in the corrections to E^*^^ rather than E , one 
rewrites equation (75) as ^ 


AE 

a 



a a 


(76) 
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where AE is defined as the correction to the ground state energy due to the 
a 

interactions. Henceforth, this will be the only quantity of interest. To 

( 2 ) 

evaluate the terms in equation (76) , one considers first , given by 


E 


( 2 ) 


a 


H 


( 2 ) 


a 


(77) 


In terms of equation (3) , the above becomes 


0,0) n(i, 0,0)1 <^o. ... f f 

\ ' N Li=l Lg=l 

M=+-L 

< (I) Y (n) 

“U h2^ 


^ E r 1V-) 

[(Lj + M) ! (Lj - M) ! (Lz + M) ! (L2 - M) ij 


(-1) ^ r, ^ (47r)(Li + L,)! 

jjLj+La+l j" + 1) (2 Lj + 1) ] ^2 

'M 


17 1 ( 1 , 0 , 0 ) n(i,o.o) 


(78) 


( 2 ) 

Since the operator factors into quantities corresponding to atoms I and 

n and does not contain field operators, the above expression may be written 
as 


.(2) 


. . = ^ ^ ^ ' Y 

^ Lj L 2 M |^(2Li + 1) (2L2 + 1)] (Li + M)! (Lj - M)! (L 2 + M)! (L. - M)! 


(-1) ^ (4ir)(L, + L,)l 


X <1(1, 0,0) 


r^‘ ( 6101 ) 


1(1 


.o,oi^<^(i,o. 


0 ) 


rLay - M * 


U(1 


,o,oj^ 


Using the definition for the atomic eigenfunctions given in equation (18) , the 
matrix elements in equation (79) factor into products of the form 


1(1, 0,0) 

\ 

<(Ri,o(I) 


M 


1(1, 0,0)^ 

Rl,o(I)\ <Yo<>(I) 


M* 

(I) 

I-*! 


Yo“(I) 
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nil Mini 


Using the relations between spherical harmonics [ 12 ] , 



= (- 1 ) 


mz I ( 2 £, + l)(2iz+ 1 ) 
4tt (2ig+ 1) 


-1 1 


■/z 


C (i li 2^ ;nii,-m2, mg) C (i 1, iz,-? 3; 0, 0, 0) , 

(80) 


one obtains 







by setting mj = m3 


and mz = M. 


Substituting this, equation (79) reduces to 


,( 2 ) 


= E Z 

Li.Lz M 




L1+L2+I 


[(2L^ 


(81) 


+ 1) (2L„ + IJ (L, - M)!(L, + MM(L„ - M)!(L, +M)i 


Equation (81) shows that the only term in the series corresponds to 
the case where Lj = Lz = 0, and M = 0. This term is the monopole contri- 
bution of the electrostatic potential between the two atoms. When the charge 
distributions are neutral as is the case here, this term does not contribute. 

12 ) 

Hence, one sees that is identically zero for the case in which the atomic 

initial states are picked to be the (is) ground states; that is 


E 


(2) 

1 


= 0 


(82) 


With this result the second term in equation (76) is considerably simplified, 
since the various matrix elements involving the electrostatic operator 
(2) 

H ' will vanish when the atomic states considered correspond to the initial 

q 

state. Identifying these factors in equation (73) and setting the terms containing 
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( 4 ) 

them equal to zero, the expression for E may be readily obtained. Since 

( 2 ) ^ ( 4 ) 4 

E is zero, the resulting expression for E ' is just AE/e . Applying 

^ ( 4 ) 

the above comments and rearranging terms in the expression for E , one 
obtains 
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The terms in equation (83) are then the nonzero quantities which need 
to be evaluated to get the interaction energy of the system to the fourth order 
in the electron charge e. The interaction diagrams corresponding to these 
terms are summarized in Figure 14, in which all the possible nonzero com- 
binations are included. Referring to Figure 14 one can see the reason for 
rearranging the terms in equation (73) before writing down equation (83). 
This rearrangement yields the following groupings: 


a. 



b. 


The first diagram corresponds to the coulomb type interaction 

.;»> . 

/ ( 2 ) 

The next two correspond to interactions of the type H 






The next six correspond to interactions of the type 



d. The next six correspond to interactions of the type 




The last 12 diagrams correspond to interactions of ihe type 


H 


( 1 ) 



H 


( 1 ) 





The above diagrams in Figure 14 are now used in the calculation of the matrix 
elements of equation (83). 


Evaluation of Terms in Equation (83) 

The terms in equation (83) are evaluated by considering the terms in 
the order in which they appear. The first term is given by 



Because this term corresponds to the electrostatic interaction between the two 
atoms, it has been considered previously by many authors [6, 13]. The 
procedure used in evaluating this type of term is simply outlined, and the 
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result is given in terms of quantities used in this calculation. In evaluating 
quantities of the form given by equation (84) , various techniques [12, 14] may 
be used to obtain approximate results since the exact evaluation is quite 
difficult. The usual approximation technique is known as the Unsold 
Approximation [11], which consists of replacing the denominator of the 
quantity in equation (84) by a constant average energy e. Once this is done 
the expression* in equation (84) is evaluated by using the matrix summation 
rule over complete states. In this approximation, equation (84) reduces to 



An alternate approximation used in this calculation consists of restricting the 
intermediate atomic states to the 2-p states. Using this approximation 
(which, incidentally, is the same used by Casimir and Polder, and Power and 
Zienau) , the expression in equation (84) becomes [15] 



The factoring of the denominator from the sum over m is due to the fact that 
now the sum over m is just over the triply degenerate 2p states. Summing, 
the above expression reduces to 


AE 




Since the electrostatic interaction operator does not connect s-states. 


«“ I I “» 




vanishes when a > represents the Is ground state. 


In this case both approximations yield the same result, if one lets 

e = 2(E -E j. In subsequent terms the latter approximation will be much 

' ^ ^ ' ( 2 ) 
easier to handle; even though, in the terms involving H , both approxima- 

Q 

tions will again yield the same results with only a different constant. 
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Having introduced the approximation techniques to be used, one 
proceeds with the calculation of terms in equation (83). By substituting 
( 2 ) 

H from equation (3) and factoring the atomic eigenstates as before, the 
Q. 

preceding equation becomes 


AE 


(»D 


16tt^ 


y Y (-i) 

2/ E - E ^ ^ ^ 

^ 0 L,i,m Tyfx 


i>r 


(L + i)r(r + Y)i 


L+i+r+v+ 1+1 


I A*"' I -2^"' hl.0<°^> 


[ (2/ + 1) (2L+ 1) (2y + 1) (2r+ 1) (y -n)Hy + ii)l(T - f<): (r + !(i - m) !(i + m) !(L - m) • (L + m) ! ] 




Y™ (I) (I) 


Y°(I) 


Y°(n) 




y"(h]> . 


Using the coupling rule for Spherical Harmonics [10] and the result 



(- 1 ) 


m 


4 7T 


6. 6 

m, -ju 


the above simplifies to 


AE 


(-1)^^ {(L + i)! 


-III , , : 

^ ^ f L i m 2 - E J R 


iK^l.oW I --f I«1.0<^^ <^1 .qW| rf Ut 

|^(2i + 1) (2L+ 1) (I - m)!(i + m)!(L - m)!(L+ m):J 


2L+2i + 2 


(85) 


The radial matrix elements may be evaluated using the eigenfunctions defined 
by equation (21). The result is [16] 



( 86 ) 


The expression in equation (85) will be used to obtain the various multipole 
contributions due to the electrostatic interaction energy. For example, if one 
wishes to obtain the dipole-dipole approximation, one lets L = i = 1 and sums 
over m = +1, 0, -1. This will be done later when one has the complete 
expression for AE. 
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The next term in equation (83) to be calculated may be redefined as® 

■> 4' I I <.> 


/ \ / V /a I I a'y <a» ; 

\ O' O'* / 


(87) 


Expanding the above term, one gets 

S (l T" 

/ \ cr"^o? I - j^E^l 

X <(u 


[uomo) 

X <^(n) II(m) 


<^0 


\ 

I 


(I) + E (II) + NJi'ck + N'Ji’cK'll 
m J Jy 

N^(T), . . . N4(r' ) , . .^ I I(n) n(m)^ 
^ 1(0) U(0)^ 


(88) 


N (K), ... N' (k'), .. 


r(2) 


which can be simplified since H' depends only on the coordinates of atom I, 

( 2 ) ^ 
and Hjj depends only on atom II. 


To find the form of the intermediate photon states, one recalls that 


(2) 


H consist of four different products of field operators, as shown in 

i, 11 

equation (52). Using the relations given by equations (33) to (36), it can be 
seen that only one of these products can have a nonzero solution; that is 


a^(K) (2^,(/c’) 


N^(k), ... Nj^,(/c’), ... 


= 'Tn (k) Vn' (?) 40 
A A 


(N-l)^(/c), ... 


, (89) 


which is satisfied for the case in which N = N' = 1. Thus, the intermediate 


photon state may be given by 


l^(/<), ... • An equally suitable 


3. This type of definition will be used throughout the remainder of this 
calculation. The form used indicates the type of operators associated 
with the various interactions. 
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solution is given by ... . which is obtained by inter- 

changing the respective photon parameters. Thus, the most general solution 
is (2) I + j . This intermediate photon state picks 

out one term corresponding to k,X and /c’,A* from Having established 

the form for the intermediate photon states, one evaluates the matrix products 
in equation (88) as follows: 


( 2 ) 

Taking terms for ^ that give nonzero results and substituting into 

the respective matrix products, one performs the indicated multiplications by 
using the approximations indicated before. Having done this, and having 
operated with the field operators, one obtains the following for equation (88): 






^i(K+K’)-Pl 1,0 


H(0) 


-i(/C+K’)'P2 


11 ( 0 ) 


(90) 


After rearrangement of the quantities and substitution of and Pz in terms 
of the relations given in equation (47) and (48), one obtains 


JxJj,' .^'‘'<'^ + '‘9 


“i{K+K^) ■ R 


1(0) e 


i(K+/f') • rj 


I( 0 )> <. n( 0 ) e 


-i(K+/<*) • T2 


U(0) 


(91) 


Before one can simplify the above result, one needs to replace the sums over 
K and k' by using the well-known substitution [17] 


Z- / 


Vol. 

(i^ 


= l2^ /dfi^ 


where (27 t) (Vol. ) is the density of individual states in /c-space. Using 
the polarization relations given in Appendix A, equation (91) becomes 


= -(iW?) _ KK'(K+ o' _ ( 92 ) 


<^( 0 ) I 1 ( 0 ^ 


1(0) ><n(0) I n(0) 


-> 
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The preceding equation could have been obtained much easier by considering 
the interaction diagram in Figure 14 which corresponds to this term. If this 
had been done, steps required to show the form of the intermediate photon 
states could have been omitted since the interaction diagram shows that the 
intermediate state consists of two photons and that both atoms are in their 
respective initial states. Nevertheless, the answer would have been off by a 
factor of 2 since we would have neglected the possibility of interchanging k' 
and K.^ This remark applies also to the third diagram of Figure 14 and to 

( 2 ) 

the next six diagrams that follow, since they too contain the ^ operator 

which gives rise to instantaneous emission and absorption or absorption and 
emission of photons. Hence, each one of these diagrams corresponds to two 
possible modes which are indistinguishable and must be taken into account by 
taking linear combinations for the photon states as was done for the term just 
considered. A variation of this situation was discussed in conjunction with the 
remainder of the diagrams in Figure 14, in which there are two possible non- 
zero diagrams for each of the terms in equation (73) belonging to two possible 
unique ways to go from the initial state to the final state. Let us proceed to 
obtain the rest of the terms in equation (83) using a more direct approach. 


The next term to be evaluated is very similar to the term just con- 
sidered. This can be seen by referring to the third diagram of Figure 14 
which corresponds to this term. The result is given by 



^ 2® ir* ^ ^kk'(k + k') 

<n(0)| |u(^<<(0J( e 


+ i{K+K 


*)’R 2 /^ 

1 1 + COS^(/C,K*) 


~i(/c+/c') ■ Ti 


1 ( 0 ^ 


(93) 


which is the same as the previous term, except that the exponentials have been 
replaced by their complex conjugates. Before computing the remainder of the 
terms in equation (83) , the terms just evaluated are combined into a group, 
(The reason for this will be made clear later. ) Combining equations (92) and 
(93), one defines 


X(l) 



9 


4. This point was discussed by correspondence with Prof. E. A. Power, 
London College. 
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where 


\) (94) 

^ ^+i(k+k»)-R/+i(«+k»)- r2\/-i(K+K')‘ r^X I . 


The next term to be considered corresponds to the fourth diagram in 
Figure i4. This term is evaluated using the previous results and making use 
of the information contained in the interaction diagrams. 


This term is defined as follows: 



By referring to the interaction diagrams and substituting for the operators 
and only those terms which give nonzero results, one obtains 




‘III 


K\ »f'X' m 


Me (K + k') [(Ej - Eg) +ifc« 




X 


X 


^(0) II(0)j <o| |l(0) U(0^ 

<^(0) U(0)| <^X:K'X'j |k^ |i( 0) U(m^ 

<M 0 ) U(m)| <^X I (-7 'x | ^ | 


(96) 


Operating on the photon states, substituting for and P 2 and rearranging 
terms, the above equation becomes 



(97) 
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The momentum operators may be replaced by position operators by making use 
of the well-known approximation [18] 




JJL 

li 


(E. - E ) r 
i m 



(98) 


Making this substitution, replacing the polarization products using results in 
Appendix A, summing over the atomic states, and replacing the sums over 
K and K*, equation (97) becomes 



f dK f d7< 


-i (k+k^)' R 




E^) +hcKj 


X 2 I ~ + (k) . • K'} I 

X <I(0) I i(0p> <4(0) 


e 2 (r 2 )^(r 2 ). 


n(0) 


(99) 


The next term to be evaluated corresponds to the fifth diagram in 


Figure 14 and is defined by 


II ^ I ^ I 


One notes 


that this term is identical in form to that given in equation (99) with atoms I 
and n interchanged. If the coordinates centered on atom I were replaced by 
those centered on atom II in equation (99) , the result would not be correct 
jmless the sign of the exponential factor containing the internuclear separation 
R is also changed. The necessity for this change is that interchanging atoms 
also requires the change of the direction of the vector R, because the assumed 
convention calls for R to be directed from atom I to atom II. Hence, this 
term is given by 



(E,-E^) 


+i ( k+k ') • R 


J / I kkMk + k’) [(E^ - E^) +/«ck] 

2 j " + (ic) . (fC*)j(K . !c*) I 


X 





i(o)> <n(o) 


i(K+K') • T2 


U(0^ 


(100) 
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The next term to be evaluated corresponds to the sixth diagram in 

Figure 14 and is defined by . Using the 

interaction diagram corresponding to this term, and performing similar 
operations as for previous terms, one obtains 


^VnV-yiV-yn “ZEE 


(2 7tK)2 

f . .1 1 

’voi7ifpP“ 

j kk' j^(E^ - E^) +-KcKjpE^ - E^) +ilc/£*j 


X <^(0) n(0)|<^ I SjK') I K'^ 1 1(0) n(m)^ 

X <^(0) n(m) j <^k'V I S^MS ^.(k"') e***'" (it) • e^, (It') j j 1(0) II(m^ 

x<^(0) U(m)j 


I q-Mk*P2) ^ /T- 


e^(«)'P^(r2) I 0> 1 1(0) D(0) 


( 101 ) 


The form of this term is very instructive in that one sees quite clearly that a 
term with k,X and k\X' interchanged is also a possible combination. Thus, 
even though the intermediate photon states are not made up of linear combina- 
tions of K,X; k\X\ one still has two possible modes due to the instantaneous 

( 2 ) 

emission and absorption process associated with the operator. Per- 

forming similar operations as before, equation (101) becomes 




2^7? uc^ 


fXC'- J 


f ^ ^ r dp 
J K J k’ 


-i(K-K')-R 


j^(E^ - E^) +^1 ck] [^(E^ - E^) +ilc/c’] 

^ ja.. - (k).(k)^ - (k').(k')^+ (k).(k').(k • Ic‘) I 

<^(0) j 1(0^ <n(0) 


g-i(K-K0T2 I 


■> 


( 102 ) 
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The next term is obtained directly from this result and is given by 



■<VV’ 

2^ II* II 


r dfc r QK* 
K k' 


dK' 

/c* 




-0' pi-V 


E J +-fic(, 


X 


X 


2 |*ij - • «')| 

<(l(0 




1(0) ><n(0) 


i (k-k') • T2 



(103) 


The next term corresponds to the eighth term in Figure 14 and is 
defined by ^ . Carrying out similar operations 

and simplifications as above, this term reduces to 



2^n*HC^Yi J 


r dK r QK’ 

J /c J fc' 


?i. 


-i(K+K*) • R 


(K + /<■) [(E^ - Ejj) +hcK'] 


E l A A A A A, A A A 

|6y - (k)j(k).- (k')j(k')j+ (k).(k')j(k • k' 


X 



^i((c+/c') • rj 





(k+k')t 2 



(104) 


The above quantity could also have been obtained from equation (100) by taking 
the complex conjugate of that expression and reordering the various matrix 
products, provided one notes the interchange of intermediate states o;' and 
O'". Otherwise, the denominator of the expression in equation (104) would be 
in error. 


The next term to be calculated corresponds to the ninth diagram in 


Figure 14 and is defined by 


<A.p\<(A.P>n4-^>i . 


Comparison with 


equation (104) shows that this term may be evaluated using the previous method 
of interchanging the atomic coordinate systems. When this is done, the above 
term becomes 
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j K J k' I (<c + <c') [ (Ej - E 


■R 

E^) +J(ck'" 


<■ 


(K). - (^).(K‘). + (J).(^)^(2 


^)j 


X < 1(0) e 


-i(K+K*) • rj 


l(0)><n(0) 


n(o^ 


(105) 


Before proceeding to the next term, one notes that the basic differences 
in the last six terms just evaluated occur in the energy denominators, the 
exponential factor in R, and the matrix products over atoms I and II. Com- 
bining equations (99), (100), (102), (103), (104), and (105), one obtains 
an expression for this group of terms which is defined as follows: 



X(2) 


(E,- 

■F 


+ (K).(K^).(K- in 


x| g-i («+«') -R - - \ 

j (K + /C') [(Ej -Ejj) +hc«J X / \ (r2).(r2).^ 


+i(K+«')-R 


+ T- ^ /r, \Z+i(K+K')-r2 

(k + k') - E^j) +JicKl X 


hce 


-i(K-K»)-R 


[(E^ -Ej,) +iic/c] [(E^ - E^j) +licK’ ] 


iic e 


+i{K-Kn ‘ R 


[ (E^ - Eg) + JICK (E^ - Eg) + Kc/c'] 




+ 


+ 


g-i(^r*)-R 

(/c + K») [^(E^ -Eq) +Hck’J 
^+i(K+K*)-R 

(k + k') [(E^ - Eq) +JiCK*J 


■ (72)^ (72).^ 


( 106 ) 
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In this term, the matrix elements are to be taken over the appropriate eigen- 
states. In subsequent discussions this group of terms will be referred to 
simply as X (2) . 

The next term of equation (83) to be considered corresponds to the 
tenth diagram of Figure 14. This is the first of a group of six terms corre- 
sponding to the third row of diagrams in Figure 14. Due to the great number 
of operations involved in obtaining each term, only the first term in this group 
will be evaluated in some detail. After this is done, the results of the other 
terms within this group will be expressed as a sum of terms defined by X(3). 
Let this term be defined by 



Expanding this expression and substituting for the operators 

only the factors which give nonzero results, the above expression becomes 



Operating on the photon states, ^llectii^ terms, substi^ting for the momenta 
using equation (98), replacing pj and P 2 in terms of R, summing over 
atomic states, replacing the sum over k by its corresponding integral, and 
substituting for the polarization relations using the results of Appendix A, the 
above expression becomes 



(E -E )^ 
^ 1 0^ 

-ik'R V 1 
r dK ij 1 

1..- 



■ 

^ K 

pE^-E^) +iic/c] 

Fe -e 

L f 0. 

1 (108) 


X 




, -iK-r, 
(ri)^e ^ 
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Performing similar operations on the next five terms and using similar 
definitions, one can combine these terms as follows: 



where 


X(3) 






-i/c-R 




e “iK-r, . +iK«r2 \ 


„ -i/c-R 

/ +1 K • — 

■(E^ - Eq) +Jic«][(E^ - Ejj) +Jic«J X i 


-- H<“> (7,)^y 


-2e 


+iK-R 


(E^ - Eq) +^CK Jj^(Ei - Eq) + Kck J \ 




-i/c* R 


e 


[(E^-E^)][(E^-E^)+ifc«] 


„(2) 


Mri), e 


+1/C-R 


'(^1 “ V][^^ 


^ / -iKTt +i«^r9 ,.(2)\ 

:,-E„)-«o,] M ' "■Ij «, / 


(109) 


The above group of terms shows, in a systematic manner, the way the field 
operators and the electrostatic potential are coupled to produce nonzero 
contributions to the interaction energy. These terms will be discussed later 

( 2 ) — — 

in jnore detail, by expressing the operator H in terms of r^ and r£, and 

i«- r . ^ 

e expanded in a power series. 
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The remainder of the terms in equation (83) to be evaluated correspond 

to terms involving products of . They are illustrated by the last 12 

diagrams of Figure 14. These terms have several common factors and will 
be combined into a group as before. Since the number of steps required to 
evaluate each of the 12 terms is rather large, only the first term of this group 
will be evaluated in some detail. The evaluation of subsequent terms is done 
in a similar manner with only minor deviations. Let this term be defined by 





(110) 


Using the interaction diagrams to obtain the form of the intermediate states, 
and substituting for the operators , only those factors which give 

I, JJ. 

nonzero results, equation (110) becomes 


4-P> J <^.P> J <A.P> y <A.P> JJ 


V V V 


I. 


a- a"''” \ PS " ^^1 ^ ^0 ''•>] [ “^0 ' '-'^0 - (Ej p E^, . licK)] 

X <I(0) U(0)| <0 I I h-) nto> 

X <^(m) H(0)| I ^^,(F-,.Pj(7., ] I k-x,k-a^ |i(0) U(0))> 

X <^( 0 , II(0)| <Ca.a-A-| e-‘'^’-'’=a^,(?,.Pjj(L>] |A;y>ji(0) uouj 


X <U0, nun) I < 1 -A I (.^) e-‘''-'^^^^(K).Pjj(r,)] | o)>|jW, n,0))> . (HI) 
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Collecting terms, operating on the photon states, substituting for pj and P 2 , 
replacing the momentum operators, substituting for the polarization sums 
from Appendix A, replacing the sums over k and /c*, and summing over the 
atomic intermediate states, the above expression reduces to 

<A. P>,<A. P>j<A. P> ji<A.p> 


- v‘ ' 

r dK 

r dK* 

-i {k+k'< ■ K 

. e . . _ 


j — . 

i ~ I 

[(E^ - Ep) +Kck][/c + k'J [(E^ -E^) +Hcx] 


x I |«ls - {■’jJ - } 


x<I(0)|e*' ' ‘Mri).(7i). 1(0)^ 


i(K+K')-ri 


xc n(0) ( 7 ^)^ n(0) 


(112) 


The remainder of the terms in this group are obtained in a similar manner 
with only minor modifications. Defining this group of terms as X(4) , one 
can write the results as follows: 


X(4) = 


(V^ 


/ ^ E S 






J. 


[(E^ - E^) +Kc/cJiiCK + hc/c'J|^(E^ - E^) + HckJ 

e <e (ri)j (d) <e 

- Eq) +hCK'J[^hCK +4icK'][^(E^ - Eq) +11ck] 
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(T.), (7,), (7,)^ 

* [(E, - E,,) [(E, - E„)+(E^ - E„)][(Ej - E„> *Kc«] 

^ pE^-E^)+Kc/c] |-(E^-E^)+(E^-EQ)+JfcK+JlCK' ]p^ -E^)“+ilCK'] 

■" [< - ®o) -*<=«][(E,- V* 1=1 - V][l^ - V *“’] 





S' • ‘ (t,), (r,)^)> 


i(K'+K')-r. 


j^(E^ - Ep) +^CKjj^ifc«- + .«ck'J[^(E^ - E^} + Jic/cj 


(113) 


The above terms are listed in the same order as the corresponding diagrams 
of Figure 14. The reason for expi*essing the denominators in this manner is 
to indicate explicitly the intermediate state energy denominators and to keep 
track on the various terms. A check may be made on the exponential factors 
and energy denominators by comparing the first and last terms in equation (113) . 
One should be able to obtain the exponential factors of the last term by taking 
complex conjugates of the imaginary factors in the first term. The remainder 
of the terms may be checked in a similar manner. 
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Combining the above results, the interaction energy correction AE/e^ 
is given by 


AE/e^ = AE(H ) + X(l) + X(2) + X(3) + X(4) 

q 


(114) 


where the above quantities are defined by equations (85) , (94) , (106) , (109) , 
and (113). Note that to this point the only approximation one has made involves 
restricting the intermediate atomic states to the 2p levels; otherwise, the 
above result is good to all multipole orders both in the electrostatic potential 
and the radiation field. The various approximations will be examined in 
subsequent discussions. 


D I POLE- D I POLE A P PR OX I MAT I ONS 


Approximations to the Interaction Energy 


The results obtained in the preceding discussion are used in recombining 
the individual groups given by equation (114) into a form suitable for subsequent 
approximations. The dipole approximation will be applied and the results 
compared to those of Casimir and Polder. The results will be given in terms 
of corrections to the interaction energy due to only the electrostatic interaction. 
In the dipole-dipole approximation this will consist of taking the terms 


corx'esponding to = L 2 = 1, and M = 0, -1 


in the expression for 
ire- r 


H 


( 2 ) 


in equation (3) and setting the exponential terms e‘“ ‘ equal to imity in the 
various matrix elements in equation (114) . Since higher approximations will 
be considered later, the initial results will be put in a form suitable for use in 
subsequent applications. 


The terms in equation (114) will be evaluated in the order in which 
they appear. The first term is already in the form desired; hence, it need 
not be considered. The next term X(l) is defined in equation (94). Inter- 
changing the coordinates for atoms I and II, and substituting for the volume 
element d« in terms of fc^d(cdf2 , X(l) becomes 

K 


/ \ ? , r 7 2 r {1 + cos2(^ K')l 

^ ' 0 47T 0 47t 


(115) 
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The above expression may be factored by noting that the wave vector fc is 
given in the same' coordinate system as the electromagnetic v^tOT potential 
A(p). Since R is along ^ (Fig, 2), the exponential factor «.R is just 
kR cos (k,R). Additional terms depending on the angular coordinates of k 
are contained in | 1 + cos^ {«,«') } and the exponential factors in the matrix 
elements. Hence, one needs to decide on the degree of approximation to be 


i /c • I* 

used in the retardation factor e before performing the integrals over 

i /c • r 

df2 and dS2 , . If the dipole approximation in e is used then the only 

K K 

factors depending on the angular dependence of the wave vector k are the 
remaining two factors above. To determine the number of terms to be included 
i K • I* 

in an expansion of e , note that the coefficient of the term in equation (115) 


is proportional to 




If the highest approximation is to correspond 


to quadrupole-quadrupole orders, then the terms proportional to 


or 



which correspond to dipole-octupole 


orders in r 2 and rj and to octupole-diple orders in r 2 and rj must be 
included. Henc^ one needs to retain the first five terms in a power series 
i /c ‘ I* 

expansion of e ; that is, 


iK* r s 


1 + iK* r + — (k-v)^ 


(k- 1’)® + («• r)‘ 


(116) 


Letting (k+k') - w and using the above expansion, the matrix products of 
equation (115) give 



^i (k+k') • r^< ^-i(K+K') • V2y 


1 + i(co.ri) (cj.ri)' 


(o)-ri)' 


+ (wri) 


1 


1 - i (cu- rj) - — (cu- r2)^ + — (w rj)' 


+ (cu-rj)** 


(117) 
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Neglecting combined powers of r greater than the fourth power in r^ and 
the above becomes 


</gi {7+7 ' ) . 7^y {7+7 )>7^y 

= <(^'y - -J + 'J ^ (w. rj)^ j ^ 

+ i<^-ri”^<^jl - i (w.^g) --y + y 

- ^ ^ (w- rj) - -j j ^ 

- ^ ^ ^ 1 ^ 

+ ^<(^J-7.)‘><{lj> . U18) 

These matrix products simplify considerably since 

<w-7i y=Yj ® » 

i ^ 

<({7.7i)y= Yj ° ’ 

ij,« 

where (rj) . are the Cartesian components of the position vector r^ . The 
first term is readily shown to be zero using the definitions of (r^)^ in 
spherical coordinates; <\(co- rj)®^ is shown to be zero since 
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In addition, all the various combinations of the (r)^ components are expressed 

in terms of Spherical Harmonics® which are then integrated as sets of three 
products each. Incorporating these results, the nonzero terms of equation (118) 
are given by 

<^gi {«+«’) • ri^<^g-i(K+/c') • r2^ 

= 1 ‘-T E “t“s 

I ts 

ts qn 

ts qn ^ L ^ 

The above sums are evaluated using the results below, obtained for the matrix 
elements over (r^)^ components.® 


5. 


(ri) = - 
1 


(rj) = 
3 




^ r. y.« 


6. The terms <^(r) . (r) ^ (r)^ (r)^^ are evaluated using the Spherical 

Harmonic Addition Theorem twice on each set of two, and using the fact 


that 
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where ^ and ^ ^ are matrix elements taken over the radial portion 

of the atomic eigenstates given by equation (86). To obtain results in terms 
of K and k' and their corresponding angular coordinates, one replaces w 
by {k+k') at the end of each calculation. For instance, the second term of 
equation (119) may be expanded using the results of equation (120) in which 
case one gets 


E "s > “ T > E “t "s 

ts t, s 


= <^rj^ + 2 KK* cos (k,k')J . (121) 


Following the same procedure, the terms of equation (119) are expressed in 
terms of k , k ' and cos ( k , k '). Substituting these results in equation (115), 
one obtains 


X 


cos^ («,«')/ e~ 


\ / 0 47T 0 47T 

X I 1 - {k^ + k'^) + 2kk' cos (k,k') J 


kk' (k + k') 


* ("t) [ (i") ' + (i) (-i) ( <^^2 >)] 

X ^{/C^ + k '^) + 4/Ck' ( k ^ + k ’^) cos (k,k') + 4/C^K'^ COS^ (k,k')J I 
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Defining cos ( k , k ’) = cos 0 


equation (122) may be written as 


X I (1 + cos^ e) - [ (f<^ + {1 + COS^ e) + 2 kk* (cos e + cos^ 0)J 

X + k’^) (1 + oos^ e) + 4kk' (k^ + k'^) (cos e + cos® e) + 4k®(c'® (cos® e + cos* e)J I . (123) 


In the above result, the various approximations are separated as follows: the 
first factor, containing only (1 + cos^ 0), corresponds to the dipole-dipole 
approximation; the second term containing a more complex dependence on 0, 
corresponds to the dipole-quadrupole approximation; and the last term 
corresponds to the quadrupole-quadrupole approximation. In subsequent 
calculations the various X(i) terms will be put in this format. 


The next term to be considered is X(2), defined by equation (106). 
X(2) may be recombined by interchanging the coordinates of atoms I and II in 
the second, fourth, and fifth terms; this results in the following expression: 


X(2) = 


2(E^-Ep)® 

2Vmc“X 


/ ^ ^ 2 ja,.- 

^i(r+r!)-^Xi(7+r')-T2 - 

(K + «'} [(E^ - Eg) + Hck] \ / \ (^2>i 

[(E^ - Eg) + JIckJ |^(E^ - Eg) +iiCK'J X / \ (^2)i (i-2)j^ 

)-H y-i(7+7’>-T\y+i(7+7’)-T2 - - \ I 

;« + K')[(E^-Eg) +JfcK'] X /X ■ 


(124) 


Noting that the coefficient of the above expression is proportional to , 

and that the factors ( 1 * 2 ) ^(r 2 )^ in the above matrix elements give results 
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proportional to ^ one needs to retain the first five terms in the 

i K • r 

expansion for e in order to obtain results accurate to quadrupole- 
quadrupole orders. To simplify this expression, one needs to show that the 
i K* r 

expansions due to e for each of the terms in equation (124) yield similar 
results. The matrix product (letting co = k+k') corresponding to the first term 
in equation ( 124) is given by 



( 7 + 7 ') 


7,y^-i(7+7).7, 


+ <^jl + i(W.7i) (W.7l)3 + ^ )^<^i(W-^2)(^2)ify2)j^ 

+ ^{ 1 + i(w- rj) - -y (^- r^l)^ - (w- 7i)^ + (w-7i)^ <^- (w-72)2(^) . 

+ <^{l + + ^ (“-1^2)^ (l^2)ify2)j^ 

+ <^jl + i(w-7i) (w-7i)2 --i- (7-7y + -^ ("•i^i)^|^<^-^(^-^2)^i^2)j(i^2)j^ 


Retaining only the appropriate terms and neglecting the terms proportional 
to <^(7).^ and <^(7).(7).(7)^^ , one gets 




(125) 


7. 


1 /c* r 

These factors alone, when e 
approximation. 


1, give the dipole -dipole 
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All but the last two factors in equation (125) have been evaluated previously. 


The remaining terms are evaluated as follows. Take (w.r^)3(r^).(^^).^; 

this term is proportional to , where i, j, t, s, 


and q run over 1, 2, 3, corresponding to the x, y, z components of r. 
Expressing the above factors in terms of Spherical Harmonics and using the 
Spherical Harmonic Addition Theorem on each set, one shows that all the above 
terms vanish. Incorporating these results in equation (125) yields 



Letting (k-k*) = v , one obtains a similar expression as given above for the 
matrix products of the second term in equation (124) . The last term in 
equation (124) gives the same results as given by equation (126). This is 
because of the alternation of plus and minus signs in the various products in 
equation (125) . The terms having opposite signs are those which do not 
contribute. By combining the first and last terms in equation (124) , X(2) 
becomes 
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X(2) 




2(E^ - E,) 

~2V 


^-(k+k').R ^+iU+K»).R 

(k + k') pEj - Eq) +*c/fj ^ (k + /£') ^ (E^ - Eq) +ifcK'J 


(t) *ij ■ T 2 ( (i)' “t "s ^s *ij * “t “s 

+ T 2 E ((t)(t) <’^>'> ^"t"s"q"h <<“>'t<~‘^s<'‘'q<“>»h> 
+ (t ) <-^.'> “t "s “q “h < > 

* (T)“t“s“q“h 


He 

■" [(E^ - Ejj) +«0(t][(Ej - Ej|) +«ok'] 

(t) <’^0«ij-T 2 {(t/ <f-*><‘^^'>''t''s^s'*ij'"''t''s 

+ T E eJ ( i) ( t)< * ij n ''s ''q ''h > 

^ (t) <(■?;) 


(127) 


The above expression may be simplified further; but for the present purposes, 
this is adequate since at this point one wishes only to separate the various 
terms corresponding to the various approximations.® 


The next term to be evaluated is X(3) , defined by equation (109). 

This expression may be recombined using the following procedure. Interchange 
coordinates of atoms I and II in every other term in equation (109) and 
rearrange the quantities within the matrix elements, remembering that the 
order of the terms does not matter. After expanding exponentials, one obtains 


8. Note that the sums over i, j are more complex and that one needs to 
calculate <^72) . (r" 2 ) . (r^) ^ ( 1 ^ 2 ) g (^ 2 ) g (^ 2 ) • 
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X(3) = 


2^ 


f ^ 2 |6... 

J K ^ ( y 1 3 J 


(E^ - Eq) [(E^ - Eq) +*c/tj [(E^ - Eg) +^c«] j^(E^ - Eg) +« c«] 

" (<«r (r,)i(r,).^+l (ri),(r,).(7,).^- (^)j(7^)^(7;) 

-T E«t''s[<»r (rO,(rO.(T,)^a).y 

+ T EZZ [<(C (7,)^(70.i7,)^(7,) 


+ (7i)j(7i)^(ri) - 


i(r2)i(7,).^- (7>)j(7,)^(7,)^(7,).(7,).^] 


+ Z Z 

ts Ih 


«[<«r (ri)j(r2)j(72)g(72)^(r2)j^(72).^- 4 (ri) j(7i) . (rj) ^ (72) ^ (r2)|^(r2) , ^ 

+ 6 (7i) j,(7i) ^ (7i) . (72) ^ (72) (72) - 4 (7i) ^^(7i) ^ (7i) ^ (7i) , (72) (72) 

(72)t(72)s(72)j(72)h(72).(72).^] j . 


(128) 


X(3) is obtained by expanding both exponentials in K*r^ and k- V 2 and 

neglecting the matrix terms proportional to the power of r greater than six. 

Neglecting the higher powers is necessary since if one sets exponentials equal 

/ (2) -* — \ (2) 
to unity, the leading term is The lowest term in 

corresponds to terms proportional to (ri)^(r 2 )^ which makes the last terms 

in equation (128) proportional to ^r® r 2 ^^. These factors then correspond 

to the desired octupole-dipole orders. The above results are further simplified 

( 2 ) 

by making use of the r^, r 2 dependence of Hq , For instance, the sum 
S () S ^q^^ adds up to 
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zero, even though each individual term does not necessarily vanish. For 

example, in evaluating the term (r^) (r^) (r 2 ) one replaces the 

components (rj). by their respective Spherical Harmonic equivalents, and 

then substitutes for in terms of and r/"^, and and . 

Factoring the matrix elements into factors depending only on atoms I or II, 
one finds that the conditions for nonzero terms resulting from the Spherical 
Harmonics require the sums over Lj and L 2 to terminate at a given Lj, L 2 
value. Hence, choosing the lowest L, one sums over the appropriate range 
over M. In this case the term gives a nonzero result which in turn is 

cancelled by ^H^^^ (r^) (r 2 > (r 2 ) similar procedure is used to show that 

the fourth factor in equation (128) , consisting of four terms of the form 

(r^)^(7i).(r^)^(r^)^(T2) adds up to zero. This may be verified 

quickly by taking each corresponding set of terms and interchanging the labels 

( 2 ) ( 2 ) 

on the quantities, and, since H (r^, V 2 ) {V 2 , r^) , the sum of terms 

q q 

vanishes identically. Thus, the nonzero terms which contribute to equation 
(128) are given by 


X(3) 


2 ^ c 


J /cdK J 


dn 


0 


z 

ij 


6 .. ■ 
IJ 


(fc)j («) . ^ < 


-iK* R 


1 + Z± - 

’ ^0^ [^^1 " ^0^ [(E^ - E^) +.«c/c] [(E^ - E^) +-Kck] 



(ri).(r2) 




t s 




+ ^ Z Z «h 

ts Ih 

«[<"f (7l).(Tj)^{rj)^(r2)^(72)^(rj).^-4<H^^' ^ (r2)^(r2> 


( 129 ) 
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Later on, subsequent simplifications will be performed when each of the 
preceding terms is considered. 

The next group of terms to be considered is defined by X (4) and is 
given in equation (113). X(4) is simplified using the same techniques as 
before and results in the following expression: 


X(4) = 


2(E^ - 




|1:E^ - E^) + JIckJ ||-Rck +JicK'J ^ (E^ - E^) +^ckJ 

1 

j^(E^ - E^) +*c/fj pCK +i^CK' J j^(E^ - E^) +«C(C'J 

[<®1 - Eq) (E^ - E^)J [-(E^ - E^) 


[(E^ - E^) +ifc(c] [|{E^ - Eq) + (E^ - E^j) +-KC/C + JfcK'j [^(E^ - E^} +XckJ 

1 

^ [(Ej. - Eg) +«ck] [(E^ - Eq) + (E^ - Eq) +Mck + Kck'J [(E^ - E^) + JIck'] 

^ [(E^ - Eg) +*ck] [(E^-eJ+(E^-E^)] [(E;-E^)+«ck'] 


(130) 


Expanding the exponentials as before and retaining the appropriate terms, 
we now have 
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X( 4 ) = 


- E„) 

2' 


■ilc(K + k') 

2(Ej - Eg) 


2 Z {«is - I 

^-1(^7) -R r j 1 

^ic(K + /c*r I [(E^ - E^) +<ficK ■*■ [(E^ - E^) +*Cfcj [(E^ - Eq) +^ck'] 

(i) <-^>(f )<-^^> ^is -T Z > 

tq hf L 

+ 6<^(i^l)t('ri)q(7l)j(7l)j^<](“2)jj(72)f('r2) ^ 

+ <^(i^i)i(“i)j^<1“2)t(“2)q(i^2)b(^2)f(^^ 

^-i(r-7)-R I , 

^ |-<E^ - Eg) +*okJ 2(E^ - Eg) |-(E^ - Eg) +*CK'J 

1 r 1 + 1 1 1 

^ 1^2 (E^ - E^) + ^ck+^ck'J [(E^ - Eq) +'fiCK (E^ - E^) +iIcK' J | 

(t) (t) %^is-i z [<(ri),(' 2 )q(?;)iw;><((^;),(' 2 ),> 


tq hf L ^ 

+ 6 <((^l) t (“l) q (^l) i <^<~2) h <^2> f (^z)i (”^2) ^ 

+ <(("l)il“l)j)><((K)t(“2)q(“2)h<^2)fU2)i I 


(131) 


The above expression has two distinct terms with coefficients 


-i(K+K») . r 
e 


and e ^ ^ ^ ^ ’ as well as various denominators. Subsequent calculations 

show that the two terms may be combined into a simpler expression. This 
will not be done in general, rather each one of the above terms will be con- 
sidered as corrections to various multipole orders of the electrostatic inter- 
action energy given by equation (85) . Once the X(i) terms are defined, as 
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above, one can write the interaction energy AE in terms of corrections 
resulting from various approximations as 


AE/e^ = AE, ,+AE, + AE 

d-d d-q q-q 

Dipole-Approximation AE^_^ 

The terms corresponding to the correction energy AE have been given 
as sums of terms corresponding to the various approximations to be considered 
in this calculation. To evaluate the first term AE^ which belongs to the 

dipole-dipole approximation, one collects all the first terms in the X(i) 
expressions and adds them to the electrostatic dipole-dipole interaction term 
given by equation (85). These terms are obtained by letting Lj = L 2 = 1 and 
m = -1, 0, +1 in equation (85); this yields 


/h<2)\ = y . 

(-1)^ ((2)-)^<'Ri_o(I) rf '•2 

\ ' d-d m=-l 

2(Eq - E^) r8 [(3) (3) (1 - m)'.(l + m) ! (1 - mj!(l + m) !] 


(132) 


Summing over m, and redefining this quantity as X(0) one obtains® 


X(0) 


( 1 ) 


<r,®> 

3R® (E^ - E^) 


(133) 


The corrections to this term, due to X(i), where i = 1, 2, 3, 4, are given by 
the first term in each of the expressions given by equations (123) , (127) , (129) , 
and (131). Combining these leading terms, the interaction energy AE to first 
order (dipole -dipole approximation) is given by 


9. 


Note that is a matrix taken over only the radial portion of the eigen- 

function given in equation (86). 
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d-d 


• <r?Xr}> 


3 R® 


-> <r.^> 1 

- V J 


{ \ r r { 1 \ -i 

\ 2* 7J^ C^ / J K K* \k + k'/® 


-'{'‘■^'‘’'•^(l + cos ^ 0 ) 


_ 2“ 


/f / ^ (J'),(J').. (J),(^)^(^-?')) (i) 




-i(K+K')-R 


+i(«+K').R 


Jjc e 


-i ( k-k') • R 


(«+ k’) [(E^ - Ejj) +Kck] (K + K') |^(E^ - Eg) +*(CK'J ""[(E^ - Eg) +J<CK^ [^(E^ - Eg) +iICK'] 






X 


(E, 


[( E ' - Eg ) . ifc ] " [( E ^- Eg ). J ( C « i [( E ^- Eg ) 


+ Hck 


t) 




F /¥ 2£J^.-<£,,,;i,)[v-i.V''<)(<.t)‘<-.’>«> 


6.- <5. 

ij is 


I ^-i(K+r').R 

1 1 J 

Kc(k + \ 

1 iic{K + k') 

||(Ej - Eg) +J«ck] ‘ ^(Ej - Eg) +J<okT [(Ej - Eg) +2icK'j 1 

1 ^ ‘^1 - V 1 


-i(ff-K»)*R 
" [( El - Eg )^^ CK ] 


2 (E^-Eg)[(E^-Eg)^iiOK-] 


^2(Ej^ - Eg) +XCK +-Bck'J 


( E ^- Eg )‘ 


+ C K 


(E^ - E^) +}iCK' 


( 134 ) 

This expression may be readily simplified by summing over the i, j and i , s. 
These sums contract easily due to the 6-functions included in the respective 
terms. When evaluated, the sums in question give 
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= 2 2 { ) «ij «is 


ij is 
1 + cos^ (k,k') 


( 135 ) 


Further simplification is attained by expressing all the coefficients in terms of 
and <^r 2 ^^. This is accomplished by using the sum rule [1, 19] 


1 ^ 1 

^ <a [P(l,]. m><m 
■ ■ 

[p(D], 

+ <a [p(d], |m)><^m 

1 1 



2 IX 

(E^ - Eq) 


(136) 


and the approximation [15] | P |m > =-«- ("'i - ^m)<' I r I m />, 

given in equation (98) . The sum over I is over the electrons of atom I and 
the sum over m is over the 2p atomic states. Hence, the quantity 1/2 ^c^ 
<rO(E^ - E^) 

is replaced by — r; — rs , which is obtained by carrying out the indicated 

3 (-lie) 

operations above. Using the results obtained in Appendix B for the integral 
over solid angles dS2 and dJ2 , and using the definitions kR = b, k'R = ^ , 

K K 

the second term in equation (134) reduces to 


y 2 \ / 2 \ 

,,, 2^ftc(E, -E^)^<^r “ 


BdB (F(b) FW) + 2G(b) G(g)} 

(b + /3) ’ (137) 


10. 


Reference 19 gives this as 47T/L£ v{a',a) Ka I 7 I = 3X , 


a ' 


where (E ,-E ) = 2 7rhp' (o:’,Q!) 
Ck' Q! 


74 



where 


Fsinb ^,,.1 ^/i.v / cosb sinb\ 

I.(b, ,[__+G(b)J , G(b) . 

where one defines the above expression as in X(0) the superscript 
denotes the degree of approximation as before. Letting (E. - E )R/Kc = a, 

^ ^ (1) 

the coefficient of equation (137) may now be expressed in terms of X(0) as 
follow: Taking 




TT^ 3^ (Kc)^ 


8(E^ - Eq)® R® 

3 7T^ (-K c ) ® 


X 


/ 2\ / 2\ 
~ XIX X2X 

3 r6 (E - E ) 


equation (137) becomes 


X(l) 


( 1 ) 


8a^ X(0) 
37t2 


P J /3d/3 {F(b) F(/3) + 2 G(b) G(/3) } 

J J (b + /3) 


0 


(138) 


Applying similar modifications to the third term in equation (134) and 
remembering that the sums over i, j result in (1 + cos^ 0) , one obtains 


X(2) 


( 1 ) 


- V O 

2*/ (c-h)“ (3)2 


OO oo 


f / K-dK- (^) f 


X { 1 + COS^ 0 


-i(/c+/c*) • R 


+i (k+/c*) • R 


[(El - E^) +Kck] [(El - E^) 


„ / ,V - i ( K - K')-R 

lie (k+k') e 


[(E^ - E^^) +-Kck] [ (E^ - E^) +AicK< ] 


(139) 
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Evaluating the integrals over dn using results in Appendix B and using the same 
definitions, the denominators in equation (138) may be modified to give 


X(2) 


(1) _ 8a^ X( 


(1) oo oo . . 

^ J bdb / [F(b) F(/3) +2G(b) G(/3)] 


(-i- 

1 (a + 


(b + g) 


b) (a + g) (a+b)(a + g) 

The next term to be evaluated is given by 


(140) 


X(3) 


U) 


2 7T^iiC 


f k6k 
0 


(E^ - V [(E^ - Ep) +Xck] [-(E^ - Ep) +;fcK] 

X / I {a... (;?).(^).}<h(2) Q,(73).^ . 

To consider the sum over i, j, one defines I 5. . - (x) . (/c) . \ a.s !P .. . 

(2\ I ^ J i 

Substituting for h'^ from equation (3) this term gives 


(141) 




IJ 




- Z E E E V 

ij L 2 m R ‘ 2 1 |(2 Li + 1) (2L2 + 1)| 




{ (Lj + m) ! (Lj - m) ! (L2 - m) ! (L2 + m) ! } 








(I) 


(rj) 

J J-,2 


(U)^ 


(142) 


The matrix elements containing (r) . and Y are evaluated by taking each of 

1 Xj 

the components (r^). expressed in terms of Spherical Harmonics; the results 
are: 
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J 


(ri) Y^' 



i=l 


1 

4tt 


m.' 





i 

47T 




-1 ^m,+l} ^Li,+1 


A)i< 



i=3 


1 

47T 



(•^M.O \,+l) 


and similar results are obtained for (r 2 )^. Noting that the L^, L 2 
are restricted to Lj = L 2 = 1, equation (142) reduces to 




(-2) (47t) <Cr|^ 

3 R® 


+1 

z 

-J m=-l L 


(1 + m) ! (1 - m) ! 




Performing the above sum over i, j and the matrix multiplications, 
obtains 



(-2) (47t) <r<^ 

3 B? (47r)=‘ 


-(■5-) -‘j 


(143) 


values 


(144) 


one 


(145) 


77 


Substituting this in equation (141) and redefining the constant terms yields 


X(3) 


(1) _ a^X(O) 


(1) 


(3) ( 2 t ^) 


/bdb [ 
n L 


a(a + b) (a + b) 


■ / 


dJ2 e 

K 


iK-R 


f 3 cos^ d 


( 146 ) 


Performing the integration over 9 , cj) , using results in Appendix B the 

K K 

above equation simplifies to 


X(3) 


(1) ^ 4a^ X(0) 
Stt 


( 1 ) 


I 


b^db 
a(a + b) 


j{F(b) + 2 G(b)/ 


(147) 


The last term in equation (134) is simplified considerably since the 
sums over i,j and i,s collapse to (l + cos^0). Carrying out the integrations 
over do and modifying the results as in previous terms, one obtains 


( 1 ) 


X(4)^^^ = J bdb f j3d^ {F(b) F(^) + 2 G(b) G(^)} 


1 r 1 

(b + ^) [ (a + 


^ (2a + b + /3) 

b)2 (a + b) (a + /?) (2a) 


(a + b) |_ 2a (a 


(2a + b + |3) 


+ 13) (a + b) (2a + b + j3) (a + j3) 


. (148) 


The above result may be expressed in a more symmetric form by rearranging 
the fractions as follows: 
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X(4) / bdb / /3d^ { F(b) F(^) + 2 G(b) G(^) } 

^ I (, a + b } [(b + /3) (a + b) ^ (a + b) (a + /3) j 


la+b ^^(b + /3)( 


(a + /3) a (a 


fei-,] 


(149) 


Collecting terms, one obtains the following expression for the interaction 
energy AE/e^ to first order: 

AE^ ^ = X(0) + X(l) + X(2) +X(3) + X(4) , (150) 

d-d 


where the above quantities are defined in equations (133), (138), (140), (147), 
and (149). This expression gives us the interaction energy to first order 
expressed as groups of terms corresponding to various types of interactions. 
The reason for writing equation (150) in this form, as mentioned earlier, is 
to pinpoint the sources of the interaction energy. For instance, the first and 


fourth term 


X(0) X(3) 


depend on the electrostatic interaction 


( 2 ) ( 2 ) 
through its operator H ; whereas, the other terms are independent of H 

^ ( 1 ) 

In addition, one notes that X(3) is due to a mixture of electrostatic and 
field operators and the remainders are due to either the electrostatic operator 


H 


( 2 ) 


or the field operators A.P and/or A* A. These terms are discussed 


further in the following paragraphs. 


Comparison with Previous Results 

The expression given by equation (150) corresponds to the results 
reported by Casimir and Polder [1] , and Power and Zienau [4] . Part of these 
results agree with those reported by Leech [2] in his attempt to solve this 
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problem to first order, bx particular, the term corresponding to X(3) ^ ^ 
agrees with equation (29) of Leech's paper [2]. This may be seen by making 
appropriate substitutions of variables and constants as follows; x -* b, 

. 3 cos X 3 sin X \ 

Xj — X(0) ^ and |F(b)+2G(b)f -*-(^inx + - — ^ J . 

Agreement with Leech's work is also found for X(i) to within a factor of 

2 and with various terms contributing to X(2) and X(3) No further 
comparison is possible since he arranged his expressions differently than is 
done here. Leech's inability to verify Casimir and Polder's results is because 

of this error in X(l)^^^ and because he obtained the limiting cases before 
doing the integrations over /3. These integrations are required to show that 
equation (150) is indeed Casimir and Polder's result. This is done below to 
show the correctness of equation (150). 

To perform the integrations over /3 in equation (150) one considers 
once more each one of the terms containing /3 integrals and then recombines 
the results to obtain an expression involving only integrals over b. The first 

element to be considered is X(l) defined by equation (138) . Rearranging 
factors yields the integral factor given by 


7 bdb F(b) 7 + 2 7bdb G(b) 7 

0 i ib + p) J J 


gdg GW) 
(b + ^) 


(151) 


These terms may be evaluated using the results given in Appendix C. For 
example, the first term above transforms into 


OO OO 

J bdb F(b) J 
0 0 


F(ff) 
(b + /3) 


OO OO 

2 f bdb F(b) f 
0 0 


F(/3) 
(^ + b)(^-b) 


(152) 


A similar transformation may be made in the integral containing G(b) and 
G(/3). Using results given in Appendix C for the /3 integral in equation (152) 
and for the term containing the G function gives the desired results for the 

expression in (151). Substituting in X(l) yields 
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X(l) 


( 1 ) 


8 a® X(0) 



— r 


db 

~V 


|(b® + Sib - 3) - (b® - Sib - 3) e j 


2 ib 


H f I + 2ib® - 6b® - 6ib + 3) 

- (b^ - 2ib® - 5b® + 6ib + 3) | 


(153) 


The next element to be considered is X(2) given by equation (140). 
To perform the /3 integration it is necessary to rearrange the fractional factors 
as follows: 



1 

(a + b) 


+ 


1 

a + /3 


, ■■ (b + ^) _ ~ 

(a + b) (a + ^) _ 


I (a + b + a + /3) ^ 1 1 

I (b + ;8) (a+ b) (a+ /3) (a + b) (a + /3) J 

/ 2 \ ( r 1 1 1 b(a + b) +^(a + g) | 

\ a+ by I L (b + /3) (^ - b) J (b - )3) (b + /3) (a + ;8) j 


(154) 


With this modification, X(2) becomes 

4 (i) °° 

X(2)^^^ = / bdb / )3d/3 { F(b) F(/3) + 2 G(b) G(^) J 


( ‘ 1 

f ' . ' 

1 b(a + b) + ^ (a + ^) ) 

( (a+b)' 

[b + /3 j3-b ^ 

' (a + b) (a + /3) (b + /3) (b - )3) j 


(155) 


By arranging the terms in this manner, one readily obtains the J3 integral of the 
second fraction. Analysis of this fraction shows that it is symmetric in b and 
j3 in all the factors except the one containing (b - /3). Interchanging b and j3 
results in a change of sign. This feature is utilized by taking 
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/ bdb { F(b) FW . 2 G(b) aw I ^ J 


and separating the above quantity into two parts given by 


OO 00 OO oo 

J . . . db / . . . dj8+-|- J . . . db J . . . diS 
0 0 0 0 


then interchanging b and /3 in the second half to obtain 


OO OO OO oo 

— / • • • db J . . . d/3 J • • . db J . , , d/3 
0 0 0 0 


Hence, one sees that this part of equation (155) gives a vanishing contribution. 
Using results in Appendix C, the other portion of equation (155) is integrated 
over /3, resulting in 


X(2) 


(1) 16a^X(0)^^^ 

37.^ 


^ f (a + b) i + 3ib - 3) - (b^ - Sib - 3) 


TT r db I , A <3 n 

4l J 9 (a b) HbH2ib3-5b^-6ib+3) e' 


' (a + b) 


2ib 


- (b^-2ib^-5b^+6ib+3) e 


(156) 


X(3)^^^ is already expressed in terms of an integral over b, hence, one 
only needs to replace | F(b) + 2 G(b) | by its definition. Thus, one obtains 

( 1 ) 


(1) _ 4a^X(0) 


JL f 

’ bdb 

r - ^1 

2i J 

b^ 

|_a (a + b)^ 


X (b^ + 3 i b - 3) e^ ^ - (b^ - 3 i b - 3) e ^ 


(157) 
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The next term to be considered is X(4) defined by equation (149) . 
Before performing the /3 integration it is necessary to rearrange the 
denominators involving b and /3, in a manner similar to the one used in 

simplifying X(2) In doing this, one is guided by the symmetry displayed 
by the quantities involving b*s and /3*s. Considering the fractions in equation 
( 149) and rearranging factors yields 



1 

(b + /3) (a + ^) 


(a+ b) (a + /3) 


+ 



1 

(b + /3)(a + /3) 


+ 


1 

a(a+ /3) 



( (2a + b) 

1 1 

^ a(b + )3)+)3(a + /3) + b(a+b) j 

1 a(a + b)^ 

L(b+ /3) ^ (^ - b) _ 

a(a + b) (a + jS) (b + /3) (b - /3) j 


(158) 

The last term in the right-hand side of equation (158) is symmetric in b and ^ 
in all terms except (b - /3) . Performing the same operations leading up to 
equation (156), indicates that the second term in equation (158) does not con- 
tribute to the /3 integration of equation (149). Hence, considering only the 

first term in equation (158), X(4) becomes 


X(4) 


( 1 ) 


4a^ X(0) 
3 7T^ 


f bdb f |8d/3 [p(b) F(/3) + 2 G(b) G(^)] 
0 0 


I (2a + b) r 1 1 ~| 1 

|a(a+b)^ |_b+/3^/3-bJ j 


(159) 


Recombining the above factors shows that the ^ integrals are the same as in the 

first term of equation (155) . Using results obtained in Appendix C, X(4) 
becomes 
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oo 


(2a + b) bdb 
a(a+ b)'** b** 




/ 


X I [(b^ + 3ib-3) - (b^-3ib-3) e 

+ -^ [ (b^ + 2ib^ - 5b2 - 6ib+ 3) 

- (b^ - 2ib^ - 5b^+ 6ib+ 3) } 


(160) 


With all the terms in equation (150) given in terms of integrals over b and 

X(0)^^^ as coefficient, the interaction energy correction AE^ ^ may now be 

expressed as a function of these new terms by separating the integrals over 
±i b j ±2 i b f 1 1 

e and e as follows: 



db 

, _ 2 a P / b 

\ a / 2a + b \ 

TP 

(a+ b) 2 ^a(a + b)^ 

) 2 V (a + b)^ } 

X [ (b^ + 3 ib - 

3} - (b^ - Sib - 3) e 


4: a? 7 db 

3 7t ib^ 

0 

r_L 

a 1 /2a^ + ab\ 



” (a + b) ^ 4 1 (a + b)^ y 



[(b^ + 2ib* - 5b^ - 6ib+ 3) - (b^ - 2 ib’ - 5b^ + 6ib + 3) e'^ 




( 161 ) 


Note that the first term of the above equation results when only the electrostatic 
interaction is considered; this yields the proper limiting case. By recombining 
the fractional quantities of the above terms, one obtains 


/ 1 2 a b2 

b 

1 

j~2 a^ + a b 

1 - 

( (a + b) 2 

i_ a(a + b)^ _ 

2 

L(a + b)2 J 

1 \2a J 
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Substituting into equation (161) yields the following expression for ^ : 


AE, 



oo 

f 

0 

bdb 

(b^ + Sib - 3) 

(b^ - 3 ib - 3) 1 

1 + 

ah'* 

2i 

2i j 


4^ P db 
35t J Tb*' 


A. ^ 

a a(a 


-b^ 1 1 (b^ + 2ib^ - 5b^ - 6ib + 3) 

rr^J 2i 


(b'^ - 2i b^ - 5b^+ 6ib+ 3) e 
2i 


-2 ib 


(162) 


This is essentially the Casimir and Polder result expressed in terms of 
complex quantities. The singularities appearing in the above expression present 
no problem, because when using stationary state perturbation theory to 
approximate physical situations, one takes the "principal value" of the sums 
or integrals appearing in the results. To transform equation (162) into real 
quantities, one considers a contour in the complex plane which includes the 
interval of integration in equation (162). This is done by taking appropriate 
contours which include the real axis. The first integral term in equation (162) 
is integrated by first considering its factors, as follows. Taking 


I 


db (b^ + 3 i b - 3) e 
a b^ 2 i 


ib 


(b^ - 3ib - 3) e 


-ib 


2i 


and letting b — - b in second integral, as well as exchanging appropriate limits 
of integration, results in the following expression 


f oo o i b 

f db (b^-^ 3ib - 3) e 
■J ab^ 2i 


The principal value of this integral is obtained by using the values given 
in Appendix C; the result is 


T3 T7- r db (b^ + Sib - 3) e^*^ 

•^* J ab^ 2i 
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The results for the other integral terms in equation (162) are 


db [~ (b* + 2ib^ - 5b^ - 6ib+ 3) (b* - 2ib^ - 5b^ + 6ib+ 3) 


/ Hi? _ 

'(ir) ■ 

oo 

p V r - ^ . 

• • J 2ab(a+ b)2 


2i 


2i 


(b^+ 2ib^ - 5b^ - 6ib+ 3) _ (b^ - 2ib^ - 5b^+ 6ib+ 3) e~^ 

2i 2i 


-] 


= -(-f^)+ / — — T <y^ + 2y*+ 5y2+ 6y+ 3) ; 

' ' 0 (a^ + y2) 

I o £ 

where iy = b in the integration containing e , and -iy = b in the 

2 i b 

integration containing e . (See Appendix C for details. ) Combining 
these results indicates that the contributions of the residues at the b origin 

cancel the contribution due to the electrostatic interaction X(0) . Hence, 

AE, , becomes 
d-d 


AE 


^ 4a^ X(0) p dy (/ + 2y^ + 5y^ + 6y + 3) e~^^ 

i-d Stt J / 2 

0 (a‘‘ + r) 


(163) 


Having AE^ ^ in the above form allows one to compare results directly with 

those of Casimir and Polder as given by Power and Zienau in equation (28) of 
Reference 4. This is accomplished by letting y = uR, (E^ - E^)/hc = E and 

recombining the coefficient into dipole moment matrix elements. Making these 
substitutions, as well as using the definitions for a, equation (163) becomes 


AE 


d-d 


4(E^ - E^)^R^X(O) 
(lic)^ (Stt) 


(1) 


f 


R (uR)^ due 


R* 


(El - E^) 


He 


12 TT 


+ u“ 


, , 2 5 6 3 

^ ^ ^ (uR) (uR) 2 TuRp" TuRp" 
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Replacing X(0) in terms of the expression given in equation (133) yields 


AE 


d-d 


y 2 \ / 2 \ 

4 (E - E Sr ><r > «o 4 -2Ru 

1 0 \ 1 / \ 2 / r u^due 

J T 


3^ 7tR^ (^c)^ 




TT 

2 


^ I ^ (uR) (uR )2 (uR)^ (uR)^ j ’ 


The integral form is in agreement with previous results. Comparison of 
equation (164) to equation (28) of Reference 4 shows that Power and Zienau 
omitted some factors of Jic. The reason may be that in part of their discussion 
they let /tt = c = 1. Analysis of the rest of the factors in equation (164) is 
simplified by writing down equation (28) of Reference 4: 


* T.M 4 

AE = - — 

7T 


<q(l)> ^ <q(2)> 2 e2 CO 4 - 2 uR 


^cR2 


I 


u’ e 


d u 


1 + 


(uR) (uR)' 


(uR) 


r + 


[E^ + u^]2 
^ 1 


(uR)^ 


(165) 


A comparison of this equation to equation (164) shows that if one lets 

<Ej - 

= E, (Kc) is left out in the fraction inside the integral. The 

2 / 2\ 

dipole matrix elements ^q(l) ^ correspond to r — . Power and Zienau 

define this matrix element in equation (23) of Reference 4 as 


2 


(166) 
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where they take the sum over the intermediate 2p atomic states. This 
operation is equivalent to the relation used here, given by 






j, 

3 




where m is summed over the 2p intermediate states. By letting q h e r 
and making the above comparison of matrix products, complete agreement is 
reached between the results obtained here and those given by Power and 
Zienau [4] . 


Since the main objective in this calculation is to obtain the form of the 
interaction energy to various multipole orders, a discussion of equation (164) 
will be given after the higher multipole order corrections are evaluated. 

There is a twofold reason for showing the detailed correspondence between the 
results obtained here and those of previous authors: (1) it shows that the 
overall procedure used here is correct, and (2) it provides the steps to be 
followed in subsequent calculations where a detailed discussion is not practical. 


DIPOLE-QUADRUPOLEAPPROXIMATIONS 


Introduction 


The dipole-quadrupole order corrections. 


AE , , to the interaction 

d-q 


energy given by AE are evaluated in this discussion. 


The term AE , 

d-q 


consists of the dipole-quadrupole order terms in the electrostatic interaction 
energy and the second set of terms in each of the X(i) expressions given 
previously. 


Starting with the electrostatic interaction energy given by equation (85) , 
one obtains the dipole-quadrupole order terms by setting = 1, L 2 = 2, 
m = 0, ±1 and Lj = 2, L 2 = 1, m = 0, ±1. Performing these operations and 

( 2 ) 

defining these terms by X(0) yields 
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X(0) 


( 2 ) 


(-1)4 [(l+2).p<R^_^ 

2 

h 

1,0/ \ 1,0 

4 

^2 

“i,0> 

R® |(3)(5)(2 

- 2EJ1 
0 1 ) 



[ (1 - m) ! (1 + m) ! (2 - m) ! (2 + m) ! ] 


(-1)2 j (2+ 1)'. 

iXvo 

4 

^1 


2 

^^2 


R» ((5)(3)(2E^ - 2E^)} 


X 



1 

(2 - m) ! (2 + m) ! (1 - m) ! (1 + in) ! 




(167) 


After performing the summation and collecting terms, one obtains the follow- 
ing expression 


X(0) 


( 2 ) 


/ 2\ /4\ 
< r X < r y 
\ 1 / \ 2 / 

(E^ - E^) 


(168) 


( 2 ) 

The next term corresponding to X{1) is obtained from the second 
term in equation (122). After factors are rearranged, this term becomes 

=[- (3)(2)Vm^c^ ] 4D") f I (t^) 


J J { 1 + cos^ ® + 2 K k’ cos 0 I e ^ ^ ^ 


(169) 


Separating the coefficient of the above term and replacing { — I ( — | 

r/2\ .ir /2\ T 

Vi><^i- V V2><^1- V 


L JL 

expression for the above coefficient: 


, one obtains the following 



(he) V . 

(3)^ (2)^ (Kc)* 


(170) 
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Using equation (86), one obtains the relations 


2 

<r2><r^> = -g- ; 

which may be used to simplify the coefficient given in expression (170) , as 
follows: 


(2) (E^ -E^)2 

(sTTsPWTlh^ 




By using results in Appendix B for the integrals over dJ2 and substituting for 
K and «' using previous definitions, equation (169) becomes 


X(l) 


( 2 ) 


2^a^ X(0) 
(5) (3) 




(b2 + /32) [ F (b) F (/3) + 2 G (b) G (^)] 


- 2b/3 


E?(b) ?(/3) + 4 F,(b) F.,(;S)_ 


h(3 


(171) 


In Appendix C the /3 integrals for this expression are evaluated as follows 
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_ _2ja^ /„ r . 3 .. r /3^d/3 [r(b) F(P) + 2 G(b) G(0)1 

- - (5)(3)«7r2 P J b db J (/3+b) (^-b) 

+ 2 r bdb r ^dgCF(b) F(g) + 2 G(b) G(g)] 

(fl + b) (fl - b) 


3 Y /n\ 


2^a^ X(0 ) 
TkTT5T» 


(5) (3)" 7T 


_4 r bdb r /3^dg[?(b) ?(^) +4F.(b) F.(b)] ] 

^ (/3 + b)(;8-b) j 

J bdb I ^F(b) + 2 G(b)] + •^ [ ?(b) - 4 F3(b)J | 
+ 2i f b^db|[F(b) F'^(b) + 2 G(b) G'^(b)] 


where 


[ J(b) J (b) + 4 FgCb) F 3 '^(b)]| 


1 ( [b^ + ib - 1] e^*^ [b^ - ib - 1] e 

F(b) = [ 


2i 


2i 


) 


( 172 ) 


+ ^ 1 ([b==+ib- 1] ,e\^, [b^-ib- 1] e~^^ 

^ ^ ^ \ 2i 2i 


^ , 1 ( [b^ + Sib - 3] e [b^ - Sib - 3] e 

= -P V 21 - 


(■ 
(■ 


2i 


-ib ^ 


ib 


F+(b) = 4 [ibi±3i^Lz.3L^ ^ Jb^-3ib-3] e 




2i 


ib 


2i 


-ib ^ 


,,(b, . ^ I WJ- lie 


2i 


2i 


-) 
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r+/hl - 1 ( lib - 11 [-lb - 1] ) 

° 21 21 / 

0 


, 1 I [lb® - 2b® - 31b + 31 (lb® + 2b® - 31b - 3] e 

J(b) =-j5sl 


-ib ^ 


1 ([lb® - 2b® - 31b+ 3] e'*” [lb® + 2b® - 31b - 3) e 
J (b) =prl n 2^ 


( 2 ) 

Using these definitions, the result for X(l) is quite compact. In addition, 

the products F(b) F^(b), J (b) 5 (b) , . . . etc, are just the difference of two 
squares, since the binomials merely change sign in the second term. 

(2) 

The next term to be evaluated is X(2) ' given by the second group 
of factors of X(2) defined in equation (127). Expressing the coefficient of 

this expression in terms of and ^ 2 ^ > one obtains 




^-i(x+xV*R («:+«*)• R 


(K + K>) [(E^ -E^) +ttCK] {K + Kn [{E^ +*CK’] J 


fic e 


-i(K-K*)*R 


^1 [(El- [(Ei-E„)+«c«'] j 


(173) 


The sums over i, j and t, s are evaluated using the same techniques as 
outlined previously. The results for the various terms are as follows: 
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E { ) (-t) 2 (t)' ^ij 

= (-"^)(“|') I ((C^ + k' 2) (1 + cos^ e)+ 2 kk' (cos e+ cos^ e) I 



- (ic).(K). 


(«')j(«').+ 


(K).{K').{K 




E 




" (■■2')(~3')('5’) + ) • 174) 

The terms containing the v , v factors have a similar form with only a minus 

V S 

sign in front of the factors 2 kk' (cos © + cos® 0) and 4 kk' cos® 0. The 
evaluation of these sums is very lengthy and will not be given here. The 
procedure is simply to carry out the sums over t, s first and then over i, j. 

When matrix elements of the form ^(r 2 ) ^(r 2 ) ^(r 2 ) ^ (r 2 ) ^ ^ are included, 

the best procedure is to evaluate the various combinations, using the fact that 
the result is invariant under interchange of t, s,i, j. Hence, the number of 
matrix elements that need to be evaluated is reduced considerably. For 
instance, in this case only 15 of the 36 possible terms need be evaluated. Of 
these, only 6 are nonzero and are given by 


<( ( r) ^ ( rO J ( rO J ( r") . ^ > 


i = 1,2,3, 


- (t)(t) • 

(175) 

Having evaluated the matrix elements, it is much easier to calculate the sums, 
since all the zero terms can be systematically excluded, allowing one to obtain 
the above results. This procedure is used in the equations to follow. The 
integrations over dn are done using the results in Appendix B and using the 

fact that the integrals over (1 + cos® 0) ^ give the same result. In 

addition, one finds that the integrals over cos © and cos® © are related as 
follows: 
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( 176 ) 


f J (cos 0 + cos® 0) e 


3 ±i(/c+/c‘)*R 


d£2 dfl , 
K k' 


= - f f (cos 0 + cos® 0) ^ ^ ^ di2 dJ2 , 

K k' 


The significance of this result cannot be over estinaated since the negative 

sign introduced by the factors containing v. , v is exactly offset by the 

t s 

integrals over dS2 , dS2 , . Incorporating this simplification in equation (173) 

K K 

and using previous definitions as well as performing the integrations over dS2 , 
one obtains the following for equation (176). 


(2) <>^* > 
X(2)' = j r ~ - 


(3) (2)^7^«cR‘ 


/ bdb / 136^ (-^) 


I i 1 (b+ /3) 

^ I (a + b) ^ (a + /3) ^ (a + b) (a + /3) 


(- - y ) (-j)^ <ri^><r2^> ^(b^ + (2®^^) [F(b) F(^) + 2G(b) G(^)] 

+ 2b^ (-2®7r2) r?(b) ?(^) + 4 F,(b) F3(^)] 


(- 4'X'3‘)('5'X''2^> f(b® + |8®) (2®jr2) [F(b) F(/3) + 2 G(b) G(/3)] 


+ 4b/3 (-2*7 t2) 


4^2. ?2(b) 52 (/3) +6F.,(b) F,(g) 

L bj9 


(177) 


where 


Fa(b) =F(b) + 2G(b), 5(b) =b®G(b) - Fg(b), 52(b) =b®G(b) -2 Fg(b) . 
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( 2 ) 

To express the above results in terms of X(0) , one uses the relations 

between and ^r^^ given previously. Hence, the above coefficients 

may be rearranged as follows: 


a® (2® ir^) 

(3) (2)27r«R^ (lie) 

a® <yi*> (25 7 t2) 
(3) (lie) 


U -L\ (~L\ (J-\ /r4\ 1 - (2)^a^X(0)^^ 

[v 2/\3/\5/\2/J (5) 


Using these coefficients, equation (177) may be written as 


X(2) 


(2) (2)^a^ X^^^ 


(^Tb) 

/ 1 1 \ b(a + b) + j3 (a + /3) ~1 

_\b + /3 p-h) (b -/3)(b+/3)(a+/3) J 


X I 5(b2 + ^2) [F(b) FW) +2 G(b) G(^)] 

- 2 [2 y(b) y (^) + 3 ?2(b) j2(/3) + 26 Fg(b) F3(i3)] ) 


(178) 


The )3 integrals are now evaluated using the results of Appendix C . Since the 
second fraction involving b's and /3*s has been shown not to contribute, only 
the first factor in equation (178) need be considered. Expanding, this quantity 
becomes 
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oo 


X(2) 


( 2 ) 


X(0) 

3^(57t2) 


( 2 ) 


b^db 


^ / (b/g) - b) ^ ^ 


10 / 


bdb 


J 7b/oU. :bl [F(b)F(g)+2G(b)G(g,] 


^ (a + b) J (b + ^) (^ - b) 


_ 8 r f° — ^ — rj(b) 5(^)1 

^ (a+b) J (b + /3)(^-b) 




- / (b/m^b) 


• (179) 


By using previous results and values, obtained in Appendix C, the above 
expression reduces to 


X(2) 


( 2 ) 


X(0) 
(3=*) (5) 7T 


5 F(b) + 10 G(b) + ^ 7(b) + ^ 72(b) - 


+ i / [lO F(b) F'"(b) +20G(b) G'^(b) _ 

(2) (26) Fs(b) Fs(b) 

V 



r bdb 
J a+b 


F3(b) J 

6 7,(b) 7,'(b) 

b^ 


(180) 
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where 


a, 1 [(ib® - 3b2 - 6ib+ 6) e^^ (ib® + 3b® - 6ib - 6) e"^*^ 1 

•^2W - 2i 2i J ’ 


[- 


1 |(ib®-3b®-6ib+6)e (ib® + 3b®-6ib-6)e 
‘^2(b) = -jgj ' ^ - 


2i 


2i 


-]. 


The next term to be considered is given by the second group of terms in X(3) 
defined by equation (129) and written here as 


X(3) 


( 2 ) 




J K d K f dn 
0 0 






-ifC» R 
e 


(E^ - E^) j^(E^ - Eq) +*ckJ [|(Ej - Eg) +*ck]‘ 




(181) 


The above terms are evaluated using the same procedure used to obtain X(3) 
except here one has no 6 functions in the sums over t and s, to shorten 
the work. Before evaluating the terms in equation (181), one notes that the 

( 2 ) 

first and last terms involving H' may be combined by using the symmetry 

<1 

of the system. Interchanging the coordinates of atom I by the coordinates of 
atom II within the matrix elements of these two terms and expanding the sums 
over i, j and t,s, it can be shown that 


IJ ts ^ 

= Z ?ij Z . 

IJ ts ^ 
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Hence, one needs to evaluate only two terms of equation (181). Rewriting 
these terms yields 




( 2 ) 

Expanding the first term using the definition for H , one gets 

Q. 


y a >.. Y K K y (Lj + ^ 

ii fa * ® Li,L 2 R^l+^2+i j(2Li+ 1) ( 2 L 2 + 1) I ^2 


X 


2 

m I (Li + m) ! (Lj - m) ! {L 2 - m) ! (L 2 + m) ! | 


(182) 


Since the summation indices i,t,j,s go over 1,2,3, one needs to evaluate only 
six matrix elements of the form '^(ri)^(rj)^ interch an g in g 

i and t does not affect the result. Expressing the (r)^ components in terms 
of Spherical Harmonics, one obtains the following relations: 




(r) 




'\W3<“>3'^L ^ J 


1 

/— 6 

4 7T ^ 

/ 15 y m 

1 


4 7T 

V 15 y : 

2 

/T2) ( 27 t ) 

4 7r N 

' (3) (15) 

i i 

^(6 

4 7T V 

15 \ m. 


+ 6 


-JT 


6 6 
m,0 / L,+2 


+ 6 + 
,+2 m,-2 


<5 . 6, 


/X, 

n/ 3 m,0 


L,+2 


.+2 
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(183) 


“ -h JM ^V-i • V+i) \ 
<®3«a^r>= 4^7¥(V->^V^i)\ 


+2 


+2 


The significance of the above results is found in the 6 factor, which 

applies to both of the (r^)^ and (r2)j type terms. This type of matrix 

element selects only the terms having = L2 = 2 in the electrostatic inter- 
action, and thus this term does not contribute to the dipole-quadrupole order 
interaction (Lj = 1, L2 = 2 or = 2, L2 = 1). |^e reason for including it 
here is that it is of the right order in powers of (r)^ and has the same 

general form as the other terms in equation (181) . Since one is interested 
in the dipole-quadrupole type interaction, further evaluation of 

(ri),.(ri) . (r2) (r2) . / is postponed until quadrupole-quadrupole order 

interactions are evaluated. 


The next term to be considered is given by 


T cr,.. Y K K y (~^)^^ (4^ (Li + U)'- 

” ^ ® T T - I ^2Li + 1) ( 2L2 + 1)| 


L1L2 


<^(?.), y!^; (1)^ , Y-f (3))> 


2 , 

m I (Li + m) ! (Lj - m) ! (L2 + m) ! (L2 - m) ! | 


V3 


(184) 


The matrix elements ^ have been evaluated and are given in 

equation (143) . They are proportional to 6 , which in turn picks out 

L.1,+1 

the lowest L value and the corresponding range in the sum over m in 
equation (184) . The other matrix product has three (r2). components, each 

of which can take on three values; hence, one needs to evaluate every combina- 
tion. The terms required are as follows: 


99 



- T /t (V +1 - "m,.l) ' 

*~4ni^Y') J~ ('®m,+l ■ ^m,-l) '’ lj.+I 


i / i \ / 2ir 


47t y 5 3 \''in,+l ^ ^m 


<(r,^(r,i(r)3Y^-*^ = ^ 


1 / 1 


27T 
(15) (7) 

42T 


-i)'^L2.+i 

^m,+2 ^ ^m, -2 } 5 V 7 ^m,0 


L,,+3 


47T V 5 y 'v/ 3 ^m, 0 ^Lj, + 1 


<(r)^(r,2(7,^y-”*>= [y¥’(V^3-V-3) ^ T ^ 

^^(t)(7?)(v^i-V-i)': 


m, -] 


L,, +1 


r; .r; 


<(r)^ 









^ 4ir [y 35 ('^m,+3'^ * 111 ,- 3 ) 5 J 21 (®m. -1 ^ ®m,+l) 

•i7(^)/f'(v..*V-.)‘L,... • 


’'O' -Tt[yTI^ (*».«•*»..-!) -T TT Vo 


L 7 , +3 


1 / 2 




4ir \ 5 / V 3 m, 0 Lj.+l 


47T 


a + a 

) <5 

m,+l 

m,-l 

; l2,+3 

5 +6 

1 ^ 

m, -1 

nij + 1 } 

' L2,+1 


<- 


i"3'n.<«oV“* 


>■^(4)7? 


a a 

m, 0 Ij2j"J"3 


47T 5 /'7 3 ^ra,0 ^Li,+1 


1 / 6 


L2,+3 


^L2,+3 


)_ *1o2,+3 


L 2 , +3 


(185) 
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Analysis of the preceding matrix elements shows that L 2 can be either 1 or 3. 
When referring to equation (184) , one notes that the terms corresponding to 
dipole-quadrupole order are those for which = Lj = 1. These terms give 

results proportional to as required; hence, setting = L 2 = 1 

in equation (184), one obtains the desired expression given by 


V „ y , « (-l)(4ir) (2)! <r,»><r,*> 1 

i ^ “t si 3 rS 

+ 1 <A). (2)> 

(l + m)!(l-m)! 



Yj 9- ■ Z 

ij ts ^ " 


X Z 


1 

(1 + m)! (1 - m)! 




(186) 


The term corresponding to = 1, L 2 = 3 corresponds to quadrupole- 
quadrupole order interactions and will be considered later. Evaluation of 
equation (184) is rather involved since one needs to consider all the terms in 
the sums over i,j,t, s and m. Since the expression (185) consists of 10 

< — m* \ 

(ri)i (1 ) y , 

one must consider 30 different products in equation (186). To show the 
procedure used in evaluating these terms, consider the matrix product 
corresponding to i = t= s = j = l in equation (186) which is given by 


<X(r,) ^ y“’' (1) ^ 1 (Y2) ^ Y-" 


( 2 ) 


> 



= (3) (27t) 

(47r)^(5) 


■ ^m, + l) ( 47 t) ( 5 ) J 3^ 


- 6 

m, - 


.) 
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This result is obtained using equations (143) and (185). Substituting these 
values into equation (186), one obtains 


m 


^11 

(1 + m)! (i - m)! (4ir)*(5) (3) 


(Vi 


m,+l/ V m,+l m,-l/ 


(5) (3) R* Ml ''^1 '^1-' u (1 + m)! (1 - m)! 





1 - («) 



(187) 


Performing this group of operations for each matrix product and defining this 

product as <i><t, s,j)>, one obtains all the elements of equation (186) . Once 

all the terms are known, the sum over i, j, t, s may be performed and the 
final result recorded. An alternate method is afforded by expanding all the 
sums and then evaluating each product. This is not the shortest way since a 

number of matrix products <i> <t, s, j> prove to be identically zero, making 

the expansion of the sums much easier. After evaluation of all these factors, 
one finds that the following are nonzero elements. 


a> <1, 1, i> 

<ri2Xr?^> 
(3) (5) B? 

<1> <1,2,2 > 

i_ 


3 

(3) (5) E“ J 

<1> <1,3,3 > 

1 

r.<ri^<<>i 

3 1 

(3)-(5) R’ J 

<T2> <2,2.2,> 

_ J_ 

3 

r<r,0<r,0l 

(3) (6) R“ ] 

<2> <2, 1,1,> 

_ 

3 

■ <r?> <r» 1 
_ (3) (5) r3 j 
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<2> <2,3,3.> 


t 


<3> <3,3,3,> 



6 r <r?y<v.^y 

3 (3) (5) . 


<3> <3,1,1,> 




(3) (5)R3 


<3> <3,2.2,> 


J. <rj‘><r|> 
3 (3) (5) R* 


(188) 


With these results the sums over i, j,t, s, in equation (186) are now evaluated. 

After combining the various components and (/c)^^ unit vectors, one 
obtains^^ 


L' ^ I E '‘t “ 


(-1)(4 7t)( 2!) <r,^Xr,^> 
3 R’ 


<(r>). Yf (1)> Y-f (2)> 


(1 + m)! (1 - m)! 


2| KiKi5’ii<1> <1.1. 1> +H 2«2 5>22 <2> <2.2,2> + /C 3 1 C 3 5>33 <3> <3,3.3> 

+ {K2 K 2 5^11 /Cl i?22 + ^1 *^2 ^12) 

+ («3 K3 ?ii + /C3 /C3 ^22 + 2 fci K3 ^?l3 + 2 fC2 /C3 ip 23J 

+ (/Cl /^l 5^33 + ^2 ^2 iP33 + 2 /Cl K3 ^13 + 2 fC2 K3 fP 23) ^ 3 ^ ^ 1 , 1 , 3 ^ | 


2^ 

3 


Ts) (5) _ 


(189) 


(1-3 cos^ 6 ) 


11 . 


The Kj denote unit vector components; the parentheses have been left 
out for convenience. 
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where Q is defined through (k )3 = cos Q , and Q is the angle between 

K K K 

K and the ^-axis (Fig. 2). Including this result in equation (181) and 
incorporating all previous simplifications, yields the following for equation 
(181): 


X(3) 


( 2 ) 


(2)^a^ X(0) 
(3)^ (57 t) 


(190) 


Hence, one obtains the dipole-quadrupole interaction energy due to the mixed 
terms. This equation is the counter part of equation (147); and by adding 
these two equations one can obtain the interaction energy correction, accurate 
to dipole-quadrupole orders, resulting from mixing of the field and electro- 
static operators. 


The last term in the group being considered is X(4) 
second set of terms of X(4) in equation (131): 



given by the 


X(4) 


( 2 ) 


2 (E^ - 

^-i(K+r')-R 

(fC + K*) 


^-i(K-K')-R 
[(Ej - Ej|) 4«c/c] 


/ KdK / K'dK' f dn^ f dn^.X 2 Z {*1J- I 

[(E^-E^)^hcK] [(E^-E^)^4 ^ck'] 2(Ej-E^)| 

+ <(<“2) t (Fj) ^ ^ ^ j 

1 1 

2 (Ej - E|j) J-(E^ - Eq) +JiCK'] [2(Ej -E^) +2 ^c(k + k')] 

^ ^ 1 1 

|(Ej-Eq)+*ck (Ej - E^j) +«ck' j 


(191) 
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The evaluation of the preceding factors is accomplished by following the 
procedure outlined before. The matrix elements of interest are ^ 

and These have already been listed in equation (120) . 

The nonzero combinations are: 

, ij = 1,2,3, 

<r^> , i = 1,2.3, 

<WlWl<“)2<~>2> = ■ 

<<“>1 W, W3W3 > = (i)(i) ■ 


Defining the expressions containing 1 6. - (k).(k) | by lP(i,j,f,s), one of 

I IS 1 S' 

the sums that needs to be evaluated is given by 


ijis tq J/ 


r. K<il> 


3~^ ^i,s 


+ 




Since the second factor is obtained by an interchange of r^ and r 2 , only one 
term needs to be evaluated in detail. Substituting for co^, the quantity that 

needs to be calculated (defining the matrix element as t, q, i, j^) is given by 
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(193) 


■J <^2> Z Z 6^^ 

ij is 

X y. (k. K +k'k'+2k k’) <t, q, i,j> 
f-' tq tq tq ^ ' 


By using the results given by equation (192) , the factor containing /c^ k 
becomes 


q 


(“I") Z | 3 ’(l.l.s.s) [2(ci/ci + /c.r]+ 3 >( 2 , 2 ,s,s) [2k2K2+k-k ] 

s 

+ !P( 3 , 3 ,s,s) [2 K3 K3 + ic^^] + [ !P(l,2,s,s) + !?(2,l,s,s)] [2 kj /cj] 
+ [ y(l, 3 ,s,s) + y( 3 ,l,s,s)] [2 icj Kj] 

+ [ S*( 2 , 3 ,s,s) + 3’(3,2,s,s)] [2 «2 /C3] I 

Using the definition for !P(i, j,i, s) and summing over the above result 
transforms into 


(“I-) <r2^>(-|“ )(-!■) (l + cos^e) . (194) 


The term containing /c* k* in equation (193) gives the same result; only 

t q 

is replaced by k*^. When the last factor that contains in equation (193) 

is expanded, it gives the following result: 


("i) (■j) (“s') 2 2 '‘1 *^1 ] 

+ !?( 2 , 2 ,s,s) [ 2 Kj k2 + «• k' ] + !?( 3 , 3 , s, s) [2 /c3«3+k*k'] 
+ [ ?(l, 2 ,s,s) + S’(2,l,s,s)] [ki K2 + K2k1] 

+ [ y(l, 3 ,s,s) + ^(S.l.s.s)] [/Cl /C3+ /C3 /c{] 

+ [ S'( 2 , 3 ,s,s) + 5’(3,2,s,s)] [«2 K 3 + «3 K 2 ] I • 
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By using tiie definitions for ? (i, »s) and summing over s, the above 

expression reduces to 


(t) <r|> (t )("i) ("s') 

Collecting these terms, the sums in equation (193) become 
|(k^ + k’^) (1 + cos^ 0) + 4 K k' cos^ 0 I 


(195) 


(196) 


Analyzing the expression in equation (193) indicates that when 




is replaced 


by Vj 


only the sign of 2 k' 
t q 


changes; hence, only the quantity in equation 


(195) changes sign. By using these results to evaluate the second term in the 
sums of equation (191) and performing integrations as before as well as using 
the results of equation (176), equation (194) reduces to 


X(4) 


( 2 ) 


4 a? X(0) 
(3)^ (5) 7T^ 


/ bdb J /3d/3 


I (2a + b) / 1 1 \ 

I a(a + b)^ \ b + /3 ^/3-b/ 

a(b + j3) + j8(a + jS) + b(a + b) 
a(a+b)(a + ^) (b + ^)(b-/3) 


xj (b2 + /s2) [F(b) FO) + 2 G(b) GO) ] 

+ 6 F,(b) F^W) 1 
^ (biS) 


(197) 
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The /3 integrals are evaluated using previous results and values 
obtained in Appendix C. Since the second fraction containing b and j3 does 
not contribute to the integral over /3, only the first term, needs to be evaluated. 
To obtain the /3 integrals in equation (197), one proceeds as follows: Take 
the first term, given by 


I [F(b) F(«+2G(b) Gtf)] 

= [F(b).2G(b,] 


+ 2i7T f [F(b) F^(b) + 2 G(b) G^(b)] 


(198) 


and the second integral of equation (197) , given by 



combine these terms and obtain the following expression for equation (197): 
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X(4) 


( 2 ) 


4 a® X(0) 
(3)^ (5) tt" 


X 


oo 


3 J 


b(2a+ b) 
a(a + b)^ 


db 


F(b) + 2 G(b) + 52(b) - F3(b) 


36 


OO 

+ 61 / 

0 


b^(2a + b ) 
a(a + b)^ 


db 


I F(b) F'^(b) + 2 G(b) g"^ ( b) - 
6 Fa (b) F^ (b) I 

^ j 


3^2(b)5?~(b) 

( 200 ) 


( 2 ) 

Further simplification is possible for each of the X(i) ^ terms, by applying 
the definitions for the various functions used so far. This may be seen by 
recalling that these functions are defined in terms of F(b) and G(b) , used 
throughout this discussion. For instance =t>^ G(b) - 2 F 3 (b) , 

5(b) = b^ G(b) - F 3 (b) , ... etc. 

Dipole-Quadrupole Interaction Energy 

( 2 ) 

In the previous discussion, each of the elements X(j) ' is given in 
terms of an integral over b (b = kR) and with coefficients proportional to 
( 2 ) 

X(0) . Combining the expressions in equations (168), (172), (180), (190), 

and (200), one obtains the dipole-quadrupole order correction to the interaction 
energy AE^ , given by 
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AE , ^ = X(0) 
d-q 


( 2 ) 


+ 8 i / b*db {(FF'^ + 2GG^) -p ( 55 ^ + 4 F 3 F^) |j 


[ 2 a s 

bdb 1 

[ 5 (F + 2 G) + [ 12 5 + 36 5 2 

- ( 3 ) (52) F3 ] ) 

(a + b) 1 

+ 4aiJ 

b^db 
(a+ b) 

{ 5 (F f'*’ + 2 G G^) (255 

■ + 3?3?2 +26 F3F3) )j 




[ 3 “'/ I (y2-3F3)j 


{ <r'F‘"+2GG'") (y2?2 +6 FjFt) 


( 201 ) 


By writing the above equation^^ as a sum of terms, one can see that the result 
reduces to the electrostatic interaction case when the radiation field is omitted. 

( 2 ) 

In addition, each term X(i) , i = 1, 2, 3, 4, is separated into two integrals 

over exponentials e"^^ ^ and e*^ ^ The reason for doing this may be seen by 
referring to the dipole- dipole approximation, where, inequations (162) and 

(163) , one notes that only the integrals over e*^^^ contribute to the final 

ii b 

result given in equation (164). The integrals over e and portions of the 

results from the integrals over combine to cancel X(0)^^\ One 

would expect the same behavior for higher approximations, so the results 
given by equation (200) are already in a suitable form for further analysis. 


Q UADR U POLE-QU A DRU POLE A P PR OX I MAT I ONS 


Introduction 

The presentation of the quadrupole-quadrupole order approximations 
is much more condensed than previous discussions; however, all the steps 


12. In the above expression the explicit functional dependence is omitted. 
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necessary to obtain a given result are either mentioned explicitly or the inter- 
mediate results are given. To have a systematic presentation, the procedure 
followed here parallels that used in the immediately preceding discussion. 

The quadrupole-quadrupole order interaction energy is given by 


AE = X(0) + X(l) + X(2) + X(3) + X(4) , (202) 

where the various X(j) j > 0 , elements are given by the last group of 
terms in equations (123), (127), (129), and (131), The electrostatic 

(3) 

interaction energy X(0) ' is obtained from equation (85) by letting 
Li = 2; M = -2, -1, 0, +1, +2; and Lj = 1, L 2 = 3; Lj = 3, L 2 = 1; 

M = -1, 0, +1. Hence substituting these values in equation (85) , one obtains 


X(0) 


(3) 


(-1)'^ 1(4)!} VriO<r,0 
2(E^ - E^)Ri»(5)(5) 


+2 


Z 

_ m=-2 


1 

(2 - m) ! (2 + m) ! (2 - m) ! (2 + m) ! 


+ 


(-!)« { (4)1 } (-1)^ {(4)!} ^ <r^g><'r,^> 

2(Eq - E^)R“(3) (7) 2(E^ - E^)r1«(7)(3) 


^ ^ (1 - m) ! (1 + m) ! (3 - m) ! (3 + m) ! j 

= (E^ - e|^)r1«' { 'T + T <^1 ><^'> j 

This result may be expressed in terms of using the following 

relations, [see equation (86)]: 


<r«> = <r2><rO . 


<r^ 


( 2 ) 


<r2 


Xr2> 


Thus, X(0) 


(3) 


may be written as 


X(0) 


(3) 


(5) (7) <riO<r,0 
32 (E^ - Eq)R^® 


(203) 
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The next term to be evaluated is given by the last group of elements 
inequation (123): 

^ [ (t) (t) ' <■•■'> <‘'^'> (m ) (t ) (<■■*> * <‘'*>) . 

X I (k^ + k'^) + 4 k k’ (k^ + k’^) cos 0 + 4 k'^ cos^ 0 } . (204) 

The coefficient of this term may be expressed in terms of ^ by 

using the relations previously obtained. Combining the r matrix elements, 
and substituting for (l/2pc^), one obtains 


(2)=*(3)^(7r)*(5) (^ic)3 


Again substituting for k and k* and using the above coefficient, as well as 

(3) 

expressing the quantities in terms of a and X (0) ' , equation (204) becomes 


X(i) 


(3) 


2^ a? X(0) 

(7) (3)2(5) V 




(b + /3) 


X j (b2 + /?2) [r(b) FW) + 2 G(b) G(/3)J 

_2b/3(b2 + /32) IZlhIJM^4F,(b)_F,_(^_)] 

+ b2^2 1^ {F^(b) F 2 O) + 2 G(b) G(^)) 


+ -]^{f(b) f(/3) + 12 ?4(b) ?4(^) + 24Fg(b) F3(/3)}j 


(205) 
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where 


F 2 (b) = F(b) + G(b), f(b) = b2 F 4 (b) - 8 F 3 (b), = b^ G(b) - 4 F 3 (b) . 

The /3 integrals are performed with the aid of Appendix C . The results for 
each of the above integrals are as follows, (after the (3 integrations, the 
explicit functional dependence is omitted) : 

J bdb / /3d/3 ^ (b2 + ;S2)^ [F(b) F(/3) + 2 G(b) G(/3)] 

oo oo 

= 7T J b3db{F + 2 G} + 4i7T J b®db { F f'^ + 2 G g“^ } 

0 0 

= - 3 7T j bdb I?- 4 Fgl + iTT J b^db 1 + 4 F 3 fJ j 

0 0 

f bdb f { [F 2 (b) F 2 (j 8 ) + 2 G(b) G(^)] 

+ -^ [f(b) f(/3) + 12 54 (b) + 24 Fg(b) F3(^)] I 

OO OO / 

= - 24 IT J { f + 6 ?4 - 3 F 3 } + iTT J b2 db j b^ F 2 F^ + 2 b^ G G^ 

0 0 

+ ff +125454 + 24 F 3 F3 I , 

Using the above integrals, the expression in equation (205) becomes 
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X(l) 


( 3 ) 


X(0) 

(7) (3)2(5)2 7T 


x| J db|^b3(F + 2 G) + 12b(5- 4 Fg) - ^ (f + 6 J 4 - 3 Fg)] 

00 

+ 4i J db 
0 

b® + + 

+ — (Fg Fg + 2 G G ) 


b® (F f"*" + 2 G G*^) - b^ ( + 4 F; 


Z^t) 



(f f + 12 ?4 + 24 Fg Fg ) 




(206) 


In the above expression, the explicit functional dependence in the various 
functions has been omitted. This will be done in the future after doing the yS 

( 2 ) 

integrations. Comparison of equation (206) to the expression for X(l) ' ' 
given in equation (172) shows the increased complexity in the higher approxi- 
mations, even for X(l) which is the simplest of the electromagnetic 
interaction terms. 


The next term to be considered is obtained from the last group of 
elements in equation (127); it is given by 


X(2) 


(3) 


2 (E^ - E^)2 




^-i(/c+K')-R ^+i(K+K')-R 

(k + «') [(E^ - E^) +ficK] (k + k’) j^(E^ - E^) +ficK'J 


T E E { (i) (t) “t ‘"s ‘"q “h 

L ts qh I X 'i ^ 


(t) "t "s "q “h <(^2)t(r2)^(r2)q(73)j^(72) . (^3) . ^ I 


114 



+ 


e ' 

[(E^ -E^) +fic/c ][(E^ -E^) +iicK«] 


X 



I {(t) (t) w w > 


(■r)’'t ''s ''q ''h <^<^>,<~2)sS>q‘^>h<“‘>i<'>>j> ) 


(207) 


To evaluate the terms in equation (207) , one needs to consider three 
different types of sums, as given in the above expression. Since the evaluation 
of these sums is very involved, only the intermediate results are given in each 
case. Denoting the sums over unit vectors (k). by ^(i.j), the first sum 
to be considered is 



3’(iJ) S Z 

ts qh 


6. . CO CO CO, 

1 ] t s q h 




Using the matrix values for 




obtained before, the 


expansion over t, s, q, and h is simplified considerably. Furthermore, the 
function allows one to collapse the sum over i, j. After this is done, the 

above sum reduces to 


cos 0 


(^ )(^^) ("i) ("5 («' + 

+ 4 k'^ cos^ 0 j (1 + cos^ 9) . (208) 


The next sum to be evaluated is 


' * ij ts q^^ 
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Summing over _t, first, and then over i, j, q, h, the above expression 
becomes 


(4) (3) 




<r}> 


(3) (5) 


+ k’^) (1 + cos^ 0) + 4 K k' COS^ gJ 


= ^ 3 ) ^ 5 ) + cos^ 0) + 2 K /c’ (k^ + (cos 0+3 cos® 0) 


+ 8 /c® k'® cos^ gJ 


(209) 


Evaluation of the last sum in equation (207), given by 


ij ts qh M N n j / 

is carried out by first replacing o) and w’ with k and k' in the above 
expression, and then performing the multiplications indicated. Five different 
sets of (k) ('f') j combinations are obtained; for instance, one of these sets 

having no primed components is (k^ ^q ^h^ ' other groups have various 

combinations of (k ) and (k ’) . The matrix elements over (r). are 

t t X 

evaluated in the same manner as those for dipole-quadrupole approximations; 
the nonzero combinations are as follows: 

<((^)i(^)iU).(^)j(^)j(^)i)> = , i = 1, 2, 3 , 

, i.j = 1. 2, 3; j . 

<(WiW.(V) = TshIW’ j.‘. = *. 2. 3i i’* J’* t • 

(210) 
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Using these results, the term being considered reduces to 


(“^) Z Z Z <t,s,q,h,i,j> 

ij ts qh 

X I K. /C K K, + k[ k' k' k' + 6 K k' K* + 4 K K k! + 4 K* K* K,* I 
[tsqh tsqh tsqh t sqh tsqh) 

< r® > I 2 

= (8) (3) (5) (7) I 0) + 4 (-cos^ © + 3 cos^ 0) 

+ 8 K k' (k^ + k’^) cos® 0 I 


( 211 ) 


Analysis of equation (211) shows that when the sums containing (v^) 

in equation (207) are evaluated, a sign change occurs in those terms containing 
cos © or cos® 0. Thus, using previous results, one notes that the two groups 
of terms in equation (207) combine after the integrations over do are 
performed because of a sign change introduced in these terms in these operations. 
[See equation (176). ] Thus, the equivalent of equation (207) may be written as 


X(2) 


(3) 


2 (E^ - E^)2 
2^ tT* uc^ii 


J KdK J K'd«' J J" e 


■i(K+K') • R 


(K + K') [(E^ - E^') +*ckJ (k + k') [(E^ - E^j) +ftOK'] 


c 


||(E^ - Ejj) +«c»c] [(E^ - Eq) +*ck'] 


<r,^Xr,«> 


(24) (3) (5) 


I [ (k^ + + 4 k k' (k^ + k’^) cos e + 4 k' k’^ cos^ e j X (1 + cos^ e) I 


('^'1 

'<r,2><r^> 1 

w / 

(3) (3) (5) 


[(k' + k'") (1 


+ 8 k'^ cos* ej 


<r«> 

(8) (3) (5) (7) 


[(k^ + k'^) (1 + COS^ 0) + 4 k’^ (-COS^ 0+3 cos^ 0) 

+ 8 K k' (k* + k'^) COS^ 0J I , 
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where one combines the terms in equation (207) before the integrations over 
dJ2 are performed. An additional simplification to equation (207) is accom- 
plished by expressing (ju) in terms of and the resulting factors in 

terms of ^r/^^r 2 ^^. After performing these modifications, expressing k 
and k' in terms of b and /3, and listing the coefficient in terms of a 

(3) 

and X(0) ' , the preceding equation becomes 


tm5)U7) I I I 


a^ X(0) 


-i(/c+/c‘) • R 


2*7T* (3) (5)^(7) 


/ ^ / 1 ^ 1 \ 
|b + /3\a+b ^ Si+ p J 


• I ■ ■ ■ “ 

a + 0 / (a + b) (a + 0) 


X I 5 (b^ + 02)^ (1 + cos^ 0) + 8b0(b^ + 0^) (cos 9 + 4 cos^ 6) 
- 4 b^0^ (cos^ 0-11 cos^ 0) j 


( 212 ) 


By performing the integrations over dfi the above expression simplifies to 


(3) (7) (5)^ TT^ J J (a + b) 


(3) (7) (5) 

x| 5(b2 + /3^)^ [F(b) F(i3)+2G(b) G(/3)] 
- 4b/3(b^ + /3^> 


1 ^ 1 ^ b(a + b) + /3(a + fi) 


(b + 13) (/3-b) (b-/3)(b + /3)(a + /3) 


bG(b)^G(^) .4 

u p 


- 2 b^^^ 


F2(b) F 2 W) + 2 G(b) G(^) - 11 


(213) 

[f(b) f(g) + 12 7j(b) 7^(13) + 24 Fa(b) F.,«3)] 

1 , 2^2 


Since the second fraction involving b and 0 does not contribute to the 0 
integral, one needs only consider the following 0 integrals: 
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/ + [F(b) F(^) + 2 G(b) G(P)] 


= ^j <^ + 2G) +41. 7 {fF-*-+2GgM . 

7 bdb 7 -)(rir^7^ -j [b/3(b2 + /3=“)]x I bG(b) fiG(fi) 

= - <8)(3.) 7 {y2-3F3} +2.1 7 I Gg"' + -^^ 


tfdb + 4?,yr , 24F,F3 

[a + b) * b* b* 


eo oo . 

/ bdb f pdp [-^] j bV [F2(b) F3(^) +2 G(b) G(^)] 

- 11 [f(b) UP) + 12 y4(b) yi(/3) + 24 F3(b) Fj(^)] j 

oo OO A 

= (2)(11)(12). / j f+6y,-3F3} +1. / j(F3Ft + 2GG^ ) 

- ^ (£ r + 12 ?4 y; + 24 f, f, > 

Incorporating these results into equation (213), one gets 

V£91<3) _ (2)^ a^X(O) 

^ ^ (3) (5)^7) (F) 


J j 5b3 (F + 2 G) + (4)(8)(3)b (52 - 3 F 3 ) 


(2) (2) (11) (12) 


(f + 6?4 - 3 F 3 ) 


+ i / -(Tfly i (5) (4) (F + 2 G G'*’) 


- (4) ( 2 ) I b“= G G + 


2 ^[f 2 fJ + 2 G G"*"] 


+ 4 5, 5, 24 F, F, 


f f + 12 54 54 + 24 Fg Fj 
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nil mil 


Using the definitions for the various functions, one notes that the preceding 

result is separated into two terms proportional to and as done 

previously. 


(3) 

The next term to be evaluated is X(3) . This term consists of the 

last group of elements in equation (129) and the terms of quadrupole- 
quadrupole order in the second group of this equation. [See equations (182) 
and (183) . ] The reasons for including all these terms here were explained 

(3) 

in the preceding discussion. Collecting these terms, X(3) ' ’ is given by 


X(3) 


(3) 


4 (E^ - E^)2 
2“ir^4ic 


/ /<d/c J Z {^ij - 


-i/c- R 


X 


X 


1 <^l-V [(eJ- ) 



Z Z «h 

ts Ih 

) ■ ,2,5. 

Expressing the elements involving ( k ) ^ as , the sums needed in 
the above equation are evaluated as follows: The first term to be considered 

IS 

gives 


13. Using equation (183) one picks out only the terms for which = L 2 = 2. 
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2 % 

ij ts 

= 2 ^Ij Z <-.V<r.‘> 


ij ts 
+2 


2 [( 2 - m )!( 2 . m ) T ] 


m=-2 


{ ^ - ^0 COS^ H. 35 cos* 9^ . 


where the following nonzero matrix elements have been used: 


(216) 


<^h^^'(7)j(7).(7).Wi)>= (^) ■ 


(r) . (r) . (r) . (r) . 
\ q 1 'i 1 1 



i= 1,2 


i = 3 


<h<2)(7)j(7)j(7,,Wj>= . u = i,2. i#i 

<h<2>W,WjWj(T)j)>= . i,i = i,2, i^j 


= 1,3, i# j , 


<h®Wj(T)jW^(Dj>= (-i|)(^-^^^) . M 
<h<2'(T)iWjWjWj^= (-n)(^^^R^) ■ = ‘'‘J • 
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<CWiWiWjW)>= (- 


4 \/±<Il>^ 


, i,j = 2,3, i^j 


, i,j = 2,3, i#j , 


The remainder of the matrix elements are obtained by interchanging t and i 
in (r^) and j and s in (r 2 ) . Recall that the evaluation of each element in 
equation (217) follows the same path as that followed to obtain equation (187) . 

The next sum to be evaluated gives 


ij ts ^ ^ •’ 


- y S’ 7 K K r <r^><rj’> 

U « te ‘ “[ [(3) (7)] 

+1 < (l)><(f2)j Yp* (2)> 

' Ij ; — 

m=-l I (1 + m)! (1 - m)! (3 + m)! (3 - m)! } 


The above results are obtained by taking Lj = 3, L 2 = 1 in the expression for 

( 2 ) 

. The choice of L values is made using the results of equation (185). 

The various matrix products needed to evaluate equation (218) , [defining 
quantities as in equation (188)] are: 


< 1 , 1 , 1 > < 1 > 


(5) (7) R® 


= < 2 , 2 , 2 > < 2 > 


< 1 , 2 , 2 > < 1 > = - — 


3 (5) (7) R“ 


= < 1 , 2 , 2 > < 2 > 


<1,3,3> <1> = 


3 (5) (7) R® 


= <2,3,3> <2> 
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<3,3,3 > <3> 
<1,1, 3> <3> 


_8. <T^> <r,^ > 

3 (5) (7)R® 

± Xr,0 

3 (5) (7) W 


= <2,2,3> <3> 


(219) 


Using these results equation (218) becomes 


- E, *'ii S 

1,J ^ t,s ^ •’/ 


Jl. 

3 


<v^ > <r|> 


(3) (7) rS 


(1-36 cos^ 0+23 cos^ 0 ) 

K K 


° - - iM0 ~ { ^ °°°' ^ } ' 


( 220 ) 


where the last step is obtained by substituting for <^r® > and <r^ > in terms 
of <r^>. 


The next sum to be considered is given by 
Z \ (-t) Z Z «,«s“q«h 


Ij 


ts qh 


I f.. E Z «,« « «t Z ^ 

ij ts qh ^ ^ Li, L 2 2 { (2 Lj + 1) (2 Lj + 1) } 


E 

m 


<((ri)t(ri)^(r ,) . (1)^ <C^^2)g(r2);^(r2) . (2)^ 


I (Lj + m) ! (Lj - m) I (Lz - m) ! (Lj + m) ! | 


I I /a 


( 221 ) 


The matrix elements required to evaluate the above sum have been evaluated 
and are listed in equation (220). The quadrupole-quadrupole order terms are 
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obtained by choosing only the terms for which Lj = L 2 = 1 in equation (221). 
Letting 


(^i>t<^i>s^^i^i^^ 2 )i(^ 2 )h(^ 2 )j^ = <t,s,i><l,h,j > 


the nonzero matrix products corresponding to equation (221) are as follows: 


< 1 , 1 , 1 > < 1 , 1 , 1 > = 


V 5 / \ (3) (5) r2 / 


<^ 2 , 2 , 2 > < 2 , 2 , 2 > 


<1, 1, 1 > <1,2,2 > ( 5 ) (^ (3) (5) r2 } > <1.3, 3> , 

<1,2,2> <1,2,2> = = <1,2,2> <1,3,3> , 


< 1 , 1 , 3 > < 1 , 1 , 3 > 
< 2 , 2 , 2 > < 1 , 1 , 2 > 
<1,3,3> <1,3,3> 
<3,3,3> <3,3,3> 
<3,3,3> <1.1, 3> 
<2,2,3> <2,2, 3> 
< 1 , 1 , 2 > < 1 , 1 . 2 > 




= <1,1, 3> <2,2,3> 

2,2> <2,3,3> 

= <2,3,3> <2,3,3> 


= (-!)(- = <3.3,3><2.2.3> . 


<1,1, 2> <2,3,3> 
( 222 ) 
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Using these results, the sums in equation (221) are evaluated and the results 
are given by 



ts Ih ^ 


(4) (3)2(5)^ 


( 


1-3 cos^ 


e 

K 


) 


(223) 


This result is in itself interesting, in that after numerous operations on the 
terms in the sums, the final result can be expressed in terms of 

(fc) 3 (fc) 3 = cos^ 0 . In the evaluation of X(3) the angle 6 plays the 

K ' K 

same role as 0 in the terms having both k and k '. 

The next sum to be evaluated is 

IJ ts Ih ^ ^ 


The sum obtained from this expression by replacing (r^) 
same value and need not be evaluated. Substituting for H 
expression, one finds that matrix elements of the form 


and 

( 2 ) 

q 


(r 2 ) has the 
in this 


A)jA)jY“Ni)^ and Y-“"(2)^ 


need to be evaluated. The former group has already been evaluated and is 
listed in equation (183) . This group of terms requires that = 2 for 
nonzero results. Thus, the quadrupole-quadrupole order terms associated 

with this sum require that L 2 = 0 to get results proportional to ^i^<r 2 ^^ . 

With these restrictions, one notes that sums of this type do not contribute to 
the quadrupole-quadrupole order results. 
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The last two sums that need to be evaluated are of the form 

” tS Ih n N -a 


( 2 ) 


being obtained by interchanging (r^) and (r 2 ) . Substituting for H in 


terms of rj and r 2 , one notes that matrix elements of the form 

<^(^i)i y“*( 1)^ and <^(72)^(72)^(72)^(72 )j^(72). Yl“'‘'( 2)^ need to be 

considered. The former set is given by equation (143) , where one sees that 
only the value = 1 gives nonzero results. To quadrupole-quadrupole orders, 
then, L 2 can be only unity. Hence, this sum reduces to 


(^) S 

^ ' ij 


^ ts Ih 


t s 





E E E 


ij 


ts qh 


K. K K K, 

t s q h 


+1 

X E 

m=-l 


Yf" (l)^<^A),(rj)^A)^(r,)^^(r,) . Y~”’ (2)^ 


(1 - m)! (1 + m)! 


(224) 


The matrix elements in this equation are evaluated in the same manner as 
before; except now, the substitution of the (r) ^ components requires the use 

of up to fifth-order Spherical Harmonics. These functions may be obtained 

using the definitions for Y^ (0,0) given in equation (21) . The nonzero 

matrix products corresponding to equation (224) are as follows: 


< 1 > < 1 , 1 , 1 , 1 , 1 > 
< 1 > < 1 , 1 , 1 , 2 , 2 > 



< 2 > < 2 , 2 , 2 , 2 , 2 > 

< 1 > < 1 , 1 , 1 , 3 , 3 > 
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< 1 > < 1 , 2 , 2 , 2 , 2 > 


<1> <1,2, 2,3, 3> 
<2 > < 1 , 1 , 1 , 1 , 2 > 
< 2 > < 1 , 1 , 2 , 2 , 2 > 
<3> <3,3,3,3,3> 
<3> <1,1,1,1,3> 
<3> <1,1,3,3,3> 
<3> <1,1,2,2,3> 


/A /<r?><ri> ) . 

\ 5 / \ (3) (7) R» I 

( jJ) ( 

\ 5 / \ (3) (7) R® / 


<1> <1,3, 3,3, 3> 


= <2> <1,1, 2,3, 3> 


= <2> <2,3,3.3,3> 


( J_\ ( \ ^ 

V 5 / \ (3) (7) bF J 


< 2 > < 2 , 2 , 2 , 3 , 3 > 


(- 2 ) 


(3) (7) r3 j 


( 2\( <r,^y<rj> \ 

\ 5 J\ (3) (7) R“ / 


'<r?y<rn \ 

(3) (7) j 


h: 

( (- — / <^i><^2> \ 

\ 3 I \ 5 ; (3) (7) r3 j 


<3> <2,2,2,2,3> 
<3> <2,2,3,3,3> 


(225) 


Using these results, the sums in equation (224) combine to give 
(^) 2 ’'ij 2 E 

' ' ij ts Ih \ J /- 


where the last results is obtained by replacing <r® ^ <r^ in terms of 
<r‘*^<r‘*^. The sum obtained by interchanging (r^) and (rj) inequation 
(224) is the same as in equation (226). Substituting these results in 
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equation (215) , letting kR = b, and factoring common factors, this equation 
becomes 




Using the results of Appendix B to perform the integrations over dS2, the 
above integrals reduce to 

f dn [13 - 234 cos^ 0 + 197 cos^ 0 ] 

= f dn [l3(l - 3 cos^ 0^)- 195 (cos^ 0 - cos^ 0) + 2 cos^ 0 ] 

^ K ^ K K K 

(13) (-2^ (F3) + (195) (2^) + (2) 
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f dS2 Q~it>cos ^ 2<7 ^^^2 0 _ 43 cos^ 6 ] 
j K K K 

= I (7) (-2^it) (Fg) - (48) (2M -||- j 
Using the above results, equation (227) becomes 


X(3) 


(3) 


a^ X(0) r bMb 

' (7) (3) (5)“7T J a(a+ b)^ 

X I (3) (2)2 b2 ^ (13) Fg - 195 


+ 


b^ 7 F, + 



(228) 


The last term to be evaluated is given by the last group of terms in 
equation (218) , and is given here by 


X(4) 


(3) 


2 (E^ - 


f KdK f K*dK* J f dl^^, 


ij is 
e-i(K+K’) 

■ftc (k + i 


b-R ( ^ 2 (Ej - Ep) +lic (/c + k') I 

( [(E,- ^ '[<^1 - ^ fic']Y(Ej“E„)+«CK'][2(E^-E„)] j 

>< { n 2 E “t “q “h "f [ s} 

+ <^(Fj) .(Fj) <(F2)^(F2)^(F2)^(F2)^(F2)^ (7,)^y j 

^-kF-Fv-r 1 1 

[ (E^- Ej|) + Kck] I 2 (Ej - Eg)[(E^ - E^) + Kok’] 


2 (Ej -Ejj) +ftc (k + k') 


[2{E^ -Ejj) +^c (K + «•)] [(E^ - Ep) +«ck] [(E^ -E^) +J1ck’] 
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X 


+ <^ (^t) i (^l) j t ^^2) q (i^) ji (1^2) f (^z) i (1^2) ^y j 


(229) 


The evaluation of this equation requires evaluation of six different sums 
involving co^ and v^. Later on, one will see that it is necessary to do only 

the sums involving co^. Taking the first sum, given by 




tqhf 


one evaluates it by using the results of equations (120) and (210), where the 
various matrix elements needed have already been listed. The results are as 
follows; 


{~^} (“s') (“ 3 ')(“ 5 ‘ )('f ) 1 ^ ■*■ 

+ 6 ( k ^ k '^) (1 - cos^ 0 + 6 cos^ 0) 


+ (4) (3) (2) (k® k' + k'^ k) cos® 0 } . 


(230) 


The listing of terms reflects the expansion of 


(co. o) CO, CO J into various 
tqhf 


(k)^ and (k*)^ components, as was explicitly indicated in equation (211) . The 

sum obtained by interchanging (r^) and (r 2 ) yields the same results and need 
not be evaluated. The next sum to be evaluated is given by 


E Z (i) Z 2 1”® )> • 

i,ji,s t,qhf^ L \ ^ J/N /J 


The evaluation of this term is accomplished using the results given in equation 
(175) for the matrix elements. Since t, q, h, and f appear in different 
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I 

/ 


matrix elements, one must be careful when performing the expansions and 
contractions of the various (k)^ and (k')^ components. In addition, since 

i, j,i, and s are coupled in both matrix elements and the polarization function 
?(i, j,i, s) , the sums are not straightforward. Applying the various techniques 
developed in evaluating previous sums of this type, one finally obtains the 
following result for this sum: 


(■^ j (■f ) ("s) l"3’) ("s') 

+ 2 k'^ (1 + cos^ 0) 

+ 4 (1 - 3 cos^ e + 4 cos^ 0) 

+ 4 (k®/c' + k'^k) (2 cos^ 0) } 


(231) 

Analysis of the results in equations (230) and (231) shows that, when the next 
group of sums in equation (229), containing (v ) instead of (w ) is evaluated, 

I L 

the elements undergoing a sign change correspond to terms containing cos^ 0. 
This property allows one to combine the two groups of terms in equation (229) 
because the integrations over dS2 of cos 0 and cos® 0 give an additional sign 
change which offsets the previous variation in sign of the terms in question. It 
is of interest to note that while this feature of the calculations appears throughout 
the various approximations considered, it is necessary to keep the various terms 
separated until the final sums and integrations over dfi are carried out. In 
all previous cases this has been done, except in equation (212), where the 
various terms in equation (207) were combined before performing the integrations 
over dl2. This will also be done here to shorten the discussion. Hence, 
incorporating the results given in equations (230) and (231), and combining 
the various elements in equation (229) , one obtains 
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X(4) 


(3) 


2 (E^ - 




^-i(K+/c')«R 


( _R® _ 

1 2a+(b + 0) 1 

1 («c)*(b + |3) 

(a+b)* 2a(a + b) (a + )3) J 


R-* 


(*c)''(a+ b) 


' ^ 

2a(a + 


2_^ 

b) ■*■ (2a + b + 


+ (b + g) 1 ) 

/j) (a + b) (a + |3) J | 


I + ^^)(1 + cos* e + 2b*^*)(l - cos* 0+6 cos^ 0) 

+ 8 (b*/3 +/3*b) cos* 0 ] 


+ 


(4) (3)* (5)^ R* 


[(b^ +;8^) (1 + cos* 0) 


+ 2b*/3* (1 + cos* 0>i-4b*^* (1-3 cos* 0 + 4 cos^ 0) 

+ 8 (b*/3 + j3*b) cos* ©j] | • 

' (232) 


Expressing and in terms of r^ )> and rearranging terms, 

equation (232) becomes 


X(4) 


( 3 ) 


-2 (E^ -E^)* 
2^7 t^ (iic)®R*’ 


<r^><TrO 

( 4 ) ( 3)3 ( 5 )^ 


/ bdb / 


( (2a+ b) / 1 1 \ a(b + jS) +g (a+ ^) + b(a+ b) 

^ \ a(a+ b)2 V'b + i3 |3 - b / a(a + b) (a + /3) (b + ^) (b - /3) 

X I 7 (b^ + ^^) (1 + cos^ 0) 


+ 2b3^3 _ j^9 (jos^ 0 + 48 cos^ 0) 
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+ 8 (b^/S +P^h) (7 cos® 0) 



(3) 

Expressing the coefficient in terms of a and X(0) , and performing the 

integrations over dJ2 , the preceding expression becomes 



a5x(0)^^^ 

( 7 ) ( 3 ) 


X J bdb J /3d/3 


(2 a + b) 
a(a + b)^ 



+ 


1 

^ - 


b 


) 


a(b + /3) + /3(a + j3) + b(a + b) 
a(a + b) (a+ (b + /3) (b 


7(b^ + /34) [F(b) F(i3) + 2 G(b) G(^) ] 


+ 2 b^ 


13 


{F(b) - G(b)} {f(|3) - G(/3) } 


{^2(b) F2(/3) + 2 G(b) G(/3) } 


48 |f(b) f(/3) + 12 3^, (b) 54 (^ 8 ) + 24 Fs(b) F,(^)} 

2 b^/3^ 


- 8 (b®/3 + /3®b) 


7. [y.(b) y,()3) + 6 Fa(b) Fa(/3)] 
2 b/3 


( 233 ) 


The integrations over jS are executed using the following results: 
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/ ii? - ; [p(b)F(»»2G(b)G(e)] 

/ “ly^y / -ijrwfrb) <'>’'^' [(F(b)-G(b)Hp((i)-G(f)i] 


/ b-j^ I wX(e-n [F,(b)F,w*2G,b)Gm] 


(|3 + b) (/3 - b) 


f ^i<f^±^{F,Ft+2GGM 

J a(a+b)2 I z z 


r b(2 a + b)db r ^ , . . 

J a(a+ b)^ J (p+b)(,p-b) ^ ^ L 


2gMg 




f(b) f(/3) + 12 y,(b) + 24 F^(b) F,(g) 


= - { f + 6 ?4 - 3 F3 } ^ i^ J - JfcVy {ff-+12J4?r+24F3F:} , 


a(a + b)^ 


r b(2a + b)db r 2gMg 3 I J ^ 1 

J a(a + b)=‘ J (p + b) (^ - b) ^ ^ ‘'H b^ J 


= / *^f(a+ b!^^ 1 ? 2 - 3 Fs } +2i5T J {3^2 3^2 + 6 F 3 F 3 } 


^ a(a+br 


Using the above integrals, equation (233) becomes 


X{4) 


(3) 


a° X(0) 
("7)(3)(5)";r 


+ (4M7}(6) b (y^-3F3) j 

' f (F + 2 G + 13 (F - GJ (F"^ - g"^) 

- (19) (Fj Fz +2GG‘^} - - (? 2?2 + 6 F 3 F 3 ) 


+ (f r + 12 3=4 3 =; + 24 Fg 


•’ll 


(234) 
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I 

i 


(3) 

The preceding result for X(4)' may be expressed in terms of polynomials in^ 

and exponential functions and by using the definitions for F(b), 

G(b), . . . f(b), given throughout the discussion. If this is done, the resulting 
expression would cover several pages; therefore, this will not be done here. 


Quadrupole-Quadrupole Interaction Energy 

(3) 

Having evaluated each of the X(j) quantities in equation (202), the 
expression for the quadrupole-quadrupole order correction to the interaction 
energy may be written down using the quantities defined in equations (203) , (206) , 
(214) , (228) , and (234) , By referring to each one of these equations, one can 

(3) 

see that each of the terms X(j) , j 0, has been given as sums of two terms 
proportional to and and that each coefficient has been expressed 


in terms of X(0) 
(202) becomes 


(3) 


Hence, substituting for the various X(j) 


(3) 


equation 


AE = X(0) 

q-q 


( 3 ) 


1 + 


(3) (5)’7ir 


/ db I b’ (F + 2 G) + 12 b ( 7 - 4 Fs) - (f + 6 7 ^ - 3 Fj) I 


4i J db b® (F F'*' + 2 G G"^) + "Y (F2 Fj + 2 G ) 


b‘ ( + 4 Fa Fj ) + -y (£ £ +12 7^ 7^ + 24 F, F3 ) 


J (F + 2 G) + (4) (8) (3) b (7a - 3 F3) 

- (2)(2MUMt2> (t, 6 7, -3 Fa)) 


+ i / {(5)(4)b«(FF*+2GG'^) -2b‘(FaFj+2GG'^) 


+ 


J a"(a + b) 2 


- (4) (2) b^ (b^ G g'*' + 4 3^2 + 24 F 3 F3 ) 

+ (2) (11) b2 (f f" + 12 ^4 3^7 + 24 F 3 F^)} 

1 12 b^[l3 Fj - (13H15) - -1^ j + b* (7 Fj + I 
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/ ( 7bMF ^ 2 G) - (t ^ 6 ?4 - 3 Fa) 

+ (4) (7) (6) b (?2-3 Fs)| 


‘ / ' ^f^a ' +^bF ^ b' (F f'" + 2 G G*) 


13 b® (F - G) (f’*' - o’*") - 19 b® (Fj fJ + 2 G G'^) 
- 56 b® {?2 72" + 6 Fj Fj) 

+ 48 b® (1 r + 12 ?; + 24 F 3 Fj ) } 


(235) 


In the above expression for AE , various groups of functions appear in more 

than one place and each group of terms consists of two integral expressions 

proportional to and e^^^^ as before. Comparison of equations (201) 

and (235) shows explicitly the increased complexity encountered in going from 
the dipole-quadrupole to quadrupole-quadrupole approximations. In particular, 


the group of terms corresponding to X(3) 

( 2 ) 

corresponding group in X(3) 

r(2) 


(3) 

is much larger than the 
This shows that if one is interested in the 


coupling between H and the field operators for approximations beyond the 
Q. 

dipole -dipole, one must go on to quadrupole-quadrupole orders before obtaining 
additional results of interest. 


ADDITIONAL RESULTS AND CONCLUSIONS 
General Remarks 

In this, calc Illation the interaction energy of a two-atom system has been 
obtained to quadrupole-quadrupole orders. The results are given in terms of 
the various approximations, to compare the first order results with those of 
Casimir and Polder [1]. Subsequent approximations are expressed in terms of 

the electrostatic interaction energy X(0) which in turn gives self- 
consistent results in the various degrees of approximation. This procedure 
also serves as a check on the coefficients associated with the different cases 
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considered. In addition, the results for each degree of approximation is given 

in terms of expressions proportional to f(b) and f’(b) where f(b) 

is a polynomial in b (b is a dimensionless parameter defined by b = /cR). Also, 
the final results are given in terms of functions of R (R is defined as the 
distance between the atoms) with an explicit R~^ dependence to show the 
modifications to the electrostatic interaction energy due to the electromagnetic 
field interactions. 

The total interaction energy of the system may now be written as 


AE = e^ I AE , , + AE, + AE 

\ d-d d-q q-q 


(236) 


where each of the above quantities is defined in equations (162) , (201) , and 
(235). In addition, these quantities are expressed as sums of terms defined by 

X(j) which correspond to various types of interactions between the electro- 
static and field operators. In general, the breakdown is as follows: The 

electrostatic interactions are defined by X(0) i = 1,2,3. The terms resulting 

solely from electromagnetic field interactions are represented by X(l) 

These terms contain the interactions due to the second-order electromagnetic 
field term in the Hamiltonian, which is proportional to A(p)«A(p). The next 

group of terms, denoted by X(2) , includes interactions resulting from 

considering both the first and second order terms in the radiation field 
perturbation given by A (p) . P (r) and A (p) *A (p) , If the problem had been 
treated using only the first-order term A(p)*P(r) in the Hamiltonian, then 

X(l) and X(2) would not have entered into the discussion. The terms 
resulting from the interaction between the electrostatic interaction operator 

H^^^ and the field operator A(p).P(r) are denoted by X(3)^^\ Finally, the 

— «• — » — ► — ► (i) 

terms resulting solely from the field operator A(p)*P(r) are given by X(4) 

If the problem had been treated considering only the first-order field operator 

A(p)*P(r) and the electrostatic interaction operator then only X(3)^^^ 

and X(4)^^^ would have contributed to the interaction energy. In this way, the 
various contributions to the interaction energy may be studied separately. 
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In subsequent discussions, the dipole-dipole case as well as the higher 
approximations are again considered. The behavior for large and small 
R (R ^ ?C) for the dipole-dipole case is considered. The dipole-quadrupole 
result shown in equation (201) is expressed in terms of polynomials of the form 

f(b) e , f’(b) e . A further transformation in terms of real quantities 


-y "2y 

f(y) e and f’(y) e is made on the resulting expressions. The results 
are then examined for the limiting case of large R. The small R case is 
immediate since the radiation interaction terms contain a coefficient a®, which 
for small R (a = R/X) causes these terms to vanish leaving only the term 
( 2 ) 

X(0) ' , which corresponds to the electrostatic interaction only. [See 

equation (246). ] Finally the quadrupole-quadrupole results are considered 
again and are expressed in terms of polynomials in b, as for the previous 
cases. 


Dipole-Dipole Approximations 


The first check on this calculation is provided by the comparison of the 
results obtained here with those given by Casimir and Polder, and others. The 
expression for the dipole-dipole approximation generated here and given by 
equation (163) is shown to correspond explicitly to the results reported by 
Power and Zienau [4] . In obtaining this first-order result by means of straight- 
forward stationary state perturbation theory, one also obtains additional 
information regarding the types of interactions involved in the overall result. 

For instance, one can show explicitly how the electromagnetic field interactions 
combine with the electrostatic interaction to absorb the factor resulting from 
the purely electrostatic interaction, thus allowing one to write the interaction 
energy AE^ ^ as in equation (163). One also can show precisely how the 

±i b 

contributions from the integrals over f(b) e are eliminated by the residues 


at b = 0 of the integrals over f (b) e 


±2ib 


which exactly cancels the contribution from X (0) 


Also, these residues give a factor 

( 1 ) 


[See equation (163). ] 


The resulting expression for AE^ ^ is then found to result only from the 
integrals containing f'(b) It is for this reason that the higher order 

results are given in terms of integrals over and 
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The behavior of ^ [defined in equation (164)] for small and large 

R is of interest. The comparison of the magnitude of R discussed here is 
made with respect to the characteristic wavelength X corresponding to the 
frequency associated with the is 2p atomic transition. Using the relations 
(E^ - Eq) =H2i:v, c = Xv, one obtains the connection between X and R, 
given by (E^ - Eo)R/lic = {2ir/X)R. In addition, if one defines K = X/2ir, then 
a = R/X. Thus, in the case in which the atoms are separated by a distance R 
such that R » X, the result given in equation (164) becomes 


AE , ,(R » X) 
d-d 


4 (El - E,)\<Tf><Tjy 
(SpffR^lficp 


X 


oo 


; 


-2uR 


u^d u e 
[(El - Eo)/hc_‘ 


T 


X 


2 5 . 6 , ^ 

^ T^) ^ ^ (uR)3^ (uR)^ 


(237) 


The above result is obtained by writing the denominator of the integrand in 

equation (164) as [( (Ej - Eg) R/hc } ^ + { uR } and then neglecting { uR} . 
This is allowed since (Ej - Eo)R/hc >> 1 and the leading contributions to the 
integral in question result from small values of { } due to the exponential 

factor e ^ . This is the same as neglecting (^ck) or (hc/c') in favor 

of (Ej - Eg) in the denominators of the form [(Ej - Eg) + HckJ found in the 
expression for AE^ ^ given by equation (134). There the major contribution 

to the integrals over k and k* comes from those values of k and k' such that 

K, k’ « (Ej - Eg) /he on account of the exponential factors e ' in the 

integrals [10]. Factoring out R"® from the integrand of equation (237) and 
rearranging factors, one obtains 
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AE^_^(R » X) 


4^c 

(3)2 ttR* (El - Eo)^ 


X J (Rdu) e { (Ru)^ + 2 (Ru)2 + 5(Ru) 2 + 6(uR) + 3 } . 

0 

(238) 

oo 

Integrating term by term, using the integral formula J x e dx = n!/t , 

0 

the integral factor in equation (238) is evaluated, and the resulting expression 
is 


4 he <Cr2^<Cr?2)> / 23 \ 
AE- (3)2 7 tR^ (El - Eo)2 \ 4 / 


^d-d' 


(239) 


Casimir and Polder [1] give this result in terms of the static polarizabilities of 

2e2 ^r^)> 

the atoms defined by o; (I) Vt • Hence, including e^ in the final 

6 (lij - I^q) — 

result for AE^ ^ and rearranging numerical factors, equation (239) can be 

shown to be equal to equation (56) of Casimir and Polder's paper [1] . Equation 
(239) shows explicitly the R~^ behavior of AE^ ^ when the internuclear 

separation distance R is large in comparison with the characterisitic wavelength 
X. Equation (239) may be put in a simpler form by expressing it in terms of 

X(0)^^^ as follows: 


^ /T, 23 he X(0) 

^®d-d^^ “ 3 7T R (Ei - Eo) 


(240) 


Letting y = uR in equation (163) , one obtains 
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AEd_d = X(0) 


( 1 ) 


4 r Rdu { (uR)^ + 2(uR)^ + 5(iiR)^ + 6(uR) + 3 } 


3 7T 


/ 


[a^+ (uR)2] 


AiI><£iX 

3 (El - Eo)R® 


(a) 


(241) 


Comparison of equations (240) and (241) shows that for large R, (a) 

decreases monotonically with increasing R (Fig. 15). 



Figure 15. Behavior of the correction factor § (R/?t) for various 

internuclear separation distance R/?t. 

The behavior of AEd d small R may be obtained by considering 

the factors of equation (164) as follows [10]. The denominator consists of two 
small quantities; thus, one cannot neglect one against the other. On the other 

— 2uR 

hand when R is small, e may be replaced by unity since the leading 

contributions to the integral in equation (164) result for small values of (uR) . 
Finally, the leading term in the polynomial in (uR) corresponds to (uR)“^. 
With these modifications equation (164) becomes 
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oo 


du 


AE^_^(R « X) 



(242) 


Writing the integral as f — ^1L_ 

[((Ei-Eo)/iic}2 + u2] 


, and using the 

2 U j 

a^dx/[a^ + x^] =7r/2a^ , equation (242) reduces to 


0 


AE^_^(R «X) 


3 (El - E«) E‘ 


To express the above result in more familiar units, one needs to include the 
factor as indicated in equation (236). Using the following relations, [see 
equation (86)]: 


/ 2i \ ^ (2S. + 2 ) ! / an \ 
^ 2 \ 2Z / 


2i 


E = _ 

0! 2 an 


and including the e^ factor, AE, ,(R«X) becomes 
— d-d 


AE^^^ = e^ AE^ ,(R « K) 
d-d 




(243) 


Note that when (Ej - Eg) is replaced by an average energy [6], the (-8) 
factor in equation (243) becomes (-6). Equation (243) is just the result 
obtained when only electrostatic interactions are considered; it corresponds to 
the London-van der Waal's interaction energy in the absence of radiation fields. 

Thus, in the limit of small R, §^^\a) goes to unity, showing that for small 
separations the interactions with the radiation field are unimportant and for 
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large separations the introduction of the radiation field gives rise to a 
weakening of the van der Waal’s force. 

The behavior of for various R values (R » a^) is 

illustrated [1] in Figure 15. 


Dipole-Quadrupole Approximations 


The detailed analysis of the higher approximations is not as straight- 
forward as in the dipole-dipole case. Before one can analyze the equations 
corresponding to the dipole -quadrupole order interaction energy ^ given 

by equation (201) , one needs to simplify this expression by substituting the 
definitions for 3 ^ 2 ^ ••• ®^c. Analysis of equation (201) shows that each 
( 2 ) 

X(j) consists of a group of elements having the same functional dependence 
as in the dipole-dipole approximation, but with the powers of b in the various 
denominators decreased from b“^ to b“^. In addition to these terms, each 


( 2 ) 


group X(j) 

of Fa, ?, 5a. 


j = 1, 2, 4, contains a set of "new" quantities defined in terms 


Using the relations given in equation (201) and the following equalities: 


{ 12 y + 36 ?2 - (3) (52) Fg I = 12 ( J - 4 Fj) + 36 (Jg - 3 Fg) 
{ 4 5 y “ + 6 ?2 + (2) (26) Fg fJ } = 4 { jy " + 4 Fg fJ } 

+ 6 { ?2 3^2 + 6 Fg Fg } 

{ y-4Fg} = j ?2-3Fg } 


14. This effect was proposed by Verwey^ Overbeck and Nes in their book. 

Theory of the Stability of Lyophobic Colloids (Elsevier, Amsterdam, 1948) 
and led Casimir and Polder to consider the problem. 
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the expression for AE^ ^ defined by equation (201) becomes 




( 2 ) 




* f *bMa+^b)2 { 4 F 3 } db + 2i / ^ { F + 2 G G'*' } db 


- 8 i / 
0 


b=i y?-.4F.Ft|db j” 


b^ (a + b) 


where the above functions are defined as follows^®: 


-ib 


( F + 2 G } = (b^ + Sib - 3) - (b^ - Sib - 3) e 

I?- 4 F 3 J = ^ ^^ 5 - j^(ib^-6b^-15ib+15) + (ib^+6b^-15ib-15) e 


-i b 


{f f' 


+ 2 G g"^ 




2 i b 


(b^ + 2ib^ - 5 b^ - 6 ib + 1) - (b^ - 2 ib^ - 5 b^ + 6 i b + 1) 


-2 ib 


{yr + 4F3F^ = 


(-b® - 4 ib^ + 14 b^ + 42 ib® - 81 b^ - 90 ib + 45) e 


2ib 


- (-b® + 4 ib® + 14 b^ - 42 ib® - 81 b® + 90 i b + 45) e 


-2ib 


{ ?2 ?2 + 6 F3 F3 } = (-b® - 6 ib® + 27 b^ + 84 ib® - 162 b® - 180 ib + 90) 


(-b® + 6 ib® + 27 b^ - 84 ib® - 162 b® + 180 ib + 90) e 


-2ibj 


15. These relations are obtained by using the various definitions for 
F, G, ... 3 ^ 2 > forming the appropriate sums and products. 


(244) 
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Comparison of equation (244) to the dipole-dipole results given by equation (162) 

ii b 

shows explicitly the terms of equation (244) having the same f (b) e , 
f'(b) dependence as in the dipole-dipole case. 


Using the techniques outlined in Appendix C, the b-integrals appearing 
in equation (244) are individually evaluated. Incorporating the results of these 
transformations equation (244) becomes^® 


AE, = X(0) 
d-q ' 


( 2 ) 




- ( — ) + 24 / ~ ---T (y^ + 6y2 + 15y + 15) e"^ dy 

0 y^a2 + y2)' 

^ /|i) _ j (4y^-2a|j (y4^2y’+5y^+6y+3) dy 

0 (a2 + y2) 

ao 

-2('^j-4 f (y® + 4y® + 14y^ + 42y® + 81y2 + 90y + 45) 

oo 2 \ 

+ s( — ) + 6 f — y (y® + 6y® +27 y^ + 84y® + 162y® + 180y + 90)e"^^ I 

^ ^ ^ 0 y®(a® + y®) j 


(245) 


The constant factors in the above equation are the residues at the origin (b = 0) 
of the various terms considered. Adding these constants, one sees that their 
sum is zero. In addition, they combine in pairs; the residues of the first and 

2 -y 4 -2y 

third integrals proportional to (y . . . ) e and (y . . . ) e add up to zero, 

3 -y 6 “2y 

and the "new" terms proportional to (y . . . ) e and (y . . . ) e also add 
to zero. Considering these terms in two groups, one sees that the group 
associated with the "new" terms has an additional y~^ factor, which gives rise 


16. The ordering of terms in equations (244) and (245) is the same to show 
explicitly the results of the various integrations and transformations. 
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I 


■ IMI Hill 


I I 




II 


II IIIIBI I 


to a lower dependence. Rearranging^^ terms in equation (245) to reflect 


the preceding groupings and excluding the constants, AE^ ^ may be written 
as follows: 


"d-q 


AE^q = X(0) 


(a f (y=.3,.3).->-dy 

(3) (57T) ^ ^ 


^ r (4y - 3a^) ^y4 2y3 + 5y2 + 6y + 3) e' 

0 (a2 + y^) 


2y 


/ 24 / (/ + 6y^ + ISy + 15) 


(3)''(5tt) 


0 (a^ + y^} 


2 r — 2 _ (y6 + ioy® + 53/ + 168/ + 324/ + 360y + 180) e 

0 y^(a^ + /) 


oo , 2 

+ 4 r — — 5 (y® + 4y® + 14y* + 42y® + 81y^ + 90y + 45) e 

• 0 y^(a^ + y^) 


X(0)^^^ + x|^^a) +X^^^a) 


(246) 


where the new functions 



and X 


( 2 ) 

2 


(a) 


are defined in the expression 


for AE . . 

d-q 


Subsequent calculations will show the significance of the above 


definitions. Comparison of equation (246) to its dipole-dipole counterpart 
given by equation (164) shows that, unlike the dipole-dipole case, equation (246) 
consists of three separate terms, with one of them being solely due to the 
electrostatic interaction between atoms and the others arising from the inter- 
actions between the electrostatic and electromagnetic fields. 


In order to further analyze the results for the dipole -quadrupole case, 
one needs to examine the large R (R » K) behavior of equation (246) . Follow- 
ing the same procedure as that used for equations (237) through (243) and 


17. The last two terms in equation (245) are combined by first getting 
2 

y2 (a^ + y2) as the conomon denominator. 
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collecting all the terms having the same power and similar exponential 
factor, the preceding expression becomes 

AE, in»K) 
d-q 


= X(0) 


( 2 ) 


>2 ^3 


1 + 


2" a’ 
(3)"(57r) 


-- / 
aj J 


(y^ + lly + 51 + 120y-‘ + 120y"^) e"^ dy 



(2y^ + 12y* + 58y^ + 174y + 327 + 360y~* + 180y~^) e 


-2y 


dy 


f (y« + lOy® + 14y^ + 42y^ + 81 y^ + 90y + 45) dy 


(247) 


Evaluation of the above integrals is accomplished in the following way: For the 
integrals having positive powers in y, one simply uses the formula 

J X e dx = (n!/ 1 ) . The integrals containing negative powers of y may 

0 

—y 

be combined by letting y 2y in the integral containing e , resulting in the 
relations 


Jy^e^dy = Jy^e^^dy , 


0 


0 


OO OO 

r -2 -y . 1 r -2 -2y 

Jy e dy = — Jy e-'dy 


Doing this, one finds that the sum of these integrals is identically zero. Per- 
forming the indicated operations and combining the resulting factors into two 
groups having a and (1/a) coefficients, equation (247) becomes 
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(3)3(57t) 




( 2 ) ( 2 ) , 7 
Using the definitions for a and X(0) , one can see that aX(0) ~ 1/R 

( 2 ) 

and X(0) ' /a ~ 1/R®; hence, the above expression may be put into the 
following form 




where 


ISSirhc 




(El - Eo) 


1854 he ^ri^Xr^^ 
135 7T (Ej - Eo)^ 


Equation (249) shows explicitly how the electrostatic interaction term B/R® 
is reinforced by a factor proportional to 1/R® and decreased by a factor 
proportional to 1/R^. To determine the relative magnitude of these factors, 
one can rewrite the above coefficients in terms of atomic units. 


Referring to the equation for AE , given in terms of factors 

Q—Q 

proportional to R~^, R“®, and R“®, it is of interest to note that this result 

should be expected from the general form of AE , [see equation (246) ]. 

Q—q 

Hence, one may write the interaction energy as follows: 
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(250) 


AE 


d-q 


(RA) 


9 ^( 0 ) 


( 2 ) 


(R/K) 


where the above functions are obtained by simply factoring out a~^, R~®, 
and a”^ from X (a ) , X(0) and X^^\a), respectively. [See 
equation (246) .] 


The behavior of the above quantities in equation (250) may be inferred 
from Figure 15 and by comparison of ' (R/?C) and g' ' (R/X) to the 
function g^^^(R/X) of equation (241) . Since the functions appearing 


in 


AE, are much more complex than g^^^R/K), a computer program would 
Q^q 

be required to numerically evaluate the integrals to construct the plots 
corresponding to Figure 15. 


Quadrupole-Qua{Jrupole Approximations 

The increased complexity of the expression for the quadrupole- 
quadrupole order interaction energy AE^ ^ given by equation (235) makes 

it necessary to give only the general expression equivalent to equations (162) 
and (244) for the dipole-dipole and dipole-quadrupole cases, respectively. 


Starting with the expression for AE given by equation (235) and 

q— q 


recombining the various terms as in the dipole-quadrupole case, AE 
may be written as 


q-q 
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^Vq " ^ (3)(5>(7f) j J ^3db [ 3a(a+b)=> J 


+ /b(J-4F3) db r. 8(6b3.ipab^-yb . ^4a3) 

J ' L a(a+b)2 


r Jf + 6 ji - 3 Fa) 

0 

r 8 (33b® + 40ab® - 40 a®b + 64 a®)! 

^ [ MiTTbP J 


7 r h2 

+ / (19b® F2 - 2 ?2 + 6 F3) db [^(77^2 


+ 1 


+ 1 


. / [r .V 3 0 o*] 


7b.db[<r-oHF^-G^,][-i2fy^; 


+ i / b®db [jy + 4 F3 F3] 1^- -7] 


+ i 7 b^ db [?3 ? 3 - + 6 F 3 F 3 ] 



32 a b* ' 

[gg J 

_ (a + b)_ 


00 

+ i J b®db Ff r + 12 ?4 + 24 Fj F 3 1 

o'- 

r 8 (5b® - 5ab+ 8 a®) 

^ L 3 (a+ b)® 


( 251 ) 
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The preceding expression is no simpler than equation (235) , but it is more 
useful since the various functions have been regrouped into sets previously 
considered, which simplifies further analysis. The "new" groupings generated 
may be expressed in terms of polynomials in ^ using previous techniques. 

The new combinations required are as follows: 

|f+6y4-3F3| |(b^+ 4ib® - 39b2 - 105ib - 105] 

- [b^ - 4ib^ - 39b^ + 105ib - 105] e 
{ f + 12 ^4 3^4^ + 24 F 3 F 3 } 



|[[b*+ 8 ib'^ - 40b® - 144ib®+ 384b^+ 768ib® - 1152b2 - 1152 ib+ 576] 
+ 12 [-b® - lOib® + 49b^ + 144ib® - 264 b^ - 288ib + 144] 


+ 24 [b* + 6 ib® - 15b^ - 18ib + 9] J 

-[ [b® - 8 ib^ - 40b® + 144ib® + 384b^ - 768ib® - 1152b2 + 1152ib+ 576] 


+ 12 l-b® + lOib® + 49b^ - 144 ib® - 264b® + 288ib+ 144] 


+ 24 [b^ - 61b® - 15b® + 18 ib + 9]] e ^ ^ 

Using the above relations and the dipole-quadrupole approximation, 
equation (251) may be further analyzed. It is not necessary to rewrite 
equation (251) explicitly in terms of polynomials in b, since only a comparison 
of this result to previous answers is given. Reference to equations (244) and 
(251) shows that the quadrupole-quadrupole result contains groups of terms 
appearing in the previous dipole-quadrupole case but with different polynomials 
in a and b. Comparison to equation (162) also shows that the quadrupole- 
quadrupole case again contains terms having the same functional dependence 
ii b 

f(b) e as in the dipole-dipole case. Thus, one can see that the results are 
systematic, in that the higher approximations can always be expressed in terms 
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of previous quantities appropriately modified, plus some "new" terms. An 
analysis can be made of the dimensionality of the various integral terms of 
equation (251) in terms of the constant a. Since a is proportional to R, 
one can see that terms proportional to R~^, R”®, . . . R~^^ are included in the 
various elements of equation (251). For example, the eighth integral of that 
equation is an example of terms having an R"'^ dependence since its value 

/ o \ 

is independent of a and its coefficient is a®X(0)^ . The rest of the integrals 

have various factors proportional to 1/a, 1/a^, 1/a^ and 1/a^, which when 

/ o\ 

multiplied by a^X(O)^ gives rise to the R~^ dependence indicated above. 
Hence, one can always express the results of equation (251) as follows: 


AE 

-q 



^ R» 





(252) 


The algebraic form of each of the above functions can be obtained by separating 
the terms in equation (251), as was done for the dipole -quadrupole approxima- 
tion. When this is done, one finds that each function is the same as in previous 

ii b 

cases with only more complex factors f(b) e . Therefore, their behavior 
for all R is expected to be like that on Figure 15, with only a more rapid 
decrease for those terms whose limiting form for large R decreases as R~® 
and R“^^ . 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, April 10, 1970 
948-70-70-27-00 
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APPENDIX A 


POLARIZATION VECTOR SUMS 


In the course of this calculation, various products involving the 

polarization vectors e (k) need to be evaluated. To calculate these relations, 

A. 

one defines a coordinate system (Fig. A-1). 



Figure A-1. Polarization vector coordinates. 


Using the coordinate system in Figure A-1 and the following unit 
vector relations; 

€i = ■^n3Cj + ij2^2'*'^13^3 > 

^2 = -^21 -^22 ^2 -^23 ^3 » 

fC = -^31 + i32 3^2 ^33 ^3 » 

e^* €3 ~ K — 0, ^2"^ ~ 


A A A A A A 

6i.€i = €2*£2 = K-K = 1 

one obtains the various sums needed. 


(A-1) 
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The first sum to be considered is evaluated as follows: 


Z 

X=1 


A *r A — * A A A A 

e^(K)*€^(K) = ei*ei + €2«€2 = 2 


(A-2) 


The next sum is given by 


2 2 






Z Z ( 

\=i x'=i ' 

( A .A, A A,, 1 2 ( A ^A A A ) 2 

= I ei(K)-€i(K') I + I €j(K)-e2(K')| 

A, A A, ) 2 I A A^ A A. ) ; 

(K)'ei(K')| + { €2 (k) • €2(k )| 


A 

+ \ €2 


In terms of the unit vectors of equation (A-1) , the above becomes 


y ( i.. il. i. . . + i.. J?'. . 

^ ( li li Ij Ij li 2i Ij 2} 


+ 

2i li 2j li 2i 2i 2j 2j 


Adding and subtracting and rearranging factors, this 

ij ^ ^ J 

expression becomes 


y i . - T .a 

H 11 Ij Ki K] ¥. li Ij 3i 3j 
13K IJ 


y 2i 2j K1 K] Af 2i 2j 3i 3j 
IJK 1] 
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Combining the sums over i, j and using the sum rule for orthonormal systems 

Z SL . i . = 6... the above quantities combine into 

/Cl /cj ij 

If 



2i 2] 



- -^o- ^o- 

3i 3j 


Adding and subtracting 


E 

ij 


3i 33 


5.. - one obtai ns 

ij 3i 33 1 ’ 



Kj 






Summing over /c in the first sum, one obtains 



which when combined with the second term, gives 



A A A A 


A ,A A 





A A 
(K).(/C). 



{k').(k'). j , 


(A-3) 


where one uses the relation /c = (/c)ifti + (k)2^2 (<^) 3^3 one in 

equation (A-1). Equation (A-3) may be expressed in terms of the angle 
between k and /c* defined as 0 . To find this relationship, one takes the 
result in equation (A-3) and expands it over i and 3 . Hence, 
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^ \ ij 3l 3j/ \ 13 3i 3j^ 

= 3-(i ,+£ „i +£ i )-(i' i’ +i' i' +i’ je’ ) 

\ 31 31 32 32 33 33/ V 31 31 32 32 33 33/ 

+ (1 L L' i' +i I i' i' +i i i' i' ) 

V 31 31 31 31 32 32 32 32 33 33 33 33/ 

+ 2fi i' i’ +i i i’ /' +i i /' i’ \ 

\ 31 32 31 32 31 33 31 33 32 33 32 33 / 


Combining the above terms, one obtains 


I ^ 2 (-^31 -^Si) ^ ("^Sl ■^3l)(^32 ■^ 32 ) (■^31 ■^ 31 ) ('^33 "^ 33 ) 

^ ("^32 '^32) (■^33 "^33) 

Analysis of the last three terms in the above expression shows that these 

elements correspond to the cross terms of the product ^ (^^3i ^3i) if 3] 

/ \ 2 ^ j 

and that the term 2 j ( "^Si ) corresponds to the sum of the squares of this 


product. Hence, 


E E I 

\=i A'=i y 






(A-4) 
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The next sum consists of polarization vector components given by [10] 


' ('i)i('i)/ ('2). ('2). 


= y i . i . 

^ K1 K} 3l 3j 


= jSy- ) . 


(A-5) 


The next term to be evaluated consists of products of the form con- 
sidered above and given by 


E E 

A=1 \'=1 






(A-6) 


The last sum to be evaluated consists of combinations of terms 
previously considered. It is evaluated by taking 


2 2 


^ \ Ik Ik li 1] Ik 2k li 2j 


+ ■«„. 

2k Ik 2i 1] 2k 2k 2i 2j / 


= js.j- (^')^(k').- (x).(x)j+ (x). (/<•). (x*k’) I 


(A-7) 
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APPENDIX B 


PROPAGATIONAND POLARIZATION VECTOR 
ANGULAR INTEGRALS 

Introduction 

In this appendix, the integrals over the angular coordinates of k and 
k' are evaluated in terms of functions depending on the magnitude of k and k\ 
In the course of evaluating these integrals, one uses the following relations: 


b = kR , 

/3 = k'R , 

A A A A 

K = sin 6 cos <p sin 6 sin cp 62 + cos 6 03 
X = ±i/cRcos0 = ±ibcos0 , 


cos 9 = cos 6 cos 9' + sin 6 sin 0* cos - <p') , 


/ 


„m ax , 

X e dx 


m ax 

X e m 
a ~ a 


r m -1 ax , 
j X e dx 


(B-1) 


The single integral functions over dJ2 are evaluated first; then using these 

K 

results, the double integrals over dfi and dS2 , are considered. Finally, 

K K 

the last section contains the relations satisfied by various terms having 
exponential factors ±i (k+k*) • R and ±i ( k-k*) • R. 
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Evaluation of 



cos*^ ed^K 


/ ±i /c* R 

e dJ2^ . 

After substituting for the quantities defined in equation (B-1) , this integral 
becomes 


7T 27T 

J dJ2 ~ f f ® 

0 0 


+ ib 


/ 2 7T \ r ^ / sin b \ , ^ „ v 

■= (vTb) / ® 


±ib 


J- 4 ~l 

e cos 0 dS2 . 

K 

Using previous definitions, this integral becomes 


/ 


-ik-R „ _ 
e cos 0 dn = 

K 


„ + i b 

2tt r X , 


-ib 


, . ( cosb sinb ^ _ , . , ^ , 
= 4 7T1 ]g— - 1^2 ) = 47rib G(b) 


Similar operations are used to obtain the corresponding integral with a (+) 
sign; the results are 


/ 


±1k-R 

e 


cos 0 dS2 

K 


+ 4irib G(b) 


(B-3) 


The next integral is given by f ^ cos^ 0 dS2 . Using the same 

substitutions as before, this integral gives 
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r. 

/• ±i/c*R ? „ 2 7 t r 9 ±x dx 

e cos^ 0 dn = ■ ■ ■ ■ . ,2 I ^ e 

±ib 


(+ib) 


Choosing the upper sign in the above expression, one obtains 


r +i K* R 2 /I 
J e cos^ 6 dS2 


^ f sin b 
= 47T j — — 
K I b 


sinb 2cosb 
+ — ? — 


2 sinb 
b3 


Similarly the other integral gives 


f -i/C'R 2 f +i«'R 2 

J e cos^ 6 dn^ = j e 


cos 6 dS2 


Combining these results, one obtains 


J ^ cos^ 0 dn = 47T (F (b) + G(b) } = 4 7t F2(b) , 

, ^ N _ sin b . ^ ^ . cos b sin b „ 

where F(b) = — g— + G(b), and G(b) = — . Further 

combination yields 


J (1-3 cos^ 0) = - (2^ {F(b) + 2 G(b)}= - (237 t) 

The next integral to be considered is given by J ^ cos^ 

Using the same techniques, one obtains 


(B-4) 


F3(b) . 

0 dn . 

K 
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J 


cos^ e dS2.. = + 47Ti j bG(b) Fg(b) = + 47rib-^ ?2(b) . 


K 


(B-5) 


/ii K • R» A 

e cos* 0 dS2^ . One can show 

that this integral is equal to 


J cos^ 0 dn^ = 47 t jFg(b) + G(b) Fg(b) j s47rb-2f(b) , 

(B-6) 

where Fg(b) = F(b) + 2 G(b) and f(b) = b^ (Fg/b) + G(b) - 8 Fg(b). 


Evaluation of JJ cos*^ 0dJ2^df2j^i 

The first integral of this type to be treated corresponds to n = 0. 
Separating the k and k’ variables, one obtains 


r r ±i(k+k')-R r r 

e df2 dS2 , = e dG e dfl , 

J J K k' J K -J k' 

When comparing this expression to the results of equation (B-2) , the above 
integral gives 

/ sinb \ / sin /3 \ 

V b A /3 y 

(F(b) - G(b)} |F(/3) - G(;8)} . 

(B-7) 


If 


±i(K+K') • R 


dS2 dJ2 , = 2^tt^ 

K k’ 


= 2*7t2 
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The next integral is given by J J e cos 0 dSl^ dQ^, . Since 

cos © is a function of both primed and unprimed variables, one needs to con- 
sider this expression as follows: Taking 



^+i ( k+k*) • R 


cos © dfi 

K 


dn , 

k’ 


= / 


-i/c R cos 6 
e dn 


J [cos 6 cos 0* 


• ^ -iK*Rcos0* ^ 

+ sm 6 sin 0* oos(<p - ^•)J e dS2 , 

K 

27T 

shows that the integral over d<^' is simply J cos (<p - <p') d<^>'. Using the 

0 

identity cos (0 - 0*) = cos 4> cos 0' + sin <p sin <p\ one readily sees that the 
integral over (f> vanishes when integrated from zero to 2 tt. Thus, the above 
integral may be evaluated using equation (B-3) , as follows: 


r r ^ cos © do dfi , 

J J K k' 

r ±iKRcos 0 „ r ±iK'Rcos 0 ' 

= J e cos 0 dn^ J e 


= - 2^ {bG(b)} {/3 G(|3)} 


(B-8) 


The next integral is given by JJ ^ ^ cos^ © dJ2^ * 

Expanding cos^ ©, one obtains three terms, one of which contains cos (<^ - 0'). 
Eliminating this term, this integral becomes 
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r r ±i(/C+/C*)*R 2 ^ 

I I e cos^ O do. dS2 , 

J J K K' 

= f f ^ ^ cos^ e cos^ 0 ' dS2 dJ2 

J J K k’ 

+ JJ ) R g gjjj2 g! gQg2 (j£2^ 


Integrating over d(p, as before and factoring the resulting expressions, the 
above equation becomes 


ff cos^ 0 dS2 df2 , 

= f ^ cos^ 0 dn r ^ cos^ 0* dJ2 , 

J K J k’ 


+ 2 71^ J e^^ ^ ^ sin^ 0 sin 0 d0 e"^^ ^ ^ sin^ 0 • sin 0 d0 


Using the identity sin^ 0 = (1 - cos^ 0) , the second term in the above equation 

is transformed into integrals over cos^ 0. Using previous results in equations 
(B-2) and (B-4) , the above becomes 


f f cos^ 0 dn dS2 , 

J J K k' 

= 2^ TT^ { [F(b) + G(b)] [F(/3) + G(/3)] + 2 G(b) G(/3) | . (B-9) 
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Combining equations (B-7) and (B-9), one obtains [2] 


r f ^ ^ (1 + cos^ ©) dS2 dfi , 

^ K k' 

= 2® 7T^ { F(b) F08) + 2 G(b) } 


The next integral is given by JJ e"*^^ ^ ^ cos® © dO^ dfi , . 

Expanding cos® ©, one obtains four terms, one of which is independent of 

cos - (p^) and the other three are proportional to cos*^ (<p - <p') > = 1> 2, 3, 

respectively. Hence, one needs to evaluate J cos® (cf> - <^’) d(p d(p' by 

expanding into products of sin (p and cos (p. After doing this and integrating 
term by term, the answer is zero. Incorporating these simplifications, the 
term under consideration reduces to 


JJ («■*■« ) Q 


K K' 


- // 


±ixRcos0 ±i/c'Rcos0‘ . 


sin d de sin 6' d0* 


X { 10 cos® 0 cos® 0' - 6 cos 9 cos® 6' 


- 6 cos® 9 cos 0* + 6 cos 0 cos 9' } 


Using the results obtained in equations (B-3) and (B-5) , the above integral 
reduces to 
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rr dn , 

J J /c /c' 

= (-4,^, j.o[bG(b,-^][^G(«-^] 

- 6 [bG(b)] 1^/3 G(^) - ^ y- ] 

- 6 [bG(b) [)3G(^)] 

+ 6 [bG(b)] [/3G(i3)] I 
„4 2 {3^2 (b) y,(0) + 6 F.,(b) F^m) 

— - ^ TT /!_ j 


where ? 2 (b) = b^ G(b) - 2 FjCb) . 

Combining equations (B- 8 ) and (B- 10 ) , one obtains the following 
relation: 


J J (cos © + cos® ©) ^ ^ 


2^7r2 


[bG(b)][) 3 G(^)] I 


+ j bG(b) ^G(/ 3 ) - 2 bG(b) 


25 ^ [y(b) y(g) + 4 F.,(b) Fa(g)' 

b^ 
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r r -i-i ^ ) * R a 

The next term is given by j J e cos* 0 dS2^ * 

Expanding cos'^ 0, excluding the terms containing cos {(p - p') and 
cos^ ip - p') , and noting that JJ cos'* (p - p^) dp dp' = 3/2 , the term 

under consideration factors into the following form; 

ff dn , 

^ ^ K fC' 

= ^ [s j 0 cos^ 9 sin e de f e*“’^ cos‘ 6 ’ sin fl- de' 

+ 24 / e±iKR cos e g e de f °°® ®' (cos== fl' - cos'* 0 ’) sin d' de' 

+ 3 J ® (1 - 2 cos' 0 + cos^ 0) sin e de / ®' (1 - 2 cos' 0' + cos' 0') sin 9 ' d0>] . 

Using previous results given in equations (B-2) , (B-4) , and (B-6) , the above 
integral is given by 


ff 


±i(K+/C*)*R 4 


cos^ 0 d^2 (K2 , 

K k' 


_ 

2 


(8) (4) {F3(b) +G(b) Fg (b)j { F3(/3) +G(/3) Fgt^)) 


+ 


+ 


(24) (16) I G(b) - I I G(^) 

(3) (16) (16) {^}{^} 


4 FM 

13^ 


= 2* 7t 2 [f(b) f(^) + 12 J4(b) ?4(^) + 24 Fg(b) Fg(^)] 


(B-11) 
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where, f(b) = F 4 (b) - 8 Fg(b) and ? 4 (b) = b^ G(b) - 4 FgCb). Combining 

the results of equations (B-9) and (B-11), one obtains the relation 


ff ^ ^ (cos^ 0 + cos^ 0) dS2 dn , 

JJ K K* 

= 2^ 7t 2 |[F2(b) Fa(^) + 2 G(b) G(^)] 


^ f(b) f(/3) + 12 5 a ( b) + 24 F«(b) F;,(^) | 

b2 /S2 / • 

Evaluation of ^ cos^ 

The first term that needs to be evaluated is given by 
f* r ii ( /c “fc ^ 

J J e dn^ the previous techniques and definitions, 

this may be written as 


II 


±i(/c-K’) • R 


dS2 


r ±1 • B. r 

dS2 , = I e dfi I 

K k' K-^ 


+iK».R _ 
e dS2 


- ff 


^±i(k+k')-R dS2 , . (B^12) 

K K’ 


The next term may be written as 


167 


rr dbi(fC— K*)’R , j 

I I e cos 0 dJ2 dfi , 

'> K k' 

r ±iK-R ^ r +!/{•• R 

= I e cos 0 dn I e cos 0' cK2 , 

J K k' 

= +2^7r2{bG(b)}{|3G(i3)} 

rr ibi(K+K^)*R _ 

= - I e cos 0 cK2 dn , 

The next term contains cos^ 0 and is evaluated as follows: 

r r k: ) R 0 (jQ dS2 

K 

r ±i/c*R 2 r +1 k**R 2 /.f 

= I e cos^ 0 dn I e cos*^ 0’ dn , 

_ 2 r k*R . 0 „ . „ c +i k**R . 2/%t - «t 

+ 2 JT J e sin*^ 0 sin 0 d0 J e sin*^ 0 ’ sin 0 ' d0 ' 

Expressing sin^ 0 in terms of (1 - cos^ 0) and using the results of 
equations (B-2) and (B-4), one finds that 

r fe*“^"^'’®oos'edS! « , 

^ ^ K 

= f f ^ ^ cos^ 0 df2 dJ2 , 

^ K K* 

The next expression, when expanded using previous results, gives 


(B-13) 


(B-14) 
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J J ^ cos® © dS2^ dn^, 


K K' 


= // 


±i/cRcos0 +iK’Rcos0' . 


sin 6 dd sin 0 ’ dfi ’ 


X { 10 cos® 0 cos® 0* - 6 cos 0 cos® 0' - 6 cos® 0 cos 0’ + 6 cos 0 cos 0* } 


Using the results of equations (B-3) and (B-5) , one sees that the respective 
integral products are just the negative of the corresponding terms in equation 
(B-10) hence, 


// 


±i(K-K')-R 3 

e cos 


0 dS2 dS2 , = - f r 


±i(/C+K*)'R 3 , J 

e cos^ 0 dS2 dS2 


(B-15) 

The next term to be considered involves cos^ 0. Referring to the 
equations leading to equation (B-11), one notes that the terms contain cos 0, 
where n = 0, 2, 4. Analysis of equations (B-2), (B-4) , and (B-6) shows 
that all these terms are invariant under sign changes in the exponential factor; 
hence, 


// 


±i(/C-k')"R 4 ^ 

e cos* © dS2 dS2 , 

K k’ 


CC ±i(/<+K’)'R 4 „ 

I e cos* © do dS2 , 

•) J K k' 


(B-16) 
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APPENDIX C 


PROPAGATION VECTOR INTEGRALS 
Introduction 

In this appendix, the integrals over ^ and b (kR = b, /c*R = ^8) 
necessary in this calculation are evaluated using various techniques employed 
to obtain the principal value of integral functions. Very few specific 
references are given since most of the material can be found in standard 
books on advanced calculus or complex variables. In evaluating integrals of 
the type considered here, one assumes convergence at infinity in all cases 
where the integral considered is evaluated over contours which include circular 
paths at infinity. This requirement is necessary to obtain finite results for 
this type of problem [1] . This aspect of the calculation could have been 
treated by including in each integral over k and k' a convergence factor 
— \y \k 

e , which would have guaranteed finite results during the various 
integrations. After the integrations, this factor is removed by letting |y| = 0 
in the final results. This divergence results when the variation of the electro- 
magnetic field within the atoms is not treated exactly; in fact, in the dipole- 

dipole approximation ( e^*^ ^ 1 1 , it is entirely neglected^®. Casimir and 

I "V I ^ 

Polder remwe this divergence by explicitly introducing the factor e ^ 
in their k and /c' integrals. In this calculation, the addition of this factor 
would only make things more cumbersome and so it is left out throughout 
the discussion. Power and Zienau [4] only mention this convergence factor, 
as is done here, to use the resulting integrals. [See equations (10) and (26) 
of Reference 4. ] 


18. The electromagnetic vector potential is then assumed constant over 
each atom. 
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Integrals 


The first integral to be considered is given by 


f F(^) 

J (P + b) (^ - b) • 


Using the definition for F(j3) in terms of exponential functions, this integral 
becomes 


7 /3^d/3F(/3) 

OO « 

r /3^d/3 

l(P^ + iP- i)e^^ - (P^ -ip - \ 

J (/3 + b)(/3 - b) 

1 (p + b) (/3 - b) 

\ 2 i/s’ 1 


r - i)e^ r d/3(/3^ - ig - l)e~'^ 

j 2i^ (^ + b) (^ - b) ~ J 2i/3 (^ + b) - b) 


Letting 13 -13 in the second integral, one notes that the integrand changes 

into the form of the first integral and the limits change from (+«>) to (-°°) . 
Interchanging limits, the above two terms may be added to give 


dj3()3^ + i/3 - l)e^^ 


/ 


_ 2i^ (l3+h){l3 - b) 
takes the following complex integral: 


To obtain the principal value^® of this integral, one 


^F(z) 


0 1 Z 

dz (z + iz - l)e 
2iz (z + b) (z - b) 



evaluated around the indicated contour. Since no poles are enclosed, the above 
complex integral vanishes; hence. 


19. The explicit notation P. V, J will be used only when there is a question 
as to the value of the integral being considered. 
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Lim 

e-K) 

► oo 


-b-€ 0 -e 

J F(z=/3)dz + J F(z = -b+ee^ )dz + J F(z=/3)dz 
-L 7T -b+e 


0 b-e 0 

+ J F (z = ee^ )dz+ J F(z = /3)dz+J' F (z = b + ee^^) dz 


7T 


L ^ 

+ J F(z=)3) + J F(z = Le )dz 
b+e 0 


= 0 


Rearranging terms and assuming convergence at L — one obtains the 
principal value (P. V. ) of the desired integral by transposing terms. Hence, 


+CO 


P.V. f dzF(z=;8) S 


j -b-e -e 

Lim I J F(z)dz+ J F(z)dz 
e— 0 [-L -b+e 

Lr—oo 

b-e L 

+ J F(z)dz+ J F(z)dz 
e b+e 


= Limj J F(z = -b + ee^®)dz + f F(z = €e^^)dz 
e-K) 0 0 


r iO 

+ j F(z = b + ee )dz 

0 


Substituting for z in each of the above integrals and interchanging limiting 
and integral operations, the above becomes 


+C50 



f ^ 

) TTib 

j i 


[ 2b^ 

r 2 

1 2ib^ 


(b^ + i b - 1) e 


+ (b^ - ib - 1) 


ib 


-ib 
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Hence, 


r e^d/3 F(/3) 

J (|3 + b) (|3 - b) 


7T ITT 

2b^ 2 


b F’^(b) 


(C-1) 


The next integral that needs to be considered is similar to the one 
just evaluated, with only F(/3) replaced by G(/3). Following a similar 
procedure, this integral is given by 


7 /3^d0 GW) 

J (13 + b) (/3 - b) 





(C-2) 


where 


G"'(b) 




(ib - l)e 
2i 


ib 


(-ib - l)e 
2i 


-ib 


The next integral expression is J /3d/3 F(|3) ^ ^ ^ ^ , 

Rearranging the factors and fractions, the above quantity can be written as 


2 TT": — • This integral is the same as the one considered 

J (i3 + b) (|3 - b) 

previously; hence, the result is obtained directly, using equation (C-1). 

The integral under consideration may also be evaluated by taking 

J |Sd/3F(/3) + — 2 ’ ) and separating F {jS} as F(yS) = , 

where I(/3) = (j3^ sin j3 + j3 cos - sin /3). Hence, this integral can be written 


as J d/3 I(/3) -p (- 


1 1 
+ 


h + ^ /3 - b 


Expanding the fraction, this integral 


becomes J d^ 1(13) j" ^ 1^ ( b + /3 ~ - ~b ) | ’ The last step may 

be checked by simply collecting terms under a common denominator on both 
fractions. Hence, one needs to evaluate three integrals given by 
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Since 


oo oo oo 

2 f l(/3)dp 1 r I(j3)d)3 1 r l(p)d^ 

'V J p J (b + /?) b^ J (/3 - b) • 

I + ij 8 - l)e^^ {p^ - ip - i)e \ r- j. ■ . i u 

I(/3) =1 2 ^ " 2 I I ’ integral may be 

separated into two parts as follows; 


CO 5 i/3 


(p^ - ip - l)e 
2ip 


dp 


Letting p -p in the second term and interchanging the resulting limits of 


+00 




integration, the two parts combine into J ^ — d/3 . Performing 

the same operations on J and using the fact that I(/3) = -!(-/?) , 

one can combine the remaining two integrals into J • Thus, 


the integral under consideration becomes 


/ pdP F(P) 


0 




_oo 


(p^ + ip - l)e'^ 
2ip 


dp 


+ 



Ki3)d^ 
(b + |S) 


, F(;S) 



Substituting for i/p and 1/ (b + / 3 ) their respective principal value 120,21] 
defined as 


P 


HP) + iiT 6 (P) 


F 

(p + h) 


ilP- (-b)] + i7r6 [|3 - (-b)] 
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where 





6(/3) = Lim 
€—0 


e 

+ €^) 


5 


then one can evaluate the integrals by using the rules [20] for integrating 
over the ^(/3) and 6(/3) functions. Taking the first integral, one obtains 


/ 


()3^ + ij3 - l)e^^ 
2i/3 


d/3 


CO 


= / dp 


-ns - l)e^^ 
2i 


^W) 


+ in 


f dp 


+ ip - l)e 
2i 


i/3 


6 ( 13 ) 



= Lim (6 
e-*0 


(z^ + i z - l)e 
2 i (z + i e) 


1 z 


+ (-^) = - (f) 


where the first integral doesn't contribute since no poles are enclosed by the 
contour^®. The next integral is evaluated as follows: 


20. Note that in this case one must use the requirement of convergence at 
infinity obtained by the implied convergence factor e 
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■ mill nil 



- 7 4 E/3 - (-b)] 

-OQ. 


CO 

+ i7T f I(/3)d^ 6 [/3 - (-b)] 

-OO 


= Lim ^ 
e-*0 


9 1 z 

(z*^ + iz - l)e 

21 ()3 + b + ie) 



- Lim ^ 
e-*0 


9 -i z 

(z^ - iz - l)e 

2 i (j3 + b + ie) 


OO 

dz + iTT J I(/3)d^ 6^)3- (-b)] 

_oo 


Since no poles are included in the contour corresponding to the first integral 
above, its value is zero. In the second integral, the contour must be taken 
over the negative half plane; hence, a pole is included within the contour. 
Thus, the value of this integral is obtained using the Residue Theorem. 
Integrating over the 6-function, the above expression becomes 


/ 


(b + /3) 


(b^ + ib - 1) e 


ib 


+ i7T I(-b) 


The integral containing the 6 function is given by 


in I(-b) = in 


(b^-ib - l)e (b^ + ib - l)e^*^ 


2i 


2i 
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Noting that the terms containing e 


ib 


add, one finally obtains 


r 1(g) d/3 

i (b + g) 


(b^ + ib-l)e^*^ + (b^-ib - l)e 


Collecting terms, the g integral under consideration is 


/ (3d^FW (-bT^+TTr)' -9 (-D+W [<b‘ + lb-l)e'' 


+ (b^ - ib - l)e 


■ib 


= j -^ + iTTb F'^(b) j 


Combining the fractions on the left hand side of the above equation, one 


obtains 2 J 


gdg F(g) 


Dividing by the factor of 2, complete agree- 


(g + b) (g - b) 

ment with equation (C-1) is obtained. The reason for including this alternate 
method of handling this type of integrals is to provide a means of checking 
some of the results and because sometimes this method is easier to apply. 

The next integral is given by 


J /3dg F(g) 


(b^ + g^) 
(b + g) 


r gdg F(g) r g^dg F(g) 
(b + j3) ^ (b + /3) 


Multiplying by (/3 - b) on the numerator and denominator and interchanging 
b and /3 on two of the terms, one obtains 


/ gdg F(g) 


( b^ + g^) 
(b + g) ■ 


r g^dg F(g) r g^dg F(g) 

•1 (b + g) (g - b) J (b+g) (g - b) • 
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The first of these integrals has already been evaluated; the other integral is 
evaluated as follows; 


/ 


-f 


/3^d0 F(/3) 


+ b) (/3 - b) 





\ 1 

/ ;S2 + i/3 - 




03 + b)(/3 - b) 13'^ 

\ 2i 

j 

t /7\ 

-p \ 

. f _y7^ 1 ^ 

' /3^d/3 FW) 1 / 

(^ + b) 0 - b) ^ \ 

M ^ 

. 2i 


-b 

+b 


J 7 ^ Li^ 


(i8 + b) ()3 - b) 


€-K) 


7T 


J F (z = -b + ee^^)dz 


+ f F (z = b + ee^^)dz 

6 


b^ F (b) 


(C-3) 


Combining terms, one obtains 


/ bdbF(b) f fidfi F(,3) (-^7^) 


= 2 J b^db F(b) + Y b F'^(b) j + 2 / bdb F(b)[^ b3 F'^(b) 


= I IT J bdb F(b) + 2 iir J b^db F(b) F^(b) | 


The integral corresponding to the above, but containing G(j3) instead of 
F(/3) , is obtained directly from the above results. 
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The next integrals of interest are given by 



/ bdb f /3d/3 (-g^) ?(/3) + 4Fa(b) Fa(^)] 


Multiplying by* (b - /3) in the numerator and denominator of the above 
expression and recombining as before, the above expression becomes 


J J (b + /3) (/3 - b) 


= 2 f bdb ?(b) J 


/3=^dgy(/3) 

(b + ^) (/3 - b) 


+ 8 f bdb Fg(b) f 


(b + /S) (|3 - b) ’ 


The /3 integrals are of the same form as the ones considered previously 
with only J (/3) and replacing F(/3) and G(j8). Following the steps 

leading to equation (C-1), one notes that since 3^(/3) and F^i^) have 
different coefficients than F(/3), only the residue at /3 = 0 changes value. 
This can be seen by considering these integrals as follows: 


r i8^d/3 5(g) 

J (b + /3) (^ - b) 



dff(i/3^ - 2/3^ - 3i0 + 3)e^^ 
2i/30 + b)()3 -b) 


= Lim 
e-^ 


J F(z = -b + ee^^) dz 

.0 



JT 10 

+ J F (z = ee )dz + J 
0 0 


10 

F(z = b + ee )dz 
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(C-4) 


Note that the sign change in the (37r/2b^) term is due to the corresponding 
sign change in the constant factors of the j3 polynomials in ? (/3) and F^OB) 
respectively. A similar integral occurs with exactly the same form as the 
above term, being considered with only ? (/3) replaced by Since 

3^2 (^) = G(/5) - 2 Fgifi) , which when expanded gives 



and noting that the residue at the origin contains the constant factor in the 
polynomial factor in )3, the result for this integral is 


oo 


/ 



i/3 + b) (/3 - b) 


(2) (37t) 


~3^2(b) 


The next group of /3 integrals are contained in the term 
/ bdb f (b2+^2^^ (F(b) F(^) +2 G(b) G(^)} 
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Expanding these factors and rearranging terms as before, the preceding 
expression becomes 


, 2 J r ^^dfflF(b)FW^2G(b>GOT| 

J J (b + /3) (/3 - b) 

+ 4/ tfdb f <3»d;s|F(b)FW4 2G(b)Gwl 

+ 2 f bdb f iF(b) Ftf) 4 2 (i . (b) , G(g) . ! .. 

(b + /3) (/3 - b) 


Analysis of the above terms shows that only the last term needs to be 
evaluated. Using the same procedure as before, 


T F(/3) ^ T 0^dl3U3^ + ig - i)e^^ 

J (/3 + b)(j3-b) 2i (/3 + b)(^-b) 



(C-5) 


The next integral appears in the term given by 
J bdb / (b2+/32) {5(b) +4F3(b) F3(/3)} 

_ 2 r b^db r {7(b) J(^) +4F3(b) F 3(^)1 

- " J ^ (b + ^)(^-b) 

+ 2 r bdb r ^(b) 703) + 4F3(b) F3(^)} 

J J ()8 + b)(/3-b) 


The first integral above has been evaluated previously. The second 
integral is obtained as follows; 
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f g^d/3 ?(/3) 
n: (/3 + t>) (/3 - b) 


= Lim 
e -*0 




ITT 

2 


b3y-(b) 


Similar operations are used to obtain 


j (/3 + b) (^ - b) 



(C-6) 


The next integral group is found in the term 


/ bdb / |b2/s2 ([F2(b) F2(/3) +2G(b) G(/?)] 




+ 24 Fg(b) F3(;S) 


])) 


= 2 / b’db / 2 I 


+ 2 


/ bab / 


■ 

(/? + b) (^ - b) 


X {f(b) f(^) + 12J4(b) J4(/3) + 24 F3(b) F3(^)| 
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Using the fact that 


F2(/3) = ^ 


f(/3) = 


W^ + 2 10 - 2)e^^ (/3^ - 2i0 - 2)e~^^ 


2i 


2i 

i/3 


()3^ + 4i/3® - 12/3^ + 24i)3 + 24)e 
2i 

(0* - 4ij8^ - 12 |8^ + 24ij3 + 24) 

2i J 


? 4 (/ 3 ) = 


(i0^ - 50^ - 12ij3 + 12)e^^ (i0^ - 50^ - 12i0 - 12)e~^^ 


2i 


2i 


one can write the results for the above integrals in the order in which they 
appear: 


2 Jb3db[^F2(b) b3 Ft(b)j + 2 G(b) b3 G'^(b)} j 
+ 2/bdb[f(b) bf^'ib)} 

+ 2 / bdb [l2 ?4(b) j - ^ b 54‘(b) j J 

+ 2 / bdb 1^24 Fg(b) { -^ + ^ b Fg (b) I 

Integration and Transformation of b- Integrals 

The b-integrals resulting from the 0 integrations are considered in 
this section. In general, two types of integrals are considered. In some 
cases the integral's principal value is obtained; whereas in other cases, the 
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integral is simply transformed into another integral expression. This 
behavior results from the fact that the evaluation of the b-integrals in closed 
form is not possible [22] . The most one can do is express quantities in 
terms of sine and cosine integrals defined by 


' I “ I -^Ob 

0 00 


The principal value of the first integral to be considered is obtained 
as follows: Taking 


oo 

r 

(b^H- 3ib - 3)e 

(b^ - 3 ib - 3)e 

A 

0 

2i 

2i 


and letting b -► -b in the second term, tiie above expression becomes 


r db (b^ + 3ib - 3)e^*^ 
J ab^ 2 i 


Using the contour indicated, this integral becomes 


oo o i b 

r db (b^ + 3 i b - 3) e 

J ab^ 2 i 


= Lim 
€-*0 


f F(z = ee^^)dz 
0 
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Using the L' Hospital rule in carrying out the limiting operations, the 
principal value of this integral is 


7 db (b^ + Sib - 3)e^^ / ir \ 

J "ab^ 2i ^ 4a/ ’ 


A similar procedure is used to evaluate the next integral, given by 


r db 

J 2ab^ 


(b^ + 2 i b^ - 5 b^ - 6 i b + 3) e 
2i 


2ib 


(b^ - 2ib^ - 5b^ + Gib + 3)e 
2i 


-2ib 


f db (b^ + 2ib^ - 5b^ - 6ib + 3)e^^^ 
J 2ab^ 2i 

— oo 



The next integral is given by 


au 

r db 

/ 2ab (a + b)^ 


(b^ + 2ib^ - 5b^ - 6ib+ 3)e 
2i 


2ib 


( b^ - 2ib^ - 5b^ + 6ib+ 3)e 
2i 


-2ib 
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Since the two parts of the above expression cannot be combined in the manner 
used before, one evaluates each of tihe above terms by considering the 
complex integral 



integrated around the indicated contour [23]. Since no poles are enclosed, 
this integral is given by 


^ F(z)dz 


L 7t/2 

Lim f F(z = b)dz+ f F(z Le^^)dz 
e -*0 € 6 


T I — ►oo 

i € 0 

+ J F(z = iy)dz+ J F(z = ee^^)dz 

iL 7 t/2 


0 


Making use of the convergence at infinity requirement the following is 
obtained: 


L 

Lim f F(z = b)dz 
e-*0 e 

Li-* 00 


Lim 

e-K) 

Ij— 


7t/ 2 . i L 

J F (z = ee^ )dz + J F(z = iy)dz 
0 i € 


The left-hand side of the above equation is just the principal value of the 
integral under consideration; hence. 
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/ 


2 ab(a 


_db I (b^ + 2ib 

a + b)^ I 


.3 cu2 


- 5b - 6ib + 3)e 


+2ib 


2i 


/37t\ i dy (y^ + 2y^ + 5y^ + 6y + 3)e 

\8^) a(iy) (a+ iy)^" 2i 


e-*0 i e = i y 

Lr-»-°o 


" (I?) \f 


dy 


4a(iy) (a + iy) 


r (y^ + 2y® + 5y^ + 6y+3)e 


—2 i b 

The term containing the e is evaluated in the same manner, but the 

corresponding complex integral ^ F(z)dz is evaluated, using the indicated 
contour: 



db 

2 ab(a + b)^ 


(b^ - 2ib^ - 5b^ + 6ib + 3)e 
2i 


( 



) 


-/ 


dy 


4(iy) (a - iy) 


j (y^ + 2 y® + 5 y^ + 6 y + 3) e' 


■2y 



Combining the above results, the integral being considered equals 



j° dy<(a-ly)»-(a4-^ly)^} ^ (^ 4 ^ 2yS+ 5y^+ 6yn- 3} . 

0 4a(iy) (a2+ y2) 
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Rearranging factors, one finally obtains 



(b^ - 2ib^ - 5b^+ 6ib+ 3)e~ 
2i 



/ {y^ + 2y3 + 5y2 + 6y + 3 |e'^^ 

0 { b ? + y^) 
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